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ABSTRACT 


Using analysis results from the POST trajectory optimization pro- 
gram, an adaptive guidance algorithm is developed to compensate for den- 
sity, aerodynamic and thrust perturbations during an atmospheric orbital 
plane change maneuver. The maneuver offers increased mission flexibili- 
ty along with potential fuel savings for future reentry vehicles. 
Although designed to guide a proposed NASA Entry Research Vehicle, the 
algorithm is sufficiently generic for a range of future entry vehicles. 
The plane change analysis provides insight suggesting a straight-forward 
algorithm based on an optimized nominal command profile. Bank angle, 
angle-of-attack, and engine thrust level, ignition and cutoff times are 
modulated to adjust the vehicle's trajectory to achieve the desired end- 
conditions. A performance evaluation of the scheme demonstrates a capa- 
bility to guide to within 0.05 degrees of the desired plane change and 
five nautical miles of the desired apogee altitude while maintaining 
heating constraints. The algorithm is tested under off-nominal condi- 
tions of ±30% density biases, mwo density profile models, ±15% aerodyna- 
mic uncertainty, and a 33% thrust loss and for various combinations of 
these conditions. Based on fuel comparisons with results of the opti- 
mization program, the guidance scheme offers a nearmoptimum solution 
without the complexity of real-time optimization. 
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CHAPTER 1 


INTRODUCTION 

A potential joint NASA/Air Force atmospheric reentry research flight 
experiment is being considered for the 1990 time frame. This flight 
experiment will help provide the technology base for development of 
future reusable ear th-to-orb i t transportation systems and Air Force 
aero/space vehicles. The experimental vehicle is envisioned to be car- 
ried aloft and deployed in orbit by the Space Shuttle. After executing 
a deorbit procedure, the vehicle will reenter the atmosphere to perform 
a variety of maneuvers and experiments and eventually glide to a landing 
site. 

One of the more interesting maneuvers is the synergetic orbital 
plane change. A synergetic plane change is a maneuver in which a change 
in the orbital inclination is accomplished through a combination of 
aerodynamic and propulsive forces rather than through propulsion alone 
(17). The guidance system presented here is designed to guide a winged 
vehicle such as the proposed NASA Entry Research Vehicle (ERV) during 
the atmospheric phase of a synergetic plane change mission. This mis- 
sion involves a deorbit maneuver and aerodynamic turn in the atmosphere 
combined with a powered cruise and exit back into orbit. The guidance 
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system takes the vehicle from entry interface (400,000 feet), through 
the atmospheric maneuvering turn and back to orbit. 


1 .1 THESIS CONTRIBUTION AND SCOPE 

This thesis present the problem analysis, development, and testing 
of a guidance algorithm for a synergetic plane change maneuver for 
NASA's ERV. The motivation for this effort was the design analysis 
requirements for the development of the guidance scheme. The major con- 
tribution is the new level of insight on how to efficiently perform this 
maneuver in the presence of density and vehicular perturbations. The 
guidance algorithm is a summation of this insight and satisfies NASA's 
requirement for a generic guidance scheme. The testing of the algorithm 
demonstrates the effectiveness of the scheme and indirectly adds support 
or proof to the analysis results. Hence the parts fit nicely together 
as a complete monograph with each part both resting on and supporting 
the other parts. 

The guidance algorithm is only one possibility of many guidance 
schemes that could and will, as demand requires, be developed to satisfy 
the mission. However, it completely satisfies the stated NASA require- 
ments at this conceptal phase of ERV development. In addition, the 
algorithm offers a basic scheme to modify and refine as the ERV program 
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progresses and a standard to measure future schemes against. In short, 
the guidance algorithm is a simple straight-forward application of the 
results from the synergetic plane change trajectory analysis. 

The algorithm testing was conducted on a tmree degree of freedom 
computer simulator developed at The Charles Stark Draper Laboratory 
(CSDL) for the NASA ERV contract. The algorithm testing is by no means 
exhaustive but for the time and money constraints, it does offer a good 
validation of the guidance algorithm. 

1 .2 Thesis Overview 

In the remaining six chapters, this thesis details the problem 
development and analysis along with guidance algorithm development and 
testing. Chapter 2 presents background information for the ERV and syn- 
ergetic plane change problem to include concept justification, vehicle 
description, an explanation of the baseline trajectories, and basic NASA 
guidelines and requirements. Chapter 3 contains a brief discussion of 
the initial problem development and formulation. Chapter 4 details the 
problem analysis with a computer optimization program called "POST". 
This analysis provides insight into the difficulty of the guidance prob- 
lem along with specific means of handling density, aerodynamic and pro- 
pulsive uncertainties. Chapter 5 presents our final iteration of the 
guidance scheme along with rationale leading to the scheme. 
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The evaluation criteria and test results for the algorithm are con- 
tained in Chapter 6. Chapter 7 summarizes the work covered in this the- 
sis and discusses what can be done to expand on the analysis insight and 
possible algorithm improvements. Appendix A provides a brief explana- 
tion of "POST" along with an example trajectory generation. A brief 
description and flowchart for the Draper Laboratory ERV simulator and 
guidance algorithm computer programs are presented in Appendix B. 


CHAPTER 2 


PROBLEM BACKGROUND 


2.1 Introduction 

During the past few years, there has been resurgency of national 
interest in vehicles that can maneuver in the atmosphere while returning 
from orbit and vehicles that can fly in the atmosphere at hypersonic 
speeds for sustained periods of time. This rebirth of interest in 
hypersonics has resulted both from the success of the Space Shuttle and 
from a developing awareness of the potential benefits that can be 
derived from vehicles with the capability to operate between the limits 
of existing aircraft and spacecraft. Examples of the types of maneuvers 
that are projected for this new class of vehicles are presented in Fig- 
ure 2-1. Most of the projected requirements are dictated by either a 
need for flexibility offered by a large landing footprint or by a 
requirement for specific entry maneuvers such as a synergetic plane 
change (4) . 

For vehicles that operate in space or in the upper reaches of the 
atmosphere, an orbital plane change offers significant flexibility for 
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Earth overflight and space operations. Specifically, the studies of the 
late 1960's indicated that the same type of high lift-to-drag ratio 
(L/D) vehicle required for the synergetic plane change maneuver could 
also perform Earth-return missions with large cross range capability, 
could perform evasive maneuvers, and could have considerable military as 
well as civilian uses (17). This new class of vehicles will dictate new 
requirements for the vehicle's thermal protection systems (TPS) and 
thermo-structures design. The plane change maneuver will require a TPS 
designed for a high heat-rate while the large landing foot print will 
require a high heat-rate and heat-load capability (6). The Space Shut- 
tle provides a good starting point but to achieve the potential fuel 
efficiency and large landing foot print will require an improvement in 
TPS capabi 1 i ty . 

Addressing the technology deficiencies for these future systems will 
provide a significant challenge for the aerospace community. The syner- 
getic or atmospheric plane change will be one of the more difficult mis- 
sions for future entry vehicles. Analysis to date have shown that 
performance of the maneuver will require sustained flight at or near 
entry mach numbers at altitudes over 200,000 feet. This is the flight 
regime where the predicted aerodynamic trim of the Space Shuttle orbiter 
was in error by as much as nine degrees body-flap deflection (4). Wher- 
eas the shuttle currently flies a tightly tailored trajectory, the 
flight envelope for the ERV will be broadened to allow the vehicle to 
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maneuver during entry. The results from Space Shuttle flights presented 
in Reference 4 have shown density variations as large as 60 percent in 
the altitude regime where the synergetic plane change maneuver occurs. 
This large density variation is random, can not be predicted and must, 
therefore, be treated as an uncertainty in the estimated aerodynamics 
and aeroheating used in the vehicle design. 

With the ERV, NASA is attempting to address these issues and develop 
a building block approach to handling the difficulties. The ERV is 
designed to demonstrate atmospheric maneuvers that are beyond the capa- 
bility of existing spacecraft and aircraft. Presented in Figure 2-2 is 
a typical altitude time profile for an entry research vehicle flight 
where first a synergetic turn is performed, and then the vehicle per- 
forms a normal entry (4). As shown in the figure, the entry vehicle is 
released from the shuttle, performs a deorbit burn, and then descends 
into the atmosphere to an altitude of approximately 220,000 feet. There 
the engines are started, and the vehicle cruises at a mach number of 25 
until its exit point. It then uses a combination of angle-of-attack and 
throttle to reboost into orbit while completing the plane change. After 
remaining on orbit for several revolutions, the vehicle deorbits and 
performs an entry and landing. 

In order to accommodate the large density variations and vehicle 
aerodynamic uncertainties referred to earlier, the new vehicle will be 
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Figure 2-1. Candidate missions for the entry research vehicle [4], 
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Figure 2-2. Typical entry profile for an ERV synergetic plane change 
followed by reentry [4]. 
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required to have a sophisticated adaptive guidance and navigation sys- 
tem. This system must be able to sense the atmospheric conditions and 
provide necessary control inputs in real-time to enable the vehicle to 
accurately achieve its end-conditions and remain within the vehicle's 
flight boundaries. The best previous example is the Space Shuttle. 
However, the shuttle worked around these uncertainties by designing into 
the guidance algorithm large margins for error. This resulted in a 
reentry trajectory with less concern for maneuvering efficiency than is 
possible with the ERV and often the solution to handling the uncertain- 
ties (e.g., roll reversals) resulted in other problems like attitude 
instability of the vehicle (8). The ERV is designed to demonstrate the 
technologies required to perform a synergetic plane change mission and 
other maneuvering entries. Thus, an acceptable guidance system must be 
able to efficiently compensate for the discussed uncertainties while 
safely executing the maneuver. 


The strongest argument for a synergetic plane change is the opportu- 
nity for a significant reduction in the amount of fuel required for an 
orbital plane change maneuver. Many studies have shown that for a large 
plane change in low-earth orbit, a synergetic plane change will use less 
fuel than an all-propulsive maneuver. The initial motivation for this 
work was the recognition that orbital plane changes required very large 
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character i si tc velocities (aV) . As pointed out by Caudra and Arthur 
(18), the AV required for a single-impulse 60 degree plane change is as 
large as that required to place a satelite in low earth orbit in the 
first place. The fuel required for the 12,000 pound ERV to perform a 
20-degree orbital inclination change is shown in Figure 2-3, and the 
corresponding velocity change requirement is shown in Figure 2-4. 
Results are shown for an in-orbit, a 1 1 -propul s ive maneuver and for two 
synergetic plane change maneuvers, one constrained to a maximum convec- 
tive heat-rate referenced to a one-foot radius sphere of 80 BTU/ft 2 -sec 
and one constrained to a maximum heat-rate of 125 BTU/ft 2 -sec. For the 
ERV to perform an all-propulsive plane change requires 7100 pounds (lbs) 
of propellant. The same plane change using an aerodynamic turn limited 
by heating to 80 BTU/ft 2 -sec requires 6400 pounds of propellant. If the 
ERV can withstand a heating rate of 125 BTU/ft 2 -sec, the fuel savings 
are even greater, and only 5500 lbs of fuel would be required (6). It 
should be kept in mind however, that significant fuel savings are only 
realized for a vehicle that already requires the added weight for TPS 

and wings. Or as Haslau concluded " the results indicate that the 

synergetic made can be more desirable than in-orbit propulsion only if 
the aerodynamic performance must be provided for other reasons..." (20). 

Further analysis done at NASA Langley Research Center (LaRC) , deter- 
mined that there is an optimum heat-rate depending on the amount of 
desired plane change. Figure 2-5 provies a curve of this heat-rate 
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Figure 2-3. Comparison of fuel requirements for all -propul si ve maneuvers. 

and synergetic maneuvers to produce an orbital inclination 
change of the ERV [4]. 



Ai, deg 


Figure 2-4. Comparison of Delta V requirements for al 1 -propul si ve 

maneuvers and synergetic maneuvers to produce an orbital 
inclination change of the ERV [4], 
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verses inclination relationship for the ERV for angles through 20 
degrees. Specifically, the results indicated that 125 BTU/ft 2 -sec was 
nearly the optimum heat-rate for a 20 degree plane change and a higher 
heat-rate of say 135 BTU/ft 2 -sec will actually require more fuel (6). 
Thus, the maximum allowable heat-rate will significantly influence the 
vehicle's fuel efficiency and, as it will be shown later, the vehicle's 
trajectory and response to perturbations. 

2.3 Vehicle Description 

The candidate Entry Research Vehicle shown in Figure 2-6 is 25 feet 
long with a wing span of 13.9 feet. It is designed to occupy half of 
the Shuttle cargo bay which is 60 feet long and 15 feet in diameter. 
The vehicle has a distinct wing-body design with a blended wing-body 
interface to increase the overall hypersonic performance of the vehicle. 
The wing has a reference area of 177.4 square feet. Conventional ele- 
vons are used for pitch and roll control, and tip-fin controllers pro- 
vide the yaw stability augmentation and control. A body flap assists in 
trimming the vehicle hypersonical ly. Reaction control system jets are 
used to control the vehicle in the upper atmosphere. The jets are 
placed in the nose of the vehicle and in the tip fins. Figure 2-7 is a 
three-view drawing of the ERV with the basic vehicle dimensions and 
relationship with the Shuttle payload bay envelope. The fineness ratio 
of the fuselage is approximately six, and its volume is 300 cubic feet. 
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Figure 2-6. Artist's concept of the ERV [5]. 


Figure 2-8 shows an inboard profile of the ERV. The payload bay housing 
many of the experiments is located in the nose of the vehicle where 
atmospheric sensing can be performed with minimal engine combustion by- 
product contamination. The avionics are located behind the payload com- 
partment. To minimize the center-of -gravi ty shift as propellant is 
depleted, the fuel and oxidizer tanks are positioned around the center 
of gravity which is located at 67 percent of the body length. The elec- 
trical system and batteries are positioned between the fuel tanks, 
approximately at the center of gravity. Propulsion is provided by three 
Marquardt R-40-B rocket motors, each providing 1100 pounds of thrust 
with an ISP of 295 seconds. The engines are aligned with the vehicle's 
longitudinal axis and thrust through the vehicle's center of gravity. 

The vehicle's launch weight from the Shuttle is estimated to be 
11,430 pounds. 6000 pounds of propellant is carried to provide the 
nearly 7000 f eet-per-second 4V required for a 20 degree synergetic plane 
change maneuver with the Marquardt engines. The leading-edge sweep 
angle is 72.5 degrees. An NACA four-digit series symmetrical airfoil 
mection is used for the wing and its hypersonic maximum L/D is approxi- 
mately two. Complete coefficient of lift (CL) and drag (CD) data is 
provided in Appendix A. The thickness ratio was varied to maintain a 
ha If- inch radius over the entire leading edge of the wing. The geomet- 
ric shape of the wing was designed with the intention of reducing the 
aerothermodynami c heating on entry. A heat pipe, in addition to new 
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thermal protection systems and materials, will be used to accommodate 
excessive temperatures (6) . 

2.4 Basel i ne Trajectory 

2.4.1 Background and Strategy 

The vehicle analysis branch at LaRC defined a baseline or nominal 
synergetic plane change trajectory. This baseline trajectory was devel- 
oped by combining the best ideas and concepts that have emerged since 
1961 when London's paper first convincingly demonstrated a significant 
performance gain with a synergetic plane change. A brief listing of 
these basic concepts are shown below as an aid to understanding the 
vehicle's design and baseline trajectory. A more detailed discussion of 
these concepts are in reference 17. 

1. For plane changes of less than 15 degrees, the a 1 1 -propul s i ve 
maneuver is generally more efficient. 

2. An L/D of at least two is required to offer a significant advan- 
tage over the all-propulsive plane change, and it is desirable to 
maximize vehicle L/D. 
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3. Turns made at maximum CL (i.e., high angle-of-attack) are quicker 
and hence more fuel efficient than the more gradual maximum L/D 
turns. This improved efficiency is primarily a result of orbital 
mechanics. A plane change made at the intersection of two orbits 
(ie. t the node) produces all inclination change whereas a turn at 
the orbit apex (90 degree from node) produces no inclination 
change, only a node shift. Hence, for maximum inclination change 
and minimum node shift, the turn should be centered over the node 
in the shortest duration possible. 

4. A large reduction in total heat load can be achieved by carrying 
out a short duration, high angle-of-attack maneuver rather than 
the slower maximum L/D turn. 

5. A thrusting turn (aerocruise) offers significant advantages over 
the gliding turn (aeroglide) when the desired plane change is 
more than 10 degrees. This relationship is especilly true when 
lower wall temperature (heat-rate) limits are imposed. 

6. Aerocruise vehicles use continuous thrusting sufficient to bal- 
ance aerodynamic drag and maintain an altitude high enough to 
allow optimum bank angles throughout the turn without exceeding 
heat-rate 1 imi ts . 
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7. For the high angl e-of-attack aerocruise vehicle, assuming thrust 
along the longitudinal axis, a significant portion of the plane 
change results from out-of-plane thrust due to the high angle-of- 
attack and bank attitude of the vehicle during the turn. 

8. For the quick, high ang 1 e-of-attack turn, the vehicle is banking 
throughout the maneuver, and the exit should be timed such that 
the plane change is accomplished just as the vehicle exits the 
atmosphere to minimize fuel. 

9. For each amount of plane change, there is a heat-rate that offers 
optimum fuel efficiency. If this heat-rate can be achieved, the 
resulting trajectory is more a combination aerogl ide-aerocrui se. 
A part of the turn is achieved on the gliding entry and the 
remaining plane change is achieved during the powered exit with 
almost no level-flight phase. 

2.4.2 Trajectory Development and Understanding 

Using these basic concepts, R.W. Powell and J.C. Naftel of NASA 
Langley Research Center developed a baseline synergetic plane change 
trajectory for a reference heating rate of both 80 and 125 BTU/ft 2 -sec. 
These trajectories were generated on their three-degree-of-f reedom Pro- 
gram to Optimize Simulated Trajectories (POST). References 1,2, and 3 
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offer a detailed description of this program while Appendix A provides a 
brief description of POST along with an example of trajectory generation 
for the synergetic plane change maneuver. 

The strategy for both baseline plane change trajectories are simi- 
lar. First a deorbit burn is performed. Once the vehicles enters the 
atmosphere and the heating rate reaches the desired cruise value, the 
engines are restarted and throttled to maintain a nearly constant vehi- 
cle velocity. At the end of this cruising phase the engines are com- 
manded to maximum thrust, and the vehicle exits the atmosphere. Bank 
angle and ang 1 e-of-attack are used during the entire atmospheric phase 
so that the required inclination change occurs during the entry, cruise 
and exit legs. Figures 2-9 thru 14 show profile plots for altitude, 
angle-of-attack, bank angle, engine throttling, heat-rate (QDOT) , and 
heat load histories for both the 80 and 125 BTU/ft 2 -sec 20 degree syner- 
getic plane change trajectories (6). Figure 2-15 shows a typical POST- 
generated altitude plot and the generic phase names that will be used 
during the plane change trajectory discussion. These profiles are pre- 
sented here as basic reference material that should be helpful during 
the discussion in the remainder of this chapter. Appendix A provides 
more details into the generation of these trajectories while Chapter 4 
will cover these trajectories and their plots in detail. 
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Figure 2-9. Altitude profiles for synergetic plane chanqe usinq 
the ERV [6], 
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Figure 2-10. Angle-of-attack profiles for synergetic plane change 
using the ERV [6]. 
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Figure 2-11. Bank angle profiles for synergetic plane change using 
the ERV [6], 
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Figure 2-12. Engine throttling profiles for synergetic plane change 
using the ERV (O-off, 1 maximum) [6]. 
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Figure 2-13. Heat-rate profiles for synergetic plane change using 
the ERV [6], 
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Figure 2-14. Heat load profile for synergetic plane change using 
the ERV [6]. 
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Figure 2-15. Post synergetic plane change altitude profile and generic 
phase names. 


Langley's objective was to obtain the most fuel efficient plane 
change trajectory (within the adopted strategy) by allowing POST to find 
optimum values for carefully selected guidance parameters (independent 
variables). The 12 independent variables used for the baseline trajec- 
tories are listed below with a brief description of how POST uses the 
parameter . 

1. Geocentric Latitude : The latitude at which a deorbit burn is 

started . 

2. Perigee Altitude : The conically predicted orbital parameter 

which determines when the deorbit burn is stopped. This determines 
the AV of the deorbit burn and thus the position and velocity of the 
vehicle at entry interface. 

3. Entry Ang 1 e-of-Attack : Angle-of-attack from atmospheric entry 

to level -off. 

4. Entry Bank Angle : The bank angle from atmospheric entry to lev- 

el-off . 

5. Cruise Ang 1 e-of-Attack : Angle-of-attack for the cruise phase of 

f 1 ight. 
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6 . 


Cruise Acceleration : 


The acceleration maintained during the 


cruise phase of flight. 

7. Exit Inc 1 i nat i on : The inclination angle where the exit is 

started. 

8. Exit Bank Angle ; The bank angle during the powered part of the 
exit. 

9. Exit Anqle-of-Attack : The ang 1 e-of -attack during the powered 

part of the exit. 

10. Apogee Altitude: 
which determines when 
mines when the vehicle 
al ti tude. 

11. Unpowered Exit Angle-of-Attack ; The angl e-of-attack after the 
engines are cutoff during the exit phase. 

12. Unpowered Exit Bank Ang le: The bank angle after the engines are 

cutoff during the exit phase. 

The constraints and targeting parameters include: 


The conically predicted orbital parameter 
the rocket engines are cutoff. This deter- 
can fly unpowered back to the desired orbital 
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1 . Maximum heat-rate 


2. Desired inclination angle 

3. Desired final orbital apogee altitude 

4. Maximum throttle setting 

It should be noted that total heat-load is not considered a constraint 
for the relatively short plane change maneuver. 

The optimizing parameter is final vehicle weight in the desired 
final circular orbit. Thus, POST attempts to minimize the amount of 
fuel used from the start of the deorbit burn until the ERV is back in 
the desired circular orbit. During this time, the 12 guidance parame- 
ters (free parameters) are adjusted to achieve the targeting parameters 
of inclination angle, apogee altitude, and maximum throttle setting. In 
addition, as a means of controlling heat-rate and achieving the cruise 
acceleration, POST can respectively vary bank angle and throttle setting 
during the level-off and cruise phase. 

To simplify discussion here and throughout the remainder of this 
thesis, an abbreviated notation will be used to represent a heat-rate of 
80 or 125 BTU/ft 2 -sec. Typically, the numbers 80 or 125 will be used 
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without units as stand-alone adjectives to describe trajectories or 
cases, or as objects (e.g., heat-rate of 80 or heating rate of 125). 
The more complete notation will be used where any uncertainty of meaning 
might exist or occassionally as a reminder. 

To provide additional understanding of the baseline trajectories, a 
more detailed explanation of the previously presented trajectory plots 
will be made. Figure 2-9 provides an altitude verses time profile 
starting and ending at 400,000 feet (maximum extent of atmosphere). For 
the synergetic plane change with a design maximum convective heat-rate 
of 80, the ERV initiates level-off approaching the design heat-rate, and 
the engines are ignited. However, for the case with the design maximum 
heat-rate of 125 BTU/ft 2 -sec, the ERV phugoids at a heat-rate of approx- 
imately 100 BTU/ft 2 -sec and then continues to glide (engines off) until 
the design maximum heat-rate is reached. Once the design maximum heat- 
rate is reached, the level-off is started and the engines are ignited. 
For the 125 case, the exit maneuver begins almost immediately. For the 
80 case the cruise phase continues significantly longer before throttl- 
ing to 100% and boosting back into orbit. These differences are the 
result of a higher heat-rate allowing the vehicle to fly deeper into the 
atmosphere, thus generating higher lift to turn the vehicle quicker. As 
stated earlier, the quicker the maneuver can be performed, allowing the 
vehicle to remain close to the orbital node, the more efficient the 
maneuver . 
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Figure 2-10 presents angi e-of-attack profiles that are reasonably 
similar for the 80 and 125 cases with length of the cruise phase making 
the biggest difference. The bank profiles in Figure 2-11 are again sim- 
ilar for the two cases with the spike just prior to 1000 seconds due to 
the ERV's heat-rate controller adjusting the bank angle to maintain the 
designed maximum heat-rate during level-off. Figure 2-12, the engine 
throttling profiles, displays the high percentage of maximum throttle 
required for level-off and cruise for the 125 case as compared to the 
relatively low throttle setting required for the 80 case. In Figure 
2-13, the "double-peaked" heat-rate (QD0T) curve for the 125 case stands 
out quite clearly from the "Table-top" curve of the 80 case. Finally, 
the heat-load profiles in Figure 2-14 demonstrate the significantly 
higher total heat- load for the 80 case. 

2.5 Sjuidel.ines._and R equirements 

NASA contracted with the Charles Stark Draper Laboratory to devolop 
a guidance algorithm to satisfy the broad requirements for a future 
entry vehicle. As stated earlier, the new class of vehicles will 
require a sophisticated adaptive guidance and navigation system. This 
system must be able to sense the actual conditions and provide necessary 
control inputs in real-time to enable the vehicle to accurately achieve 
it's end-conditions and remain within the vehicle's flight boundaries. 
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A primary objective in the development of the guidance system for 
the ERV was that it be flexible. Many past guidance algorithms have 
lacked flexibility because of severe computational constraints. With 
the availability of more advanced flight computers, the intent is to 
devise a guidance law offering more flexibility even if it is more com- 
putationally intensive. 

In addition, the guidance system should not require exhaustive pre- 
mission planning. It should explicitly calculate on-board the desired 
trajectory or guidance commands instead of relying on elaborate tables 
which are highly mission dependent. Overall, this desire for greater 
flexibility and independence is in response to the need for increased 
operational performance to go along with the improved capabilities of 
future space vehicles. 

Other general guidelines include: 

1. Explicit control of vehicle heat-rate 

2. Maximum ability to handle dispersions 

3. High accuracy 
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4. Investigate the use of both angl e-of-attack and bank angle as 
control variables. 

5. Control capability for powered and unpowered maneuvering during 
the atmospheric flight phase. 

6. Minimize fuel usage 

7. Usable in vehicles with low roll rates 

8. Transportable between vehicles. 

Specifically, the plane change guidance system should be able 
to: 

1. Achieve an inclination change of 20 ± 0.1 degrees and a final 
apogee altitude of 160 ± 5 nautical miles (nm) . 

2. Demonstrate acceptable performance for a heat-rate constraint of 
either 80 or 125 BTU/ft 2 -sec. 

3. Demonstrate acceptable performance with atmospheric dispersions 
(±30% density biases along with density shear models), aerodyna- 
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mic uncertainties (±15% for all CL and CD combinations), and pro- 
pulsion performance loss of one of its three engines. 

One caveat should be pointed out concerning the requirements. The 
longitude of the ascending node (LAN) is not specifically required to be 
controlled. It is important to remember that an orbital plane is usual- 
ly defined in terms of inclination and longitude of the ascending node. 
This defines a single orbital plane with respect to the equator and the 
vernal equinox. Thus, a change from mne orbital plane to another 
orbital plane will, in general, cause a change in both of these parame- 
ters. However, for the purpose of this thesis, only an inclination 
change will be controlled. Longitude of the ascending node will be mon- 
itored to determine the magnitude of any shift. 
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CHAPTER 3 


PROBLEM DEVELOPMENT 


3.1 Introduction 

The objective is to develop a guidance system to take the ERV 
through the atmospheric phase of a synergetic plane change mission. The 
powered flight guidance for accomplishing the deorbit and recirculariza- 
tion burns is well developed and documented. Starting with a well-de- 
fined reference trajectory, a. specific set of end-condition 
requirements, one limiting constraint (maximum heat rate), and a large 
number of control parameters, the problem of developing a guidance 
scheme seemed rather straight forward. However, it was not immediately 
clear as to which combination of control variables provided the most 
efficient solution. Some type of analysis was done for each combination 
considered. However, some of the decisions made about the nominal guid- 
ance scheme were based on heuristic arguments rather than precise numer- 
ical comparisons. Naturally, where possible, both means were used to 
arrive at conclusions, however, numerical results alone can be mislead- 
ing if they are not fully understood. Where appropriate, information 
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from other sources was used and referenced, rather than data generated 
specifically for this study. 

Still, with the nominal baseline trajectories that are presented in 
Chapter 2, the problem is quite simply: how to develop an algorithm to 
use the vehicle's control variables to achieve a 20 degree plane change 
and a 160 nm apogee altitude in the presence of dispersions? Thus, it 
was a matter of choosing a basic guidance concept and a pairing of 
available control variables with the conditions to be controlled. 

3.2 Guidance Cha ract eristics and, Concepts 

Guidance systems span a spectrum of adaptability or flexibility from 
very rigid like an open- loop guidance system to a highly flexible system 
which continually generates an optimum trajectory from the current posi- 
tion to the desired end-conditions based on the latest knowledge. As 
stated in the previous chapter, the task was to develop a guidance sys- 
tem that is as flexible as possible within the constraints of the 
expected flight computer capability. 

In addition, NASA desired an explicit type of guidance system not 
based on a reference trajectory requiring a great deal of pre-mission 
planning and loading of data. In terms of explicit guidance concepts, 
three general types of systems are available: 
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1 . Perturbation 


2. Predictor/Corrector 

3. Real-time Optimization 

Perturbation guidance generally relies on the explicit solution of 
simplified equations of motion. Typical perturbation guidance trajecto- 
ries are equilibrium glide, constant drag, and constant heating rate. 
Often a relationship between the end-condition and the type of trajecto- 
ry can be determined. For perturbation guidance, guidance commands are 
composed of two parts — nominal and perturbed. In a broader sense, the 
nominal could be a nominal reference command or a nominal reference tra- 
jectory profile. In either case, a nominal is calculated each guidance 
cycle to achieve mission objectives. The perturbed part is used as 
feedback to bring the vehicle back to the nominal. This has been used 
by the Johnson Space Center for Shuttle entry and various aerobraking 
schemes to force a vehicle to fly a drag acceleration profile to meet 
mission objectives. The advantages for this type of scheme are: 

1. Lower sensitivity to initial conditions than an open- loop scheme 

2. Guidance parameters are calculated explicitly 
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3. Low computational load 

4. Flight experience on Shuttle 
While the disadvantages include: 

1. No direct control on the end conditions 

2. A lack of mission flexibility 

3. Scheme is usually very vehicle specific 

A predictor/corrector scheme on the other hand, generally works with 
more accurate equations of motion. The guidance continuously computes a 
new trajectory to go from the present position to the end-conditions. 
The sensitivity of the end-conditions to changes in the control vari- 
ables is determined by multiple numerical i ntergrat ions of the trajecto- 
ry. Once the sensitivity is determined, the current control profile can 
be corrected to satisfy the end-conditions. Advantages for this type of 
scheme are: 

1. Sensed dispersions can be easily incorporated into the predictor 

to improve the performance during the remaining part the mission. 
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2. The scheme is less vehicle specific and thus can be transferred 
between vehicles more easily because changes in the vehicle model or 
aerodynamic data base do not change the guidance scheme. 

3. The predictor can use a conservative vehicle model to allow 
additional margin for dispersions. 

4. Very flexible in terms of using different control variables 
While the disadvantages include: 

1. A higher computational load 

2. Possible convergence problems 

3. Not especially useful for control of boundary conditions (like 
heating rate) that are unrelated to end-conditions 

4. More difficult to verify and validate 

Real-time optimization provides the most desirable guidance sol- 
ution. At each guidance cycle a trajectory is calculated to hit the 
end-conditions and minimize a cost function such as fuel usage. The 
advantages include extreme flexibility and dispersion compensation capa- 
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bility while minimizing the costs. The disadvantages are just a magni- 
fication of the predictor/corrector disadvantages of high computational 
load along with reliability and convergence concerns. 

In the final analysis, a hybrid offered the best hope of satisfying 
all of NASA's guidelines. A perturbation guidance scheme is used to 
maintain the heat-rate boundary while a predictor/corrector is used to 
achieve the end-conditions. As a means of achieving at least a nearly 
optimum solution, insight gained from our trajectory analysis that will 
be covered in Chapter 4 is used to develop techniques which allow effi- 
cient handling of perturbed situations. 

3.3 C ontrol Va ri a bles and Pairings 

As a final step of breaking down the guidance problem, the degrees 
of freedom and the constraints or end conditions that must be satisfied 
with these control variables should be discussed. The objectives for 
the plane change are: 

1. Final orbital inclination angle control 

2. Final apogee altitude control 

3. Heat-rate constraint 
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4. Minimization of fuel usage 

While the parameters available to handle these conditions are: 

1. Angl e-of-attack 

2. Bank angle 

3. Engine thrust level (adjustable throttle setting) 

4. Engine start time 

5. Engine cutoff time 

Simple accounting points out the fact that there are more control vari- 
ables than trajectory objectives. Our approach to the problem was to 
explore the use of all parameters but use only what was necessary. This 
provides a straight forward scheme that can be more easily modified for 
later developments or other vehicles. 

Based on several sensitivity studies conducted by R. Richards of 
CSDL (19) , the effectiveness of ang I e-of-attack and bank angle was 
determined to be very similar during the last two-thirds of the plane 
change trajectory in controlling inclination. In addition, either con- 
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trol variable was shown to be capable of correcting inclination errors 
of up to 30% during the exit phase. The same group of studies indicated 
that final apogee altitude could be easily control with engine cutoff 
time. It is important to remember however, that because of the rela- 
tively high angles-of-attack during the powered flight phases, a signif- 
icant amount of the turn and hence inclination change is directly 
attributable to engine thrust. Therefore, engine cutoff time displays a 
large coupling effect on final inclination. In-spite of this effect, 
the pairing of bank or ang 1 e-of -attack with inclination change and 
engine cutoff with apogee altitude offered a simple, reasonable approach 
to the problem. 

For control of heating rate and fuel usage no simple study offered a 
clear indication of what control pairings offered the best means of con- 
trol. However, based on the author's experience with POST, it seemed 
possible to gain the additional insight into the problem through this 
avenue. 
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CHAPTER 4 


TRAJECTORY ANALYSIS 


4.1 Introduction 

Although the concept of performing an orbital plane change manuever 
by using aerodynamic forces has been seriously considered and worked on 
since the early 1960's(17), this thesis offers a new level of knowledge 
for the optimum control of the maneuver in the presence of density and 
vehicular perturbations. 

As discussed in Chapter 3, it became quickly apparent that some 
method was needed to provide insight into which of the various options 
was the preferred way to guide the ERV. Having some experience with 
POST, it was felt that at least a certain amount of that insight could 
be gained by simply perturbing the nominal conditions (POST initial con- 
ditions) and studying the initial and the optimized trajectory generated 
by POST. If nothing else, we could obtain the scope of the problem 
(i.e., the net effect of the perturbations on the end-conditions) for 
the required density, aerodynamic, and thrust perturbation/uncertainty 
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criteria. Furthermore, analysis of the optimized solutions generated by 
POST could help determine a generalized method for correcting the per- 
turbation errors in an optimum manner without the complexity of contin- 
uously determining the optimum solution. That is, a means to optimize 
the guidance process without using real-time optimization. As it turned 
out, not only did the work with POST offer a simple means to size the 
problem, but it offered enough insight to logically and directly develop 
a guidance law which easily met all the required conditions stated in 
Chapter 2. 

Overall, the basic plane change strategy and the analysis results 
obtained with POST appear to be generally applicable and sound. The 
performance evaluation results of Chapter 6 strongly validate most of 
the conclusions reached during our analysis. Still, it is important to 
remember that our trajectory analysis is based on POST and the nominal 
plane change trajectories that Langley developed with POST. Thus, any 
unique properties of how POST optimized or any changes or errors in the 
basic plane change strategy could effect the results and conclusions 
contained herein. 

The remaining part of this chapter will discuss the methods and 
results of our trajectory analysis with POST. After discussing the 
methods, overall analysis results will be covered. This will be fol- 
lowed by a detailed discussion of the scope of the problem and the 
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optimization/compensation results for density perturbations aerodynamic 
uncertainties, single-engine failure and various combinations of these 
conditions on the nominal trajectory. 

4.2 Method of Anal vs i s 

4.2.1 Q,bl.s&ti-Y.S 

As previously discussed, our objective with POST was to obtain the 
type and size of trajectory errors that were caused by changes in densi- 
ty, aerodynamic characteristics and available thrust. In addition, we 
were seeking ways to optimize or compensate for these changes. To 
achieve this objective, we started with the nominal guidance parameters 
(independent variables) but off-nominal conditions and observed the ini- 
tial and final trajectory generated by POST during its optimization pro- 
cess . 

As an aid to understanding this method, a brief explanation of how 
POST operates is in order. Using POST as an optimum trajectory genera- 
tor, it goes through an iteration process to determine the best sol- 
ution. Starting with the initial conditions and guidance parameters it 
generates an initial trajectory based solely on those inputs. POST then 
tries to improve the trajectory by changing the independent variables 
based on numerically determined partial derivatives. On each successive 
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iteration, it repeats the process of changing the guidance parameters 
based on the partial derivatives and compares the resulting end-condi- 
tions with the desired conditions. This process continues until the 
desired conditions are obtained. At the same time, POST is trying to 
minimize or maximize a cost function. Thus, the final trajectory will 
not only achieve the desired end-conditions but do it i n an optimum man- 
ner . 


4.2.2 Specific Approach 

All the trajectory analysis discussed in this chapter was performed 
with the POST computer program. The program was started with the ERV in 
a 160 nm circular orbit. The specific trajectory objective was to 
achieve a 20 degree synergetic inclination change using a maximum heat- 
rate of either 80 or 125 BTU/ft 2 -sec and return to a 160 nm circular 
orbit with the maximum vehicle weight (i.e., minimum amount of fuel con- 
sumed) . 

Starting with the 12 optimum independent variables for the nominal/ 
baseline trajectories (see Chapter 2 and Appendix A for additional 
information), various amounts and combinations of density and aerodyna- 
mic perturbations and/or single engine failure were added to the nominal 
conditions. During the first cycle through (first iteration of the 
optimization process), POST would use the nominal independent variables 
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to generate a trajectory. The end-condition results of this initial 
trajectory were the "qualified" open-loop results (not a true open-loop 
as will be discussed in the next subsection) caused by the perturba- 
tions. These results provide an indication of the net effect of the 
perturbation on the trajectory. 

After the first iteration, the independent variables are adjusted to 
find an optimum trajectory for the off-nominal conditions. Analyzing 
the manner in which the independent variables were adjusted for various 
perturbations, certain trends stood out that offered insight into how 
the ERV could best be guided. 

The final step in our analysis process was to systematically "fix" 
or "freeze" some of the guidance parameters (e.g., the four entry vari- 
ables) at their nominal value. By fixing the entry variables we were 
able to gain insight into how to optimize under what we felt were "real- 
world" constraints. In other words, you should not be able to adjust 
your entry parameters before you have actually had time to sense that 
entry conditions are not nominal. Thus, insight was gained on how to 
optimize various sections of the trajectory in addition to the trajecto- 
ry as a whole. 

4.2.3 Limitations 
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Before discussing the results of our trajectory analysis, three lim- 
itations or qualifiers on POST should be covered. The first qualifier 
concerns the 12 free parameters (independent variables). The analysis 
discussed in this chapter used the same free parameters that were used 
in developing the baseline trajectories. Therefore, the plane change 
trajectories discussed are really only optimum to the degree that POST 
was allowed to adjust trajectory parameters. Still, based on the know- 
ledge of how Langley developed the baseline trajectories and the work 
done for this thesis, the baseline trajectories do represent a very good 
approximation of the "true" optimum solution. 

The second limitation deals with the manner in which POST optimized 
for off-nominal conditions. As POST iterates to an optimum trajectory, 
it has apriori knowledge of perturbations. This is in contrast to a 
real-time guidance system that can only estimate what future perturba- 
tions will be. Therefore, not all methods used by POST to optimize are 
directly applicable to a real-world guidance algorithm. The process of 
"freezing" the entry parameters as discussed in the last subsection, was 
one approach used to overcome this limitation. 

The third qualifier deals with the way the baseline plane change 
trajectories were defined and how POST handles certain constraints. 
Specifically, the strategy for the baseline trajectories was to maintain 
a maximum design heat-rate and an optimum vehicle acceleration during 
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the cruise phase. It should be pointed out that the results of numerous 
POST runs indicated that the optimum vehicle acceleration in all practi- 
cality does not vary and is almost zero (0.0036 ft/sec 2 ). Thus, a near 
constant velocity is maintained which (as it will be shown later) 
implies a nearly constant altitude to maintain a constant heat-rate. To 
handle these constraints, POST is allowed to vary bank angle to control 
heat-rate and throttle setting to maintain vehicle velocity during the 
cruise phase. The heat-rate and constant velocity (actually a near zero 
acceleration) constraints are necessary to maintain the same type of 
trajectory from case to case, but they must be considered when inter- 
preting the trajectory analysis results. Specifically, it is important 
to realize that POST will always be able to adjust bank angle and throt- 
tle setting during the cruise phase. Therefore, as suggested earlier, 
the initial trajectory generate by POST during its optimization process 
is not really a true open-loop trajectory. Still, even with these limi- 
tations, this initial trajectory offers some valuable insight on the 
effects of the perturbations, and the results will be labeled as "open- 
loop" results. 


The results of our trajectory analysis with POST are listed below as 
a means of overview. Each specific point will be discussed in the fol- 
lowing sections. It is important to remember that while the following 
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results are generally applicable to the plane change problem, they must 
be considered with the qualifiers and limitations discussed above. That 
is, while the listed results are not uniquely attributable to POST, they 
are restricted to the plane change problem as previously defined. 

1. "Open-Loop" Results : The single-engine failure case caused the 
greatest error in end-conditions, with roughly a 25% error in final 
inclination, while the ±15% aerodynamic uncertainty produced a range of 
5-20% error in final inclination. The ±30% density perturbations, sur- 
prisingly, caused almost no end-condition errors. 

2. Heat-Rate Control : Proved to be very important for minimizing fuel 
usage and density compensation. In addition, the trajectory's sensitiv- 
ity to perturbations was a function of its maximum heat-rate. 

3. An gle-of-Attack : Remained almost nominal except during the cruise 
phase where it is adjusted to equalize the thrust-drag balance. 

4. Cruise Acceleration : Does not significantly change under any tested 
conditions and is practically zero. Fixing it at zero (i.e., constant 
velocity) has no observable effect. 

5. Early Corrections : From a fuel optimization standpoint, corrections 
should be made as early as possible. However, "fixing" the four entry 
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parameters only resulted in a maximum loss of 100 pounds of fuel effi- 


ciency. 


6. Time-to-Ex i t : Significant variations in 
to compensate for the various perturbations, 
fixed, the overall trend was to adjust the 
using nominal exit parameters. 


exit inclination were made 
With the entry parameters 
exit point and then exit 


7. Bank Angle Versus Angle-of-Attack : Typically bank was a more effi- 
cient (in terms of fuel usage) means of correcting for perturbations 
than angle-of -attack. 

8. Fuel Usage : Several items significantly effect the fuel usage. 
Aerodynamic characteristics, not surprisingly, had the largest effect 
followed by the time-to-exi t, maximum heat-rate, and the transition from 
entry to the heat-rate boundary during cruise. 

9. Exit Parameters : Tended to remain at near nominal values. In fact, 
angle-of-attack and bank angle after engine cutoff never significantly 
changed from the nominal values. 

4.4 DETAILED RESULTS 
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In the following subsections, the "open- loop" and 
optimization/compensation results for each perturbation type and pertur- 
bation combinations are discussed in detail. The "open-loops" results 
consist of the "open-loop" errors and the final vehicle weight for the 
specified conditions. The "open- loop" errors indicate how much devi- 
ation (from the desired value) in final inclination, apogee altitude or 
maximum throttle setting for a perturbed trajectory using nominal guid- 
ance parameter values (see Section 4.2). As a reminder, the desired 
end-conditions values are: 

1. Final inclination change - 20 degrees 

2. Final Apogee Altitude - 160 nm 

3. Maximum throttle setting - 100% 

The final vehicle weight provides an idea of how the perturbations 
affected fuel usage. However, significant "open-loop" errors will 
affect the usefulness of this result. As reference weights, the final 
vehicle weight (in a 160 nm circular orbit) for the baseline 80 and 125 
trajectories is 5262 lbs and 5873 lbs respectively. 

The opt i mi zation/compensation results consist of the optimized guid- 
ance parameter values and the best final vehicle weight for the speci- 
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fied conditions. Recall that these optimized guidance parameter values 
are the best independent variable values that POST can generate thru its 
optimization process for the specified conditons. Using these optimized 
independent variables as POST guidance parameters or commands will pro- 
duce a trajectory to achieve to desired end-conditions with the best 
final vehicle weight. By comparing these off-nominal optimized guidance 
parameters with the nominal guidance-parameters (also optimized for nom- 
inal conditions), insight can be gained on how to compensate for the 
dispersions in an optimum manner. Comparing the final vehicle weights 
(nominal verses off-nominal conditions), the net effect of the perturba- 
tion in terms of fuel usage can be obtained. 


4.4.1 



Trajectory analysis results for density perturbations indicated that 
heat-rate control offered "automatic" density compensation. This result 
was the most surprising result of our analysis. From an initial study 
it was determined that density dispersions would cause significant 
inclination and apogee altitude errors (19). In fact, this earlier 
study suggested that the extreme cases of density perturbations (±30%) 
could result i n 4 to 6 degrees (20-30%) of inclination error. However, 
as the results in the next subsection show, the ±30% density perturba- 
tions produced also no POST "open-loop" errors. This implies that 
either the trajectory is unaffected by density perturbations and our 
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initial study was incorrect or there is actually some automatic density 
compensator in the "open-loop" trajectory. 



To obtain the scope of the problem for density perturbations, POST 
runs were made using a 30% high density bias, a 30% low density bias and 
a "worst case" density situation where the density was 30% high until 
the start of the exit phase and then it switched to 30% low. The three 
cases were run for plane change missions using a maximum heat-rate of 80 
and 125. Overall, as the results in Table 4-1 demonstrate, density var- 
iations had almost no effect on the end-condition results for the "open- 
loop" trajectory. In fact, for a heating rate of 80 BTU/ft 2 -sec with a 
30% low density atmosphere, the resulting trajectory meets the desired 
end-conditions more accurately and with less fuel than under nominal 
(standard atmosphere) conditions. Even for the "worst case", the end- 
conditions were just outside of NASA's accuracy criteria. 


One final point that warrants mention concerning the "open-loop" 
results shown Table 4-1. For two of the cases, the maximum throttle 
setting exceeded 100%. To understand the significance of this point, 
further discussion on how the throttle setting is determined is in 
order. Except during the cruise phase, the throttle is fixed at either 
100% or zero. During the cruise phase, to maintain the cruise accerler- 
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Table 4-1. Density bias "open-loop" results with post. 


Condition 

Incl i nation 
Error 
(deg) 

Apogee 

Altitude 

Error 

(nm) 

Maximum 

Throttle 

Error 

(%) 

Fi nal 
Vehicle 
Wei ght 
(lbs) 

QDOT = 80 
Nomi nal 

0.004 

-0.27 

0 

5262 

QDOT = 80 
+30% 

0.004 

-0.28 

0 

5253 

QDOT = 80 
-30% 

-0.001 

0.07 

0 

5269 

QDOT = 125 
Nominal 

0.001 

0.20 

0 

5873 

QDOT = 125 
+30% 

0.129 

0.39 

0.80 

5847 

QDOT = 125 
-30% 

-0.214 

0.09 

0 

5904 

QDOT = 125 
+30% Until 
Exit Then 
-30% 

-0.670 

-5.38 

0.83 

59 7 8 


NOTE: A negative sign indicates the error was below the desired value. 
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ation, POST has the freedom to adjust the throttle. Initially, POST has 
the freedom to exceed 100% but during the optimization process it must 
eliminate the excess (real-world necessity) as one of its end-condition 
constraints. Therefore, the maximum throttle error is similar to incli- 
nation and apogee errors — some sort of adjustment will be needed to 
compensate for the error. While in both cases shown in Table 4-1, the 
throttle error was less than 1%, it is the first hint that during 
cruise, nominal ang 1 e-of-attack and 100% thrust will not maintain a con- 
stant velocity under some conditions. 

4. 4. 1.2 Op t imizat io n/C o mpensat ion. Results 

As might be expected the optimized trajectory results for density 
perturbations do not provide much insight since the "open-loop" errors 
were so small. As the summary of the results shown in Tables 4-2 thru 4 
indicate, many of the parameters do not change at all. The few changes 
that are made do not really stand out until analyzing the data for a 125 
heat-rate with the four entry parameters frozen (Table 4-4). In these 
cases, the data suggests that the best way to correct is with a vari- 
ation in exit inclination or a "time-to-exi t" routine. However, before 
dealing with smaller details of the results, the basic question of why 
such large density biases result in such small end-condition errors 
should be covered. 
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Table 4-2. Density bias optimization results for 
QDOT = 80 BTU/ft 2 -sec. 


Guidance Parameter 
Name (Units) 

Condition 

Nominal 

+30% 

-30% 

Geocentric Latitude (deg) 

0.957 

Same 

Same 

Perigee Altitude (nm) 

30.0 

Same 

30.1 

Entry Angle-of-Attack (deg) 

32.2 

Same 

33.5 

Entry Bank Angle (deg) 

-67.8 

Same 

-63.5 

Cruise Angle-of-Attack (deg) 

30.8 

Same 

31.3 

p 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

47.4 

Same 

47.2 

Exit Bank angle (deg) 

-59.1 

Same 

-60.8 

Exit Angle-of-Attack (deg) 

21.8 

Same 

Same 

Apogee Altitude (nm) 

19 3 

Same 

192 

Unpowered Exit 
Angle-of-Attack (deg) 

8.03 

Same 

8.05 

Unpowered Exit 
Bank Angle (deg) 

-68.5 

Same 

-68.7 

Final Vehicle Weight (lbs) 

5262 

5253 

529 3 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Table 4-3. Density bias optimization results for QDOT = 125 BTU/ft 2 -sec. 


Guidance Parameter 
Name (Units) 

Condition 

Nomi nal 


-30% 

+30% Until 
Exit Then 
-30% 

Geocentric Latitude (deg) 

-0.9 58 

Same 

Same 

Same 

Perigee Altitude (nm) 

27.6 

27.4 

27.5 

27.5 

Entry Angle-of -Attack 
(deg) 

27.2 

26.8 

26.8 

27.4 

Entry Bank Angle (deg) 

-69.2 

-71.1 

-71.4 

-69.4 

Cruise Angle-of -Attack 
(deg) 

25.4 

25.3 

25.3 

25.1 

Cruise Acceleration 
(ft/s 2 ) 

0.0036 

Same 

Same 

Same 

Exit Inclination (deg) 

40.6 

40.1 

40; 5 


Exit Bank angle (deg) 

-65.7 

-66.1 

-66.4 

-65.8 

Exit Angle-of-Attack 
(deg) 

20.4 

20.3 

Same 

20.6 

Apogee Altitude (nm) 

165 

Same 

Same 

169 

Unpowered Exit 
Angle-of-Attack (deg) 

8.13 

Same 

Same 

8.12 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

Same 

-69.7 

Final Vehicle Weight (lbs) 

5873 

59 41 

5943 

5842 


NOTE: (1) "Same" indicates the parameter value did not change from the 

nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Table 4-4. Density bias optimization results for QDOT = 125 BTU/ft^-sec. 
with entry parameters fixed. 


Guidance Parameter 
Name (Units) 

Condition 

Nominal 

+30% 

-30% 

Geocentric Latitude (deg) 

-0.958 

Fixed 

Fixed 

Perigee Altitude (nm) 

27.6 

Fixed 

Fixed 

Entry Angle-of -Attack (deg) 

27.2 

Fixed 

Fixed 

Entry Bank Angle (deg) 

-69.2 

Fixed 

Fixed 

Cruise Angl e-of-Attack (deg) 

25.4 

25.2 

25.2 

2 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

40.6 

39.2 

40.2 

Exit Bank angle (deg) 

-65.7 

-67.1 

-67.3 

Exit Angl e-of-Attack (deg) 

20.4 

20.3 

20.3 

Apogee Altitude (nm) 

165 

166 

166 

Unpowered Exit 
Angl e-of-Attack (deg) 

8.13 

Same 

Same 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same- 

Same 

Final Vehicle Weight (lbs) 

5873 

589 3 

5884 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 















































The most reasonable explanation is that POST is doing something 
automatically to control one thing and ends up also compensating auto- 
matically for density variations. Closer examination and tests clearly 
pointed out that the POST "heat-rate" controller (necessary to maintain 
the same plane change strategy) was also providing a built-in density 
compensator. Examination of Chapman's heat-rate (QDOT) equation 

QDOT = 17600/Rn (p/p SL ) °- 5 (V R /V C ) 3 15 

where Rn is the vehicle's nose radius, p and p SL are the current and sea 
level atmospheric density, V c is the reference circular orbital veloci- 
ty, and V R is the relative vehicle velocity, offers an explanation. By 
controlling to a given heat-rate you are actually controlling to a den- 
sity-velocity combination. There are infinite combinations of velocity 
and density to give the desired heat-rate, but with density providing 
the predominate deacceleration force during reentry, a straight change 
in density will only change the geometric altitude where the deacceler- 
ation profile occurs and eventually you will hit the heat-rate at rough- 
ly the same velocity-density combination. In other words, the vehicle 
will basically fly through the same density profile, but at a different 
geometric altitude. Thus, the heat-rate controller insures the same 
density profile which insures the same plane change results. So, while 
providing for safety of the vehicle, you are also getting compensation 
for density perturbations. 
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4.4.2 AERODYNAMIC UNCERTAINTIES 


As discussed in Chapter 2, the ERV is expected to perform the plane 
change maneuver in the region where previous aerodynamic predictions 
have been significantly in error. Therefore, the capability to handle 
uncertainties in the predicted aerodynamics is important for an accurate 
guidance system. To obtain the "open-loop" and 

optimization/compensation results, the following four variations of 
aerodynamic characteristics were considered. 

1. +15% CL and CD (Both High) 

2. -15% CL and CD (Both Low) 

3. +15% CL and -15% CD (High Aero) 

4. -15% CL and +15% CD (Low Aero) 

where CL is the coefficient of lift and CD is the coefficient of drag. 
These four cases represent the four possible extremes for a maximum 
aerodynamic uncertainty of ±15%. 
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4. 4. 2.1 "OPEN-LOOP 11 RESULTS 


The effects of the ERV's ±15% aerodynamic uncertainties produced a 
range of inclination errors of 5-20% depending of the combination of CL 
and CD as shown in Table 4-5. Overall, two important conclusions can be 
made. First, aerodynamic perturbations, unlike density perturbations, 
produce significant open-loop errors (e.g., fuel and final inclination) 
and will require some compensation to minimize their effects. Second, 
the problem of having sufficient power during level off was again high- 
lighted as potentially serious because the +15% CL and CD case (QDOT = 
125) required nearly 120% of maximum throttle to level -off using nominal 
guidance parameters. Some adjustment in a vehicle control parameter 
will be needed to correct this problem. 

4. 4. 2. 2 QPTIMIZAIQN/CQMPENSATION RESULTS 

Analyzing the results shown in Table 4-6 thru 8, four specific 
observations or insights are apparent concerning (1) fuel usage, (2) 
angl e-of-attack (3) optimum location for corrections, and (4) sensitiv- 
ity of the 80 and 125 heat-rate cases to the aerodynamic uncertainties. 

On fuel usage, the results of Tables 4-6 thru 8 fall into two cate- 
gories. For the case where CL and CD are either high or low (no net 
change of L/D) , it is possible to compensate for the perturbations with- 
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Table 4-5. Aerodynamic uncertainty "open-loop" results with post 


Condition 

Inclination 

Error 

(deg) 

Apogee 

Altitude 

Error 

(nm) 

Maximum 

Throttle 

Error 

(%) 

Final 
Vehicle 
Wei ght 
(lbs) 

QDOT = 80 
Nomi nal 

0.004 

-0.27 

0 

5262 

QDOT = 80 
+15% CL & CD 

0.322 

7.52 

0 

5306 

QDOT = 80 
-15% CL & CD 

-0.397 

-3.33 

0 

5073 

QDOT = 80 
+15% CL 
-15% CD 

-0.043 

25.5 

0 

5996 

QDOT = 80 
-15% CL 
+15% CD 

0.0387 

-12.1 

0 

4543 

QDOT = 125 
Nominal 

0.001 

0.20 

0 

5873 

QDOT = 125 
+15% CL & CD 

1.070 

-0.50 

19.7 

5721 

QDOT = 125 
-15% CL & CD 

-1.232 

-0.86 

0 

6041 

QDOT = 125 
+15% CL 
-15% CD 

-2.411 

4.96 

0 

69 2 6 

QDOT = 125 
-15% CL and 
+15% CD 

-3.681 

-7.00 

0 

569 1 


NOTE: A negative sign indicates the error was below the desired value. 
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Table 4-6. Aerodynamic uncertainty optimization results for 
QDOT = 80 BTU/ft^-sec. 



Condition 

Guidance Parameter 
Name (Units) 

Nomi nal 

+ 15 % CL 
-15% CD 

-15% CL 
+15% CD 

Geocentric Latitude (deg) 

0.9 57 

0.956 

0.9 59 

Perigee Altitude (nm) 

30.0 

29.2 

30.5 

Entry Angle-of -Attack (deg) 

32.2 

32.1 

32.5 

Entry Bank Angle (deg) 

-67.8 

Same 

-67.4 

Cruise Angle-of -Attack (deg) 

30.8 

30.2 

30.9 

9 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

47.4 

47.5 

47.3 

Exit Bank angle (deg) 

-59.1 

-59.7 

-58.3 

Exit Angle-of -Attack (deg) 

21.8 

21.6 

22.1 

Apogee Altitude (nm) 

193 

186 

201 

Unpowered Exit 
Angle-of-Attack (deg) 

8.03 

8.04 

8.00 

Unpowered Exit 
Bank Angle (deg) 

-68.5 

-68.7 

-68.4 

Final Vehicle Weight (lbs) 

5262 

6029 

4521 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Table 4-7. Aerodynamic uncertainty optimization results for 
QDOT = 125 BTU/ft^-sec. 


Guidance Parameter 
Name (Units) 

Condition 

Nominal 

+15% 

CL & CD 

-15% 

CL & CD 

+15% CL 
-15% CD 

-15% CL 
+15% CD 

Geocentric 
Latitude (deg) 

-0.9 58 

Same 

Same 

-0.948 

Same 

Perigee 
Altitude (nm) 

27.6 

27.8 

27.5 

29.4 

28.3 

Entry Angle-of- 
Attack (deg) 

27.2 

26.2 

27.4 

26.7 

27.7 

Entry Bank Angle 
(deg) 

-69.2 

-68.5 

-69.6 

-73.3 

-70.4 

Cruise Angle-of- 
Attack (deg) 

25.4 

23.3 

Same 

25.3 

Same 

Cruise Accelera- 
tion (ft/s ) 

0.0036 

Same 

Same' 

Same 

Same 

Exit Inclination 
(deg) 

40.6 

Same 

41.5 

41.4 

44.2 

Exit Bank angle 
(deg) 

-65.7 

-65.3 

-66.3 

-69.8 

-65.8 

Exit Angle-of- 
Attack (deg) 

20.4 

Same 

20.5 

20.3 

20.9 

Apogee Altitude 
(nm) 

165 

164 

166 

166 

172 

Unpowered Exit 

Angle-of-Attack 

(deg) 

8.13 

8.12 

Same 

8.12 

8.12 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

Same 

Same 

-69.7 

Final Vehicle 
Weight (lbs) 

5873 

589 2 

5819 

6682 

4886 


NOTE: (1) "Same" indicates the parameter value did not change from the 

nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Table 4-8. Aerodynamic uncertainty optimization results for 

QDOT = 125 BTU/ft^-sec. with entry parameters fixed. 


Guidance Parameter 
Name (Units) 

Condition 

Nomi nal 

+15% CL 
-15% CD 

-15% CL 
+15% CD 

Geocentric Latitude (deg) 

-0.958 

Fixed 

Fixed 

Perigee Altitude (nm) 

27.6 

Fixed 

Fixed 

Entry Angle-of -Attack (deg) 

27.2 

Fixed 

Fixed 

Entry Bank Angle (deg) 

-69.2 

F i xed 

Fixed 

Cruise Angle-of-Attack (deg) 

25.4 

24.8 

24.6 

p 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

40.6 

41.4 

42.5 

Exit Bank, angle (deg) 

-65.7 

-72.1 

-74.3 

Exit Angle-of-Attack (deg) 

20.4 

20.1 

21.3 

Apogee Altitude (nm) 

165 

168 

185 

Unpowered Exit 
Angle-of-Attack (deg) 

8.13 

Same 

8.11 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

-69.7 

Final Vehicle Weight (lbs) 

5873 

6584 

4981 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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out a net effect on fuel usage. For the cases where you have either 
high aero or low aero, the optimized trajectory had a roughly 10-15% net 
increase or decrease in fuel usage or about 700 pounds of fuel. Neither 
of these observations are unexpected but should be kept in mind when 
evaluating the final guidance algorithm. 

Concerning ang 1 e-of -attack (alpha) control, the results of Tables 
4-6 thru 8 fall into two categories. First, for the case of both high 
CL and CO in Table 4-7, a change in alpha was able to eliminate the max- 
imum throttle error (not shown in table). Thus, angl e-of-attack is a 
good candidate for handling thrust deficiencies. Second, throughout the 
optimization process for aerodynamic uncertainties, bank angle was 
adjusted significantly more than alpha to eliminate inclination errors. 
This is especially apparent in Table 4-8 where the entry parameters are 
not allowed to change. Examination of the partial derivatives of incli- 
nation with respect to alpha and bank (not shown in tables), indicates 
that the effectiveness of alpha verses bank angle is about the same. 
However, the fuel efficiency is apparently not the same because the 
optimization process adjusts almost strictly with bank angle during 
entry and exit. For example, Table 4-8 shows that during the powered 
exit phase, bank and alpha for nominal conditions are respectively -65.7 
and 20.4, while for the high aero case, bank and alpha are -72.1 and 
20.1. In other words, a change of 6.4 degrees in bank but only 0.3 
degrees in alpha is used to optimize for the off-nominal condition. 
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This implies that bank control is more fuel efficient for aerodynamic 
perturbation corrections than alpha. 

Concerning the optimum location for corrections, one general trend 
runs through all the results thus far (density and aerodynamic). Look- 
ing at where POST modified its guidance parameters to correct for errors 
caused by perturbations, it is apparent that an "ear 1 i er- i s-better" 
logic prevails. That is, most of the changes in the independent vari- 
ables are made at the front end of the trajectory. In addition, the 
si 2 e of these changes tend to hide their true significance since the 
effectiveness of corrections (sensitivity) is much higher on the front 
end of the trajectory. Particularly surprising was how little the exit 
parameters were changed by the optimization process. This observation 
suggest that corrections for perturbations errors should be made early 
and then exit nominally. Bear in mind though that POST "knows" about 
the perturbations before hand and can correct early. This apriori know- 
ledge is not normally available "real-world", thus the results are not 
completely applicable to a real situation. Table 4-8 shows the opti- 
mized parameters after "fixing" the first four. This requires POST to 
optimize for perturbations without the ability to adjust the entry 
parameters. Just as before though, the results show most of the adjust- 
ments were made as soon as possible (i.e., during the cruise phase), 
and then the exit is made with a near nominal command profile. Thus, 
for fuel conservation, corrections should be made as early as possible. 
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As a final observation, the 125 BTL)/ft 2 -sec trajectories are much 
more sensitive to perturbations than the 80 BTU/ft 2 -sec trajectories. 
This trend is apparent in the density perturbation results, but is espe- 
cially highlighted in the aerodynamic uncertainty results. This insight 
will be further expanded in the next subsection. 

4.4.3 HEAT-RATE VARIATION 

As discussed in Chapter 2, the guidance algorithm should be able to 
handle plane change maneuvers at a heating-rate of both 80 and 125. 
During the early trajectory analysis work, both 80 and 125 BTU/ft 2 -sec 
heat-rates were examined. However, a review of Table 4-1 and 4-5 indi- 
cates the relative insensitivity of a 20 degree plane change done at a 
heating rate of 80 versus 125 to both density and aerodynamic perturba- 
tions. Specifically, for the 80 case, the perturbations have almost no 
effect (i.e., less than 2 %) on final inclination, little effect on fuel 
usage (10% in the most extreme cases) and only a moderate effect on the 
final apogee altitude. In contrast, the "open-loop" errors for the 125 
case are significantly larger. 

The only exception to this trend is with apogee altitude. However, 
the apparent exception to the sensitivity trend is not really as it 
appears. It can be explained by the large difference in engine cutoff 
altitudes for the two trajectories. Examination of Figures A-8 and A-1 
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shows that the engines were cutoff at 255,000 feet for the 80 
BTll/ft 2 -sec trajectory but at 281,000 feet for the 125 trajectory. This 
difference in cutoff altitude magnifies the perturbation effects result- 
ing in a much larger apogee error for the 80 trajectory. In other 
words, the perturbation effects are not really larger for the 80 case, 
they just have longer to propagate. While this particular difference in 
the trajectories must be considered in developing a guidance algorithm, 
the 125 heat-rate trajectory will drive most of the design requirements. 

Actually, the sensitivity trend is not surprising in light of the 
work done earlier at NASA LaRC which demonstrated that each plane change 
had an optimum heat-rate (using fuel usage as an efficiency criterion) 
(6). From that study, it was shown that for a 20 degree plane change, a 
heating rate of 125 was near-optimum for the ERV, while a heating rate 
of 80 was significantly less efficient. Thus it is reasonable that den- 
sity and aerodynamic perturbations would have a larger effect on the 
optimum heat-rate trajectory (125) than on other trajectories.. 

A more quantitative explanation can be made by comparing the dynamic 
pressure profiles shown in Figure A-11 for the 80 and 125 trajectories 
(nominal conditions). Not only does the 125 trajectory reach a much 
higher maximum value but the time- integrated total dynamic pressure is 
significantly higher (area under the curve). Considering the direct 
relationship between aerodynamic forces and dynamic pressure. 
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aerodynamic = dynamic x aerodynamic x surface 
forces pressure coefficients area 

the sensitivity differences between the 80 and 125 trajectories is 
straight foward. The 125 plane change uses less fuel because it uses 
more aerodynamic forces to assist in the plane change. Therefore, per- 
turbations that effect the aerodynamic forces (e.g., density, CL or CD) 
will cause a more significant impact on the 125 trajectory. 

The final guidance scheme must demonstrate satisfactory performance 
for both heat-rates, and Chapter 6 provides that demonstration. Howev- 
er, in the interest of clarity of discussion and because the 125 heat- 
rate is generally the more difficult case, the remaining discussion of 
this chapter will deal only with the 125 case. 

4.4.4 SINGLE-ENGINE FAILURE 

The purpose for considering engine failure is more a question of 
safety than anything else. That is, we do not expect to experience an 
engine failure as we do expect density and aerodynamic dispersions. So 
the first question is can the ERV be safely returned to orbit even if 
the desired plane change is not completed. To obtain that answer, two 
cases were considered that were thought to be the two worst possibil- 
ities for single-engine failure. The first case was where failure 
occurs at the start of the level off phase and the second case at the 
start of the exit phase. In both cases, the thrust from one of the 
ERV's three engines was eliminated for the duration of the mission. 
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4. 4. 4.1 "OPEN-LOOP" RESULTS 


As might be expected, the loss of an engine during the plane change 
has by far the largest effect on end-condition errors, but the maneuver 
can still be accomplished. Table 4-9 shows the "open-loop" errors for 
the two, single-mngine failure cases. The inclination and apogee errors 
are initially the same for both situations because the exit command is 
based on an inclination angle. However, the loss of an engine at the 
start of the level-off phase is the more difficult case because of the 
power deficiency as highlighted by the maximum throttle error of 50%. 
The other significant error is the roughly 25% error in final inclina- 
tion. The fact that too much plane change was obtained might seem sur- 
prising unless you realized that, "open-loop", the exit will be at the 
same nominal inclination. Therefore, with a slower climb back to orbit 
(less power), the extra time in the atmosphere will produce a greater 
turn and hence a greater plane change. 

4. 4. 4. 2 OPTIMIZATION/ COMPENSATION RESULTS 

Table 4-10 presents the results for the optimized trajectory with a 
single-engine failure. Analysis of these results brings out two recur- 
ring observations and one new piece of insight. 
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Table 4-9. Single engine failure "open-loop" results with post for 
QDOT = 125 BTU/ft^-sec. 


Condition 

Inclination 

Error 

(deg) 

Apogee 

Altitude 

Error 

(nm) 

Maximum 

Throttle 

Error 

(%) 

Final 
Vehicle 
Wei ght 
(lbs) 

Nomi nal 

0.001 

0.20 

0 

5873 

Single Engine 
Failure at 
Level Off 

4.651 

-2.82 

49.9 

509 2 

Single Engine 
Failure at 
Exit 

4.651 

-2.82 

0 

509 2 


NOTE: A negative sign indicates the error was below the desired value. 
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First of all, concerning angle-of-attack control, the results demon- 
strate that a thrust deficiency can be corrected with alpha. As stated 
before, some sort of scheme will be needed to correct a power deficien- 
cy, and angle-of-attack is capable of handling even this worst case of 
an engine out. Furthermore, this means of compensation for a loss of 
power is evidently fuel efficient since the optimized results for an 
engine failure (Table 4-10) show a final vehicle weight that is almost 
the same as for nominal conditions. 

The second recurring trend is the preference for bank instead of 
alpha for inclination corrections. This trend is especially evident in 
Table 4-10 when comparing bank and alpha during the exit phase for an 
engine failure at level-off. When the entry parameters are fixed, the 
optimum means of compensation is to change the bank by 3.5 degrees but 
alpha by only 0.1 degrees. Again, the sensitivity of bank and alpha are 
almost the same (the partial derivatives are not shown), but evidently, 
the fuel efficiency is not, because bank is strongly favored. 

Finally, the new bit of insight deals with determination of the time 
or point to exit. The results of Table 4-10 indicate that exit inclina- 
tion was adjusted during the optimization process as much as 2.3 
degrees. This represents a significant modification to the trajectory 
and can be explained by the basic plane change strategy. Recall from 
Chapter 2 that as much of the plane change should be done during entry 
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Table 4-10. Single engine failure optimization results for 
QDOT = 125 BTU/ft2-sec. 



Condition 

finirlanrp Parampfpr 


Single Engine Failure 

Name (Units) 

Nominal 

At Level - 
Off 

At Exit 

At Level - 
Off 

Geocentric Latitude (deg) 

-0.9 58 

-0.961 

Same 

Fixed 

Perigee Altitude (nm) 

27.6 

28.7 

27.7 

Fixed 

Entry Angle-of -Attack 
(deg) 

27.2 

25.6 

26.6 

Fixed 

Entry Bank Angle (deg) 

-69.2 

-69.6 

-68.0 

Fixed 

Cruise Angle-of -Attack 
(deg) 

25.4 

20.4 

Same 


Cruise Acceleration 
(ft/s 2 ) 

0.0036 

Same 

Same 

Same 

Exit Inclination (deg) 

40.6 

38.3 

39.7 

38.8 

Exit Bank angle (deg) 

-65.7 

-65.9 

-62.8 

-62.2 

Exit Angle-of -Attack 
(deg) 

20.4 

Same 

20.3 

20.3 

Apogee Altitude (nm) 

165 

166 

167 

Same 

Unpowered Exit 
Angle-of -Attack (deg) 

8.13 

Same 

Same 

Same 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

Same 

Same 

Final Vehicle Weight (lb) 

5873 

5827 

5827 

5732 


NOTE: (1) "Same" indicates the parameter value did not change from the 

nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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and exit to minimize fuel usage. Thus, this shift in exit just repres- 
ents the optimization process adjusting to achieve this objective. A 
similar capability will probably be needed for an efficient guidance 
scheme. 

4.4.5 MULTIPLE OFF-NOMINAL CONDITIONS 

Having discussed analysis results on each individual type of pertur- 
bation, it is also worthwhile to discuss various combinations of pertur- 
bations. Although, little additional insight was gained by combining 
various perturbations, the fact that these combinations did not cause 
any new problems is in itself, a very useful insight. 

4. 4. 5.1 ■■QPEN-LQQP 1 ' RESULTS 

Although a large number of perturbation combinations are possible, 
the following five combinations were chosen to scope the problem. 

1. 30% high density and high aero 

2. 30% high density and low aero 

3. 30% low density and high aero 
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4. 30% low density and low aero 


5. 30% low density, low aero, and single-engine failure at 

level -off 

The "open-loop" errors for these runs are summarized in Table 4-11. For 
these cases, the trajectory problems caused by the combination of per- 
turbations did not produce any errors larger than had been seen previ- 
ously and in some cases the perturbations seem to neutralize each other. 
Overall, inclination errors were within 25% of nominal, apogee errors 
within 10%, and maximum throttle errors within 23%. 

4. 4. 5. 2 OPTIMIZATION/ COMPENSATION RESULTS 

With but one exception, the insight gained from the results of opti- 
mizing the five cases are mainly a repeat of what was discussed in the 
previous sections. The new point will be discussed in detail while the 
old points will only be briefly covered. 

4. 4. 5. 2.1 HEAT-RATE CONTROL 

4. 4. 5. 2. 1.1 TRANSITION TO CRUISE 
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Table 4-11. Multiple off-nominal "open-loop" results with post for 
QDOT = 125 BTU/ft2-sec. 


Condition 

Inclination 

Error 

(deg) 

Apogee 

Altitude 

Error 

(nm) 

Maximum 

Throttle 

Error 

(%) 

Final 

Vehicle 

Weight 

(lbs) 

Nomi nal 

0.001 

0.20 

0 

5873 

+30% Density 
+15% CL 
-15% CD 

-2.362 

5.29 

0 

69 22 

+30% Density 
-15% CL 
+15% CD 

-4.471 

-13.1 

0 

5819 

-30% Density 
+15% CL 
-15% CD 

-2.523 

4.34 

0 

69 38 

-30% Density 
-15% CL 
+15% CD 

-0.166 

-2.58 

11.9 

5158 

-30% Density 
-15% CL 
+15% CD 
Single Engine 
Failure 

-4.711 

8.04 

22.8 

5915 


NOTE: A negative sign indicates the error was below the desired value. 
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Our analysis showed that the transition from entry to the cruise 
phase could effect the plane change fuel efficiency by as much as seven 
or eight percent of the total fuel required for a plane change. 

4.4.5.2.1 .1 .1 THE PROBLEM 

While analyzing the results, for multiple off-nominal conditions, we 
realized that a couple of the trajectories had lost the characteristic 
(for the 125 case) "double-peak" in the heat-rate verse time curve. 
Further analysis indicated that if the "double-peak" (see Figure 2-13 or 
A— 9) is eliminated or altered, extra fuel will be required to complete 
the mission. For example, the results in Table 4-12 show that the opti- 
mized trajectory for 30% high density and 15% low aero has a final vehi- 
cle weight at 4839 pounds. However, by slightly modifying the 
trajectory (the modification will be covered in the next subsection), 
the heat-rate "double-peak" is restored, and the final vehicle weight is 
now 5281 pounds for a net fuel savings of roughly 450 pounds. The heat- 
rate verses time profiles for the two cases are shown in Figure 4-1. 
Further study suggests two reasons for this fuel savings, both a benefit 
of the relatively high dynamic pressure during this phase of flight. 

The first reason can be seen by comparing the plots in Figure 4-2 
which show the dynamic pressure profiles for the unmodified (original) 
trajectory and the modified trajectory (both trajectories are optimized 
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Table 4-12. Multiple off-nominal optimization results for 

QDOT - 125 BTU/ft^-sec. with entry parameters fixed. 


Condition 


Guidance Parameter 
Name Hlnitcl 

Nnmi na 1 

+30% Density 
-15% CL and +15% CD 



Unmodi fied 

Modified 

Geocentric Latitude (deg) 

-0.958 

Fixed 

Fi xed 

Perigee Altitude (nm) 

27.6 

Fi xed 

Fi xed 

Entry Angle-of -Attack (deg) 

27.2 

Fixed 

Fixed 

Entry Bank Angle (deg) 

-69.2 

Fi xed 

Fi xed 

Cruise Angle-of -Attack (deg) 

25.4 

25.3 

23.1 

p 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

40.6 

44.2 

35.8 

Exit Bank angle (deg) 

-65.7 

-69.2 

-66.5 

Exit Angle-of-Attack (deg) 

20.4 

21.4 

Same 

Apogee Altitude (nm) 

165 

180 

168 

Unpowered Exit 
Angle-of-Attack (deg) 

8.13 

8.12 

Same 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

-69.7 

Same 

Final Vehicle Weight (lbs) 

5873 

4839 

5281 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Figure 4 


ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 



ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 



-1. Heat-rate profiles for +30% density, -15% CL and +15% CD. 
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as previously explained for 30% high density and 15% low aero). The 
curves are similar except the first curve does not extend above 60 
lbs/ft 2 in dynamic pressure, while the second curve has a section that 
goes almost to 85 lbs/ft 2 . Overall, the maximum and the time-integrated 
total dynamic pressure is less for the unmodified trajectory. This sit- 
uation is similar to the comparison between the 80 and 125 BTU/ft 2 -sec 
trajectories that was covered in Section 4.6. The lower dynamic pres- 
sure results in lower aerodynamic forces, hence less turning capability, 
which results in more of the plane change being done with propulsive 
forces. Thus, a similar fuel difference is seen between the unmodified 
and modified trajectories as between the 80 and 125 trajectories. 

A second explanation for the fuel difference has to do with the 
plane change strategy itself. Recall from previous explanations of this 
strategy, that during entry as the vehicle approaches the maximum heat- 
rate, a level -off is started. For the unmodified trajectory, the lev- 
el -off was started as the heat-rate curve was approaching what should 
have been its first heat-rate peak. In fact, closer examination of Fig- 
ure 4-1 (unmodified plot) shows a slight change in the slope at 115 
BTU/ft 2 -sec. For the modified trajectory the level-off was not started 
until the heat-rate curve was approaching a second peak. 

By starting the level off at the first heat-rate peak, the vehicle 
reaches the maximum heat-rate boundary earlier than necessary, and the 
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ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 


- 'M‘ 



ERV CRUISING SYNERGISTIC PLANE CHANGE 0-D0T=125 



0 350 700 1050 1400 1750 2100 

Time, sec 


Figure 4-2. Dynamic pressure profiles for +30% density, -15% CL and 
+15% CD. 
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engines are used earlier than necessary. Recall that the engines are 
started as a part of the level -off. A comparison of the plots in Figure 
4-3 shows how for the unmodified trajectory, the ERV is initially com- 
manded to increase its bank angle (driven onto the heat boundary) before 
decreasing the bank angle to join the heat-rate boundary. A comparison 
of the plots in Figure 4-4 shows how the engines are started 210 seconds 
earlier for the unmodified trajectory but only cutoff 95 seconds earli- 
er. This results in a larger time-integrated total, which translates 
into a fuel difference between the two trajectories. In essence, you 
are throwing away this phase of free, relatively high dynamic pressure 
which could be used to give a high turning rate at no cost in propel- 
lant. In addition, the higher dynamic pressures that are obtained by 
delaying your level-off (modified trajectory) are never obtained. Thus, 
the problem that was discussed in the first explanation. 

It should be pointed out that this alteration of the heat-rate curve 
also occurred for the 15% low aero case. This helps explain the large 
decrease (from nominal) in final vehicle weight that is shown in Table 
4-7. In other words, this transition to cruise problem is not 
restricted only to multiple off-nominal conditions. 

4. 4.5. 2.1 .1 .2 THE SOLUTION 




97 


ERV CRUISING SYNERGISTIC PLANE CHANGE 0-D0T=125 
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Figure 4-3. Bank angle profiles for +30% density, -15% CL and +15% CD 


ETA ETA 


ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 


Unmodified Transition 
(Lead-point = 1 15 Btu/ft^-sec) 



Time, sec 


ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 


Modified Transition 
(Lead-point = 120 Btu/ft^-sec) 


Time, sec 

Figure 4-4. Engine throttling profiles for +30% density, -15% CL and 
+15% CD (0 - OFF, 1.0 - MAXIMUM). 




The difference between the modified and unmodified trajectories was 
achieved by changing the lead-point used by the heat-rate controller. 
For a dynamic environment, to avoid overshooting the desired heat-rate, 
a controller uses a value that is different than the desired heat-rate. 
This value can be called the lead-point value or simply lead-point. For 
all previous 125 trajectories, the lead-point was set at 115 
BTU/f t 2 -sec . The only change for the modified trajectory was to raise 
the lead-point to 120 BTU/f t 2 -sec. This allowed the restoration of the 
"double-peaked" heat-rate curve and the 450 pound fuel savings. 

This simple lead-point modification could not fully restore the 
first heat-rate peak for all perturbation combinations, but it did pro- 
vide a progressive improvement in fuel usage (i.e., the more restora- 
tion, the more fuel savings). Figure 4-5 provides a comparison of a 
modified and unmodified trajectory as an example of partial restoration. 
The scales for the plots are not the same so an exact comparison is dif- 
ficult. However, it is evident that some restoration is made. A com- 
parison of the final vehicle weights in Table 4-13 indicates that nearly 
200 pounds of fuel was recovered for this example. 

There is another possible modification that could be used since the 
lead-point modification does not always work. The loss of fuel effi- 
ciency is caused by either a lead-point that is too low or off-nominal 
conditions that raise the occurrence of the first heat peak up to the 
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ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 


Unmodified Transition 



ERV CRUISING SYNERGISTIC PLANE CHANGE Q-DOT=125 



Figure 4-5. Heat-rate profiles for -30% density, -15% CL and +15% CD. 
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Table 4-13. Multiple off-nominal optimization results for 

QOOT = 125 BTU/ft^-sec. with entry parameters fixed. 


Guidance Parameter 
Nflmp (Ilni + <;1 


Condition 


Nnmi nal 

-30% Density 
-15% CL and +15% CD 



Unmodified 

Modi fied 

Geocentric Latitude (deg) 

-0.9 58 

Fixed 

Fixed 

Perigee Altitude (nm) 

27.6 

Fixed 

Fixed 

Entry Angle-of-Attack (deg) 

27.2 

Fixed 

Fixed 

Entry Bank Angle (deg) 

-69.2 

Fixed 

Fixed 

Cruise Angle-of-Attack (deg) 

25.4 

24.1 

24.0 

p 

Cruise Acceleration ( ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

40.6 

41.2 

37.2 

Exit Bank angle (deg) 

-65.7 

-67.1 

-70.9 

Exit Angle-of-Attack (deg) 

20.4 

20.6 

19 .9 

Apogee Altitude (nm) 

165 

168 

173 

Unpowered Exit 

8.13 

Same 

8.12 

Angle-of-Attack (deg) 




Unpowered Exit 

-69.8 

Same 

Same 

Bank Angle (deg) 




Final Vehicle Weight (lbs) 

5873 

5024 

5222 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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nominal lead-point. Thus, the lead-point needs to be carefully chosen 
to avoid the maximum level of the first heat peak for off-nominal condi- 
tions. However, if this is not possible, active control to maintain the 
first heat peak below the lead-point might be possible. This option 
will be discussed in Chapter 7. 

4. 4. 5.2. 1.2 HEAT-RATE UNCERTAINTY 

It has been shown several times that the efficiency of the synerget- 
ic plane change is very much a function of heat-rate and how accurately 
it can be controlled. In fact, for the latest example, a modification 
of the lead-point by only 5 BTU/ft 2 -sec offered a significant fuel sav- 
ings. Chapman's equation, defines heat-rate as a function of velocity 
and density. Therefore, to achieve accurate heat-rate control, accurate 
velocity and density measurements are a necessity. Previous schemes 
have used other control boundaries such as constant drag because it can 
be measured quite accurately and avoids the uncertainty in density. 
However, flying a constant drag level does not insure that heat-rate is 
accurately controlled or constrained. To remain within a heat-rate con- 
straint, a constant drag level still must be far enough from the heat- 
rate limit to account for the uncertainties. Therefore, no advantage is 
gained by using a controller other than heat-rate. The bottom line is 
you still must deal with the heat-rate uncertainty if you wish to 
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achieve the performance gains offered by a maneuver like the synergetic 
plane change. 



Tables 4-12 thru 15 indicate, as have the previous optimized tables, 
that there is very little variation between the optimized ang 1 es-of -at- 
tack for nominal conditions and the perturbated conditions except during 
the cruise phase. However, during cruise, the ang 1 e-of -attack is sig- 
nificantly adjusted in the cases where there was an "open- loop" maximum 
throttle error. In all such cases, enough modification in the cruise 
alpha was made to eliminate the throttle error at no cost in terms of 
fuel usage. 



The preference for bank angle over angl e-at-attack for adjustments 
during entry or exit is again highlighted in Table 4-12 and 13. In all 
four off-nominal cases, the exit bank angle is increased anywhere from 
0.8 to 5.2 degrees while ang 1 e-of -attack at most (except for one unmodi- 
fied case) is decreased by 0.5 degrees. Again, the partial derivatives 
(not shown in the tables) indicate that the effectiveness of angle-of- 
attack and bank angle are about the same. Thus, the preference for bank 
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Table 4-14. Multiple off-nominal optimization results for 
QDOT = 125 BTU/ft2-sec. 


Condition 


finirianrp Parampt.pp 


+30% Density 

Name (Units) 

Nominal 

+15% CL 
-15% CD 

-15% CL 
+15% CD 

Geocentric Latitude (deg) 

-0.958 

Same 

Same 

Perigee Altitude (nm) 

27.6 

Same 

27.8 

Entry Angle-of -Attack (deg) 

27.2 

27.4 

28.0 

Entry Bank Angle (deg) 

-69.2 

-70.1 

-69.0 

Cruise Angle-of-Attack (deg) 

25.4 

Same 

24.6 

2 

Cruise Acceleration (ft/s ) 

0.0036 

Same 

Same 

Exit Inclination (deg) 

40.6 

42.4 

43.7 

Exit Bank angle (deg) 

-65.7 

-66.7 

-66.5 

Exit Angle-of-Attack (deg) 

20.4 

Same 

20.6 

Apogee Altitude (nm) 

165 

163 

171 

Unpowered Exit 
Angle-of-Attack (deg) 

8.13 

Same 

Same 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

Same 

Final Vehicle Weight (lbs) 

5873 

6559 

4941 


NOTE: (1) "Same" indicates the parameter value did not change from 

the nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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Table 4-15. Multiple off-nominal optimization results for 
QDOT = 125 BTU/ft 2 -sec. 


Condition 





-30% Density 

! 

ftinHanrp Paramptor 


+15% CL 

-15% CL and +15% CD 

Name (Units) 

Nomi nal 

-15% CD 


1 Engine 
Out 

Geocentric Latitude (deg) 

-0.9 58 

-0.960 

Same 

-0.9 59 

Perigee Altitude (nm) 

27.6 

28.1 

27.5 

28.0 

Entry Angle-of -Attack 
(deg) 

27.2 

27.2 

26.7 

21.8 

Entry Bank Angle (deg) 

-69.2 

-69.0 

Same 

-68.3 

Cruise Angle-of -Attack 
(deg) 

25.4 

25.5 

24.1 

22.8 

Cruise Acceleration 
(ft/s 2 ) 

0.0036 

Same 

Same 

Same 

Exit Inclination (deg) 

40.6 

42.9 

40.6 

44.9 

Exit Bank angle (deg) 

-65.7 

-67.8 

-67.1 

-66.1 

Exit Angle-of -Attack 
(deg) 

20.4 

20.5 

20.5 

20.8 

Apogee Altitude (nm) 

165 

164 

168 

169 

Unpowered Exit 
Angle-of -Attack (deg) 

8.13 

Same 

Same 

8.12 

Unpowered Exit 
Bank Angle (deg) 

-69.8 

Same 

Same 

-69.7 

Final Vehicle Weight (lb) 

5873 

6474 

5124 

4882 


NOTE: (1) "Same" indicates the parameter value did not change from the 

nominal value. 

(2) "Fixed" indicates the parameter was not allowed to change 
during the optimization process. 

(3) A negative sign indicates a left bank angle. 
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angle implies that bank angle is more fuel efficient in compensating for 
the perturbations in question. 

4. 4. 5. 2. 4 EARLY CORRECTIONS 

While the results of Table 4-12 thru 15 are not as clear as the pre- 
vious results, there is still a trend of making adjustments early and 
flying the remaining trajectory with near nominal guidance parameters. 
Again the effeciveness of alpha and bank changes is much higher during 
entry so the preference for early changes is somewhat hidden by the size 
of the changes. 

4. 4. 5. 2 . 5 TIHE-TQ-EXIT 

Tables 4-12 thru 15 show a significant variation in the exit incli- 
nation. Actually, this result is not especially surprising since the 
closer a plane change is done to the node and the quicker it is per- 
formed, the more efficient it tends to be. Therefore, the variation in 
exit inclination is just a result of the optimization program adjusting 
for changes in speed and location of the plane change for the perturbed 
conditions. Thus, a guidance scheme will certainly need to make similar 
types'of adjustments to obtain a reasonable degree of fuel efficiency. 
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CHAPTER 5 


GUIDANCE SCHEME 


5.1 INTRODUCTION 

Based on the guidance concepts discussed in Chapter 3 and the tra- 
jectory analysis covered in Chapter 4, a guidance scheme designed to 
guide the ERV through the atmospheric phase of a synergetic plane change 
maneuver is suggested. This guidance scheme by design produces a tra- 
jectory very similar to the plane change trajectory generated by POST. 
As shown in Chapter 4, this trajectory provides a very good baseline 
that requires relatively few modifications to maintain its optimality 
even in the presence of significant density and vehicular perturbations. 
Therefore, the suggested algorithm mirrors the strategy used in POST 
with a few simplifications (e.g., no optimization cycle) and several 
technique refinements. The resulting trajectory can be divided into 
three phases: (1) An entry phase where angl e-of-attack and bank are 
held constant until pullout, (2) A powered cruise phase where the maxi- 
mum heating rate is maintained and the time to start the exit phase is 
determined and (3) The exit phase where apogee and inclination angle 
are controlled with thrust and bank angle respectively. Figure 5-1 
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shows an altitude profile, for a target heat-rate of 125 BTU/ft 2 -sec and 
nominal conditions, with the three trajectory phases labeled. Within 
these three phases are three primary guidance controllers that form the 
hybrid guidance system discussed in Chapter 3. (1) A heat-rate control- 
ler using perturbation guidance techniques, (2) A time-to-exit control- 
ler using predictor logic, and (3) An end-condition controller using 
predictor/corrector techniques. Fuel is not explicitly optimized. The 
remaining two sections of this chapter will discuss in detail the guid- 
ance algorithm and the rationale for choosing the suggested scheme. 
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FRAME 003 



Figure 5-1. Guidance algorithm phases. 
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5.2 GUIDANCE ALGORITHM 


This section details the guidance algorithm by individually discuss- 
ing each of the three trajectory phase. Figure 5-2 provides a summary 
of the guidance system and should be helpful as the algorithm is cov- 
ered. All the symbols used in this chart are explained in the following 
subsect i ons . 

5.2.1 ENTRY PHASE 

The entry phase is a constant attitude gliding maneuver that starts 
at atmosphere interface (400,00 feet) and ends at the pullout point for 
the cruise phase. A deorbit burn is made to achieve the desired veloci- 
ty and position at 400,00 feet. The deorbit burn, perigee altitude, 
angle-of-attack, and bank angle are determined pre-mission to provide 
the most fuel efficient entry conditions for nominal conditions. The 
entry angle-of-attack and bank angle are set at 400,00 feet and remain 
constant throughout the entry phase. The pullout point is defined as 
the point where the vehicle reaches 90% of the desired heat-rate. Heat- 
rate (QDOT) is calculated on board using Chapman's equation. 


QD0T = 17600 (P ra /P SL ) 0 5 (V r /26000) 3 - 15 
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Figure 5-2. Synergetic plane change guidance system summary. 

















where p ra is the measured atmospheric density, p SL is the sea-ievel atmo- 
spheric density, and V R is the vehicle's velocity relative to the atmos- 
phere in feet/second 2 . The density measurement is determined by direct 
measurement from an air data system like the Shuttle Entry Air Data Sys- 
tem. Velocity is obtained from the inertial measurment unit (IMU) . 

5.2.2 CRUISE PHASE 

The cruise phase begins when the heat-rate reaches 90% of the 
desired heat-rate and ends when exit is commanded by the time-to-exit 
predictor. During this phase, the ERV transitions from a descending 
glide to a powered, 1 evel -f 1 i ght, constant-velocity cruise to maintain 
the desired heat-rate. The level-off is initiated by lowering the 
angle-of-attack to the nominal cruise value and adding power as neces- 
sary to offset drag. The cruise condition is maintained by modulating 
primarily engine thrust, but also angle-of-attack as necessary to main- 
tain a constant velocity. Throttle setting is commanded such that 
thrust equals sensed drag or 

THROTTLE - DRAG/ (THRUST x COS (a) ) 

where a is angle-of-attack. If thrust is insufficient to cancel drag 
(e.g., single-engine failure), then throttle is fixed at 100% and a is 
adjusted in 0.5 degree increments until thrust is able to offset drag. 
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5. 2. 2.1 HEAT-RATE CONTROL 


With velocity fixed, altitude and hence heat-rate are controlled by 
modulating bank angle. Heat-rate is again defined by Chapman's 
equation. Heat-rate control is achieved by flying the vehicle at a den- 
sity appropriate to the current velocity. The desired density is then 

Po = p SL QDOT 0 2 / ( (1 7600) (V r /26000) 3 15 ) 2 

where the "D" subscripts represent desired values. So, given a velocity 
and desired heat-rate, a desired density can be calculated and converted 
to a desired density attitude (h p ) with a simple exponential model. The 
heat-rate controller equation is a standard second order controller 
defining a desired vertical acceleration (a v ) such that 

a v " k 1 Ah p + k 2 Ah p ‘ + * h p" 

where k, and k 2 are constants. Ah p , Ah p ', and Ah p " are the difference 
between desired and current density altitude, density altitude rate, and 
density altitude acceleration respectively. Desired altitude rate and 
altitude acceleration are determined by back differencing. The desired 
altitude acceleration is smoothed with a low pass filter. The current 
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altitude rate is the vertical component of the vehicle velocity and this 
is back differenced to determine current altitude acceleration. 

Finally the desired vertical acceleration can be converted to a com- 
manded bank angle (4> c ) such that 

*0 = Vev ■ a v (1 80 /tt) / (TAN (4> ppev ) ) a sv 

where <l> prev is the previous bank angle and a $v is the vertical component 
of the sensed nongravi tat i ona 1 acceleration (19). 

5. 2. 2. 2 TIME-TQ-EXIT 

While the heat-rate controller is maintaining a constant heat-rate, 
a predictor is run every ten seconds to see how much of the desired 
plane change could be accomplished using nominal exit controls and con- 
ditions if the exit was begun at current time. Specifically, the guid- 
ance system integrates the vehicle state ahead to 400,000 feet assuming 
nominal alpha, bank and 100% thrust to determine the final inclination. 
When the predicted inclination angle (i p ) is greater that the desired 
inclination angle (i D ) then exit is started. There is an additional 
constraint on determining the time to begin the exit. For heat control, 
any bank angle greater than the current commanded bank angle will cause 
overheating of the vehicle. Thus when determining the time to begin the 
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exit, the algorithm uses the smaller of current or nominal exit bank 
angle in the integration to 400,000 feet. 

5.2.3 EXIT PHASE 

The exit phase begins when exit is commanded by the time-to-exit 
predictor and ends when the vehicles leaves the atmosphere at 400,000 
feet. During this phase, thrust is increased to 100%, ang 1 e-of-attack 
is set at the nominal exit value, bank angle is modulated to achieve the 
desired final inclination, and final apogee altitude is controlled by 
engine cutoff. However, heat-rate control must be continued until the 
vehicle climbs to an altitude where heating is no longer a concern. 
Therefore, during the first minutes of exit, if heat control requires a 
smaller bank angle than inclination control, heat control takes priori- 
ty. Once heat-rate commands a bank angle greater than inclination con- 
trol, then inclination control takes priority (19). 

5. 2. 3.1 INCLINATION CONTROL 

During the final part of the exit phase, the inclination at 400,000 
feet is controlled by modulating bank angle, A predictor/corrector gui- 
dance routine predicts a bank angle that, if held constantly throughout 
the remaining exit phase, will achieve the desired inclination angle. 
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The predicted bank angle is determined by finding a linear relation- 
ship between the final inclination and bank angle. The predictor algo- 
rithm numerically integrates the trajectroy ahead to 400,000 feet to 
determine the final inclination (i old ) . Another intergration is made 
with a slightly different bank profile (<l> prev + &4>) and the new inclina- 
tion (i new ) is determined. A partial derivative of inclination with 

respect to bank angle (di/d<t>) is calculated such that 

9'/9* = ('new " 'oJ /8<i> 

With the sensitivity of changes in inclination to changes in bank angle, 
the change in bank angle (5<J> 0 ) to achieve the final inclination (i D ) can 
be determined as 

* ('d ~ i 0 icP/ (9i /d&) 

Then the commanded bank angle is the previous bank angle plus the update 
or 


+ 

prev 


5<i>p 


If the update is large, then non-linear effects may require the process 
to be iterated. 
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5. 2. 3. 2 APOGEE ALTITUDE CONTROL 


The final apogee altitude is controlled by engine cutoff time. The 
engine cutoff is commanded when the vehicle has sufficient velocity to 
coast to the desired apogee altitude (ha D ) . For simplicity, the initial 
apogee prediction is done assuming the vehicle is moving conically 
(point mass gravity and no atmosphere). The expression for the apogee 
of any conic orbit can be written in closed form. However, since the 
engine cutoff occurs in the atmosphere, the conic apogee (ha c ) cutoff 
altitude is biased to compensate for the residual atmospheric drag 
between cutoff and apogee. So, the cutoff time is when 

ha c = ha 0 + ha B 

where ha B is apogee bias. To compensate for perturbations, this bias 
must be adjusted using the predictor/corrector's estimation of apogee 
altitude (ha p ) at 400,000 feet such that 

(ha) € = ha p - ha D 

ha B = ha B + K (ha) e 

where (ha) t is apogee error and k is a constant. This process is 

repeated if apogee error is greater than five nautical miles. 


118 


5.3 GUIDANCE RATIONALE 


5.3.1 ENTRY PHASE 

The reasoning for going with constant attitude or open-loop is based 
on the overall results of Chapter 4. Specifically, the results in 
Tables 4-8, 4-12 and 4-13, where the entry parameters were fixed, indi- 
cated that the cost was at most 100 pounds of fuel. In actuality, the 
possibly of achieving a 100 pound savings by adjusting entry parameters 
was remote since the optimization program achieved this result with 
apriori knowledge. In fact, we could visualize situations where a cor- 
rection for one deviation is made on the front-end while it might be 
reversed on the tail-end of the trajectory. At the same time, allowing 
the errors to build up and trying to fix them at the end did not seem 
reasonable either. First, it goes completely against what our trajecto- 
ry analysis indicated was the most fuel efficient method, and second, it 
might over tax the guidance system to try and compensate too late. 
Thus, we compromised and chose to begin closed loop guidance at the 
start of the cruise phase. The final consideration during entry was the 
lead-point for the heat-rate controller. Once again a compromise 
approach was chosen. To avoid exceeding the maximum heat-rate, an ade- 
quate lead-point for the heat-rate is needed. However, as shown in 
Chapter 4, if the lead-point is too low, it is possible to cutoff the 
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first heat-rate peak which can cost up to 450 pounds of fuel. Based on 
test results, 90% of the maximum heat-rate was determined to be a rea- 
sonable compromise (19). 

5.3.2 CRUISE PHASE 

The rationale for control of heat-rate and time-to-exit are a result 
of the baseline plane change strategy and our trajectory analysis. The 
optimization analysis showed that fuel could be minimized by doing the 
pullout and cruise at maximum heat-rate. Also, the time in the atmos- 
phere could be reduced by doing as much of the plane change during the 
exit as possible. 

5. 3. 2.1 HEAT-RATE CONTROL 

The control of heat-rate during the cruise phase is very important. 
Overheating will jeopardize vehicle safety and underheating will cost 
additional fuel. In order to maintain tight control over heat-rate, 
perturbation guidance was chosen. 

Heat-rate, based on Chapman's equation, is only a function of densi- 
ty and velocity. Therefore, of the available control parameters dis- 
cussed in Chapter 3, only angl e-of-attack, bank angle and engine thrust 
are possible controller candidates. Engine thrust has no direct effect 
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on density so it could only be paired with velocity. Bank angle and 
angle-of-attack can be used to control density or velocity. Analysis of 
the POST results indicated that bank angle was always able to control 
the heat-rate and hence density but the nominal angle-of-attack along 
with maximum thrust (e.g. f single-engine failure) was not always able to 
maintain a constant velocity. ■ Therefore, it was reasonable to use 
angle-of-attack and engine thrust to maintain velocity. Because of the 
multiple influence of angle-of-attack (affects CL, CD and the component 
of thrust aligned with drag) and for simplicity it seemed more appropri- 
ate to use angle-of-attack only as necessary and primarily modulate 
thrust to maintain velocity. Therefore, the final scheme sets a nominal 
angle-of-attack and adjusts thrust as necessary to maintain a constant 
velocity. If thrust is insufficient to cancel drag (i.e., maintain 
velocity), then angle-of-attack is adjusted until drag is offset. 
Recall that this reduction of alpha is beneficial for two reasons: (1) 

The CD is lowered as angle-of-attack is reduced, thereby reducing drag 
and (2) A reduced alpha increases the component of thrust in the direc- 
tion of drag. Separately, bank angle is modulated to control altitude 
which controls density which in turn controls heat-rate. 

5. 3. 2. 2 TIME-TD-EXIT 

The rationale for the time-to-exit predictor is a result of the tra- 
jectory analysis work in Chapter 4 and an earlier study by R. Richards 
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(19). The trajectory analysis results generally showed that fuel is 
minimized for the entire trajectory by adjusting the exit point and 
using near nominal parameters during exit. The study that Richards con- 
ducted showed that fuel was minimized for the exit phase when the ERV 
started the exit at a point where it could just make the end-conditions 
using approximately 90 degrees of bank and nearly nominal angle-of-at- 
tack. However, this study assumed no perturbations were encountered 
during exit and left very little room to compensate for any future per- 
turbations. Recall from Chapter 4, optimum exit bank angles for all 
cases considered with POST ranged from 65-72 degrees. Thus, the scheme 
of exiting when the predictor concluded you could make the end-condi- 
tions with nominal exit controls and conditions seemed to offer a near 
optimum solution while providing some extra margin to compensate for 
future perturbations. 

5.3.3 EXIT PHASE 

The rationale for the end-condition controller during the exit phaes 
is a result of the discussion in Chapter 3 on control variable pairings 
and analysis of the POST results covered in Chapter 4. For simplicity, 
we decided to make a one-to-one pairing of each control variable and 
targeting parameter and if possible, handle each pair separately even 
though there are coupling effects. Likewise, even though there is one 
more control variables than items to constrain, we chose not to go with 
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the added complexity and convergence concerns unless the extra variable 
was needed or offered some clear advantage. Recall from Chapter 4, 
ang 1 e-of -attack was hardly used during exit and consistently used less 
than bank angle for fuel efficiency reasons. Thus, we chose to use bank 
angle only as an inclination controller and engine cutoff for apogee 
control. Hence, angle-of-attack is left as an additional control vari- 
able, completely free for future guidance requirements. 
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CHAPTER 6 


GUIDANCE ALGORITHM TESTING 


6.1 INTRODUCTION 

To evaluate the performance of the plane change guidance scheme dis- 
cussed in Chapter 5, many simulations are made using various operating 
conditions. Although the guidance scheme is the same, the character- 
istic trajectory for a heating rate of 80 as compared to 125 BTU/ft 2 -sec 
is quite different. Therefore, they are treated as separate cases for 
the various operating conditions. 

The entire plane change process entails using the Earth's atmosphere 
to change the orbital plane. Therefore, if the atmosphere is not a nom- 
inal atmosphere with standard density values, the guidance algorithm 
must still be able to do its job. Many test simulations are made with 
off-nominal atmospheric conditions to verify that the algorithm does, in 
fact, perform as desired. Among these tests are density bias tests, 
where the nominal density is increased or decreased by a constant fac- 
tor, and density shears, where the trajectory passes through a varying 
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density profile model. Various values of density biases and different 
shear models are simulated to evaluate the algorithm's performance. 

Another test of the algorithm is to change the aerodynamic charac- 
teristics of the ERV. This is done by using various combinations of CL 
and CD. A final test was made to determine the algorithm's response to 
a single-engine failure. This was done by eliminating the thrust pro- 
duced by one of the three ERV engines. 

After testing each off-nominal condition individually, several com- 
binations are tested as well. For example, a thin atmosphere is tested 
in combination with low aerodynamic characteristics. Some of the combi- 
nations of off-nominal conditions yield a more realistic simulation ove- 
rall, as there are bound to be multiple variations in a real environment 
at all times. However, some of the combinations probably provide a more 
severe test of the algorithm than would be reasonable to expect the ERV 
to encounter. The scheme was not specifically designed to handle these 
cases but are included to demonstrate its robustness. 

6.2 EVALUATION CRITERIA AND THE NOMINAL TRAJECTORY SIMULATION 

There are several quantities that are useful in evaluating the per- 
formance of the guidance algorithm. The main goal of the atmospheric 
maneuvering is to change the inclination of the orbital plane. If the 
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algorithm does not guide the ERV to the desired inclination, it is not 
doing the job. Therefore, one of the most important quantities is the 
amount of inclination change that was accomplished during the atmospher- 
ic maneuver. The algorithm has been designed for a nominal plane change 
of 20 degrees. If the actual plane change in a test run is within 0.1 
degrees of that nominal value, the test will be considered successful. 
The closer the inclination is to the nominal value, the better, as less 
fuel will be required to correct the orbit. 

Another important criterion is the final apogee altitude. The algo- 
rithm has been designed for a nominal target apogee of 160 nautical 
miles. If the actual apogee in a test run is within 5 nautical miles of 
that nominal value, the test will be considered successful. Once again, 
the closer the apogee is to the nominal value, the better as less fuel 
will be required to correct the apogee. 

Since the potential fuel savings is one of the main reasons for per- 
forming the maneuver, fuel use is important and is also a consideration 
in the performance evaluation. For each test run, the fuel required to 
circularize the orbit at the target apogee height is computed assuming 
that engine burns are made at the apogee and then at the raised perigee. 
Also included in this calculation is the fuel required to correct any 
inclination error. Comparing the total fuel required to obtain the 160 


nm circular orbit to the nominal trajectory simulation run provides one 
measure of the algorithm performance. 

In addition, a comparison with the fuel used in the optimized tra- 
jectory generated by POST for the same conditions provides an evaluation 
of how well the algorithm performed in comparison to the best it could 
be expected to perform. The only caveat being, the fuel results from 
POST shown in Chapter 4 and the summary of results for this chapter are 
for a 1976 standard atmosphere. The CSDL simulations used a 1962 stand- 
ard atmosphere. These models are very similar below 290,000 feet. 
Therefore, the comparison between POST and the CSDL simulation fuel 
results while not exactly comparable should offer some valid insight. 

Also computed is the fuel required to correct the perturbed orbit's 
longitude of the ascending node (LAN) to the nominal. Although no 
attempt was made to control the node shift, the node does shift during 
the atmosphere phase of the trajectory. This is because the inclination 
change was spread out across the node, and unless it was balanced evenly 
before and after the node there is going to be some node shift. For the 
CSDL nominal 80 and 125 BTU/ft 2 -sec trajectories the LAN shifted from 
zero to 3.91 and 1.44 degrees respectively. Therefore, all simulations 
will use this LAN as the final desired LAN. While the fuel to correct 
the note shift is typically small, a large value will indicate that the 
algorithm's performance could possibly be improved. As a point of ref- 
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erence, the final LAN for the POST 80 and 125 nominal trajectories (1962 
atmosphere) was 3.90 and 1.44 degrees respectively. 

The final important criterion used in the performance evaluation is 
maximum heating rate. Since there is a limit to the heating rate that 
the ERV can withstand, a run which allows the heating rate to exceed its 
limit has safety of flight problems. On the other hand, if the algo- 
rithm allows the vehicle to operated below it maximum limit, fuel effi- 
ciency will be hurt because of the correlation between heat-rate and 
plane change performance (6). Although exceeding the maximum heating 
rate limit would probably cause a much worse consequence than degraded 
plane change efficiency, for simplicity, an actual heating rate within 5 
BTU/ft 2 -sec of the desired heating rate, high or low, will be considered 
acceptable and the test successful. As mentioned in Chapter 2, total 
heat load for this type of plane change is not considered a constraint 
and is not included as a performance criterion. 

The CSDL 80 and 125 BTU/ft 2 -sec nominal trajectory simulation runs 
will be used as the primary basis for comparison. For these runs the 
inclination change is 20.007 and 20.006 degrees respectively, and the 
final apogee altitude is 158.0 and 158.5 nautical miles. The final 
vehicle weight in a 160 nm circular orbit with the complete plane change 
is 5450 and 5961 pounds respectively, and the maximum heating rate is 80 
and 126 BTU/ft 2 -sec. For reference, Figures 6-1 thru 8 present plot 


profiles for the 80 and 125 nominal synergetic plane change trajectories 
generated by the CSDL computer simulation using the plane change guid- 
ance algorithm. These figures should be useful for comparison with the 
figures of the various test runs presented in the following sections. 

6.3 SIMULATION TEST RESULTS 

A wide variety of off-nominal conditions are simulated to evaluate 
the guidance algorithm's performance. An example of the resulting tra- 
jectory plots will be presented for each type of test condition and 
where some new information is contained in the plotted data. Tables 
will summarize the test results in each subsection that follows. 

6.3.1 SINGLE OFF-NOMINAL CONDITION TESTS 

The following subsections will present the results for tests which 
involve one off-nominal condition while all other values are nominal. 
The tests include density biases, density shears, aerodynamic uncertain- 
ties, and single-engine failure performance. 

6. 3. 1.1 DENSITY BIASES 

The nominal density levels are multiplied by a constant factor in 
this set of tests to determine how thin or thick an atmosphere the guid- 
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Figure 6-1. Altitude profiles for QDOT = 80 and 125, nominal conditions 
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Figure 6-4. Bank angle profiles for QDOT = 80 and 125, nominal conditions. 
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Figure 6-6. Heat-rate profiles for QDOT = 80 and 125, nominal conditions. 
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Figure 6-7. Inclination change profiles for QDOT = 80 and 125, nominal conditions. 
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Figure 6-8. Apogee and perigee altitude profiles for QDOT = 80 and 125, nominal conditions. 





ance algorithm can handle. A -30% atmosphere is one which is constantly 
30% below normal density, while a +30% atmosphere is 30% more dense, or 
thicker, that a normal atmosphere. When the atmosphere is too thin, the 
ERV does not turn as fast, so the plane change tends to take longer and 
use more fuel. On the other hand, too thick an atmosphere turns the 
vehicle quicker and uses slightly less fuel. Overall, the algorithm 
performs excellently in atmospheres ranging from 50% of normal to 150% 
of normal. Table 6-1 and 2 contains the results of these test runs. 
Figures 6-9 thru 12 show the plot profiles for the QDOT = 80 trajectory 
with a +50% density bias. 

From these test results two points should be noted. First, a com- 
parison of the final vehicle weights for the CSDL and POST simulations 
suggests that the designed guidance algorithm performs as well or better 
(e.g., QDOT ° 80) than POST. Recall however, from earlier in this chap- 
ter, that the CSDL and POST trajectories were generated using different 
atmospheric models. For the 125 trajectories, the vehicle weights com- 
pare very close. For the 80 BTU/ft 2 -sec trajectories there is a signif- 
icant difference in vehicle weights (about 200 lbs). It should be 
remembered that the 80 trajectory stays significantly higher in the 
atmosphere (approximately 35,000 feet), and the difference between the 
atmospheric models is more significant for the higher altitudes. There- 
fore, the comparison between CSDL and POST vehicle weights may not be as 
valid for the 80 trajectories. 
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Table 6-1. Density bias test results for QDOT = 80 Btu/ft^-sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 






































Table 6-2. Density bias test results for QDOT = 125 Btu/ft^-sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 
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Second, the test results for both the 80 and 125 BTU/ft 2 -sec trajec- 
tories indicate that as the density bias increases so does the shift in 
the longitude of the ascending node. The 125 trajectory shows the larg- 
er node shift reflecting its higher sensitivity to density perturba- 
tions. Still, the node shift is relatively small considering the size 
of the tested density biases. 
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6. 3. 1.2 DENSITY SHEARS 


Density shears are introduced along the plane change trajectory to 
evaluate the guidance algorithm's performace flying through regions of 
changing density factors or levels. Two types of density shear models 
are used for testing, the density shear experienced by the shuttle on 
the STS-4 mission and a derived model. The STS-4 model is shown in Fig- 
ure 6-13 and the derived model in Figure 6-14 as a normalized density 
verses altitude plot. The density is normalized by the 1962 standard 
atmospheric density. The shear models are mechanized as a table of den- 
sity multipliers verses altitude. 

The algorithm has no difficulty with density shears. The density 
shears encountered during the descending portion of the trajectory, are 
simply flown through, as the ERV approaches level-off, the automatic 
compensation of heating rate that was discussed in Chapter 4 tends to 
automatially adjust the cruise altitude to correct any effects caused by 
the density variations. Even when the shears are encountered during 
exit, the bank controller is easily able to correct for the density var- 
iations. 

In fact, the density shears were so benign, a ±30% density bias was 
added to the STS-4 shear during cruise. The objective was to show that 
the guiance would respond to the changing density and guide the ERV back 
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to the reference heating rate. The density bias was added as the vehi- 
cle flew below 225,000 feet. This caused a step change in density and 
heating rate that does not reflect real-world conditions. However, the 
guidance algorithm very rapidly drove the ERV back to the reference 
heating rate. Table 6-3 contains the results of a sample of the tested 
cases. The heating rate values in the parentheses are for cruise equi- 
librium conditions and are a better reflection of the algorithm's per- 
formance for the shear tests with the added density bias. Figures 6-15 
thru 18 show profile plots for a 125 trajectory with a STS-4 density 
shear prof i 1 e. 

6. 3. 1.3 AERODYNAMIC UNCERTAINTIES 

This test is useful for determining how well the guidance algorithm 
can perform if the vehicle's modeled aerodynamic characteristics are in 
error. This error is simulated by inputing a straight bias for CL and 
CD to be used throughout the plane change trajectory. Guidance was 
never given any knowledge of the CL and CD shift. While there are an 
unlimited possible combinations of CL and CD errors, the testing was 
limited to the same four combinations of ±15% CL and CD used in Chapter 
4, plus two extreme cases of ±25% CL and CD error to demonstrate the 
robustness of the guidance system. It should be noted that the 15% 
cases included two cases with no net change of L/D and two cases with a 
L/D change (one increased and one decreased). 
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Table 6-3. Density shear test results for QDOT = 125 Btu/ft^-sec. 


Maximum 

QDOT 

(Btu/ft^-sec) 

lO 

CM 

f-H 

126 

130 

(125) 

128 

(125) 

126 

(125) 


cu 











c 











ft} 










O r— 

t — • 










+■> <o 

Q- ^ 










■r— 

cn 










r— +-> 

r— -O 

o 

oo 


00 


<3- 


00 


cu t_ 

ftJ r— 


f-H 




r»s. 


f-H 


3 (U 

C ^ 










U_ C 

• r— 










t — * 

E 











o 











Z 






















Z CT> 


LO 


o 


CM 


LO 



*=C QJ 


cn 


CM 


oo 


CO 



__j -a 

• 

• 


• 


• 


• 




iH 

f-H 


»— H 


f-H 


f-H 


cu 

• 










r — 

_a 










u 

o 










c 











XI *r- 

• 










CU 

— 

oo 









=> 4-> 

•r— tn 


<c 


<=£ 


<c 


■=£ 


-C 

O _Q 

oo 









= cn 

r— 

lO 

z 


z 


z 


z 


H— *i— 

E 










CO CU 

C 










O 3 











Q. 

o 










z 












r— 1 










+-> 

4-> 










-C 

•r— 










cn 

-Q 










■r- E 











Q) C 

o 










3 

x — ». 

f-H 

VO 


CM 




cn 


O 

£_ cn 

CO 

r^. 






oo 


<U VO 

03 X) 

<T» 



CT» 


oo 


cn 


r - rH 

r— r—— 

LO 

LO 


LO 


LO 


LO 


u 

3 ' 










•1- c 

CJ 










-C -r- 

C- 










CU 

•r— 










> 

O 











CU 










CU 

“O 










cu 

3 ^ 

LO 



LO 


r— H 


CM 


CT)+J E 

• 

• 


• 


• 


• 


O 

•p- c 

CO 

vo 


CT> 


o 


f-H 


CL-t-. — ' 

LO 

LO 


LO 


VO 


VO 



r— — 

i - H 

rH 









c 










c 











o 






















4-> 

cu ^ 

CO 



LO 


<2> 


o 


03 

cn cn 

O 

o 


o> 


o 


o 


c 

c cu 

O 

o 


cn 


O 


o 


•r— 

o3 -a 

• 

• 


• 


• 


• 


r— 

^ — 

o 

o 


<T> 


o 


o 


u 

o 

C\J 

CM 




CM 


CM 


c 











t— * 
















cn 


cn 


>> 





t- 

C. 

rd 

c_ 

03 


+j cn 



c 


*3 

<T3 

•r— 

<T3 

•i— 


x: m- c 



o 


CU 

CU 

^ -Q 

CU 

crS. -O 


+-> cn *r- 



•r— 


-C 

_c 

O 

-C 

O 


•r- C t- 



H-> 

r— 

to 

cn 

ro >, 

cn 

oo >, 


-a s qj 3 



•f— 

fO 



+ +j cn a> 


i +j cn 

CU 

CU "O "O 

cu 


-a 

c 


<1- 

•r- c cn 


•r- C 

cn 

> t- 

cn 


c 

•T— 

1 

i 

-C cn •«— -r- 

i 

_C cn -r- 

■i— 

•i— fT3 cn 

t— 


o 

E 

CO 

CO 

•P C L d 

co 

4-> C *- 

3 

L (DO (3 

3 


o 

o 

#■— 

H- 

T- <u 3 u 

J— 

•P- cu 3 

C- 

(Uxcn-r- 

£_ 



z 

CO 

CO 

JT3T3 a 

CO 

s -o XI 

O 

O l/> + _Q 

O 


150 


NOTE: 1. LAN - Longitude of the ascending node. 

2. The heat-rate values in the parentheses are for cruise equilibrium conditions. 

3. N/A - Not available. 
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Figure 6-16. Angle-of-attack and bank angle profiles for QDOT = 125, 
STS-4 density profile. 
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Overall, the algorithm handled the uncertainties well and exhibited 
no unexpected performance loss/gain. Tables 6-4 and 5 summarize these 
results and offer a comparison with the performance loss/gain that was 
seen with POST. Figures 6-19 thru 22 show the plot profiles for a QDOT 
= 125 test with -25% CL and +25% CD. 

Several points are worth noting: 

1. The fuel usage for both the 80 and 125 trajectories was very 
similar to that calculated by POST under the same conditions. 

2. For the performance cases with no change in L/D (e.g., +CL and 
+CD) , the final vehicle weight is almost the same as the nominal 
vehicle weight for both the QDOT 80 and 125 trajectories. However, 
in terms of LAN shift, the results were opposite. The QDOT = 80 
trajectory had a large node shift while the 125 trajectory showed 
almost no shift. 

3. For the cases with a change in L/D (e.g., +CL and -CD), the 
final vehicle weight reflected that change with a corresponding 
increase or decrease for both QDOT = 80 and 125. However, in terms 
of LAN shift, the results were again opposite. The 80 trajectories 
showed almost no node shift while the 125 trajectories showed a 
large node shift. 
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Table 6-4. Aerodynamic uncertainty test results for QDOT = 80 Btu/ft^-sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 
































Table 6-5. Aerodynamic uncertainty test results for QDOT = 125 Btu/ft 2 -sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 
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Figure 6-22. Inclination change, apogee and perigee altitude profiles 
for QDOT = 125, -25% CL and +25% CD. 










4. The 25% low vehicle performance case (QDOT = 125) lost the char- 
acteristic "double-peak" heat-rate curve as shown in Figure 6-21. 
This test and one other multiple off-nominal test were the only 
cases that lost the "double-peak" suggesting that the 90% of maximum 
heat-rate was an adequate lead-point. 

5. The 15% low vehicle performance case (QDOT = 80) failed to meet 
the apogee altitude criterion. Recall from Chapter 4 that for a 
heat-rate of 80 BTU/ft 2 -sec engine cutoff occurs much lower in the 
atmosphere than for a heat-rate of 125. This makes adequate compen- 
sation for large perturbations more difficult for a QDOT = 80 tra- 
jectory. Overall, the apogee error does not represent a significant 
fault in the guiidance algorithm but some improvement might be need- 
ed . 


6. 3. 1.4 SINGLE-ENGINE FAILURE 

In this test, the algorithm's ability maintain control of the vehi- 
cle, complete the plane change and return safely to orbit is determined. 
Two tests are conducted with single-engine failure at a different 
location for each test. We felt the start of the level-off phase and 
the start of the exit phase represented the two worst possibi 1 ties. 
Guidance was given knowledge of the engine failure in both cases (i.e., 
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the predictor/corrector used a reduced thrust model). Overall, the 
algorithm was able to compensate for the engine failure in both cases 
and still obtain the desired plane change and final apogee altitude. 
Table 6-6 summarizes the results while Figures 6-23 thru 26 show the 
profile plots for QDOT = 125 with single-engine failure at level-off. 

Three points are worth noting: 

1. The CSOL piane change algorithm was able to achieve almost the same 
fuel results as the POST simulation. 

2. A comparison of Figure 6-1 and 6-23 shows the longer mission time 
for the QDOT = 125 single-engine case. The increased time to complete 
the mission caused a larger LAN shift than with the nominal mission. 
Thus, a significant fuel penalty would be paid to achieve the same 
orbital plane as the nominal run. 

3. The QDOT ■ 80 case was almost totally unaffected by the single-en- 
gine failure. Although not shown, the mission time was nearly nominal 
and as the results show in Table 6-6 , very little LAN shift occurred. 


6.3.2 MULTIPLE OFF-NOMINAL CONDITION TEST 
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Figure 6-24. Angle-of-attack and bank angle profiles for QDOT = 125, 
single engine failure at level-off. 















600 3WddJ 



010 3Wd«3 


168 


Figure 6-26. Inclination change, apogee and perigee altitude profiles 
for QDOT = 125, single-engine failure at level-off. 







Table 6-7. Multiple off-nominal test results for QDOT = 80 Btu/ft^-sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 










































Table 6-8. Multiple off-nominal test results for QDOT = 125 Btu/ft 2 -sec. 
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NOTE: 1. LAN - Longitude of the ascending node. 

2. N/A - Not available. 
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VAR 06 = CURRENT (DEC) 



fioo 3WdUj 






03 


<U 

A 

u 

LO 

3 

CVJ 

r— 

r“H 

•n 


03 

II 

4- 

1— 

0J 

o 

c 

Q 

•r* 

O' 

cn 


c 


0J 

o 


4- 

<D 


r— 

lO 

cn 

CD 

c 

r— 

•n 

•I— 

in 

4— 


O 

* 

C- 

o 

Q. 

o 

CD 


r— 

Lf) 

cn *-h 

c 

+ 

03 



T3 


c 

C 

03 

03 


-Q 

— 1 


O 

“O 


c 

** 

03 

LO 


i— 1 


1 

O 


03 

is 

4-> 

>> 

4-> 


03 

•r- • 

1 

CO 4- 

4- 

C 4- 

O 

<u o 

1 

"O 1 

CD 

r— 

f— 

^ a) 

cn 

o > 

c 

CO CD 

< 

1 •— 


• 

00 

CVJ 

I 

co 

<D 

t- 

3 

cn 


900 3Hdyj 



172 






STfIGNRT i ON HERTING RRTE 



a; c 

JC 03 

ioo 3woyj 

-a — i 
c o 

03 



SOO 3WUUJ v 


173 
















By combining two or more of the off-nominal conditions discussed in 
Section 6.3.1, further tests of the algorithm's performance can be made. 
The purpose here is to explore the combinations which yield useful 
information in the determination of what the guidance algorithm can and 
cannot do. The results of all of the combinations tested are summarized 
in Table 6-7 and 8. 

Several points are worth noting: 

1. All the test cases that had some combination of low vehicle perform- 
ance and single-engine failure (e.g., -30% density, -15% CL and +15% CD, 
and single-engine failure) experienced bank control saturation during 
level-off. Figures 6-27 thru 30 show the profile plots for one such 
case. The "common thread" for all cases was the significant change in 
cruise angl e-of-attack . As a result, bank angle was saturated for about 
30 seconds and a heat-rate overshoot occurred. Still, the heat-rate 
remained within the performance criterion. Delaying the angl e-of-attack 
adjustment by 30 or 40 seconds after level-off would probably eliminate 
this bank saturation. 

2. The same test cases that experienced bank saturation also lost 
velocity during the initial part of the exit phase. However, this loss 
of velocity had no significant overall effect as the vehicle was able to 
complete the intended maneuver. 
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3. 


Four of the five test combination cases for a heat-rate of 80 


BTU/ft 2 -sec failed to meet the apogee altitude criterion. However, only 
two of these cases (+30% density with 15% low aero and -30% density with 
15% low aero) were significantly outside the criterion. Inaddition, the 
combination of off-nominal conditions may be too severe of a test (i.e., 
unrealistic). Still, as discussed earlier, these results point out a 
need for improvement in the apogee altitude predictor. 

4. Table 6-8 shows again how comparible the CSDL guidance system 
results (in terms of fuel) are to the POST results. 

5. Several of the QDOT = 125 test showed a significant shift in the 
LAN. For the combinations involving high aero, the final vehicle weight 
to achieve the nominal inertial plane is almost nominal. However, for 
the combinations with low aero, the vehicle weight to achieve the nomi- 
nalm orbit was significantly below nominal. 
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7.2 CONCLUSIONS 


1. A guidance system has been developed for the ERV synergetic 
plane change maneuver. 

2. The recommended plane change algorithm is a hybrid of perturbu- 
tion guidance and pred i ctor/cor rector guidance. Perturbation guidance 
is used to control heat-rate during cruise, and the predictor/corrector 
guidance is used to control the end-conditions. 

3. The guidance system explicitly controls the end-conditions and 
will guide to any specified heat-rate, inclination and apogee combina- 
tion. 


4. The guidance system provides accurate control over heating rate 
within its knowledge of density and velocity. Due to the importance of 
heat-rate control (both from a performance and a vehicle safety point of 
view), an accurate, reliable measurement/estimate of density and veloci- 
ty should receive top priority. 

5. The guidance system has demonstrated the ability to handle ±30% 
density biases, two density shear profiles, ±15% aerodynamic uncertain- 
ties, a single-engine failure, and various combinations of the above. 
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PRECEDING PAGE BUNK NOT FILMED 


The guidance system is adaptive in the sense that it will respond to 
these dispersions and maintain control over heat-rate and the desired 
targeting parameters. 

6. The aerodynamic uncertainties caused the largest uncontrollable 
performance effect (fuel usage) and therefore represent the largest 
uncertainty in efficiency of the ERV synergetic plane change maneuver. 

7. Simulated plane change missions using this guidance system have 
nearly the same total fuel usage as trajectories generated by an opti- 
mization program. 

8. Overall, with this guidance algorithm in its software the ERV 
can be expected to accurately maneuver through the earth's atmosphere 
for the programmed synergetic plane change and thereby demonstrate capa- 
bilities that will be desired on future operational vehicles. 



As with most all challenging tasks, time never seems available to 
explore all the possible options nor to deal with all the opportunities 
that surface during the investigation. Listed below are five such items 
which seem to warrant further study and/or work. 
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7.3.1 NODE CONTROL 


While the objective of this task was to provide a guidance system to 
handle an ERV during a plane change, the ultimate desire will surely 
involve a system to maneuver into a specific inertial plane. Therefore, 
in addition to inclination control, control of the longitude of the 
ascending node would also be required. While we have addressed the node 
shift in Chapter 6 by including the fuel required to correct any unde- 
sired shift, it is not clear that both inclination and node can not be 
efficiently controlled during the atmospheric phase. In other words, 
would the total fuel for handling both in the atmosphere be less than 
the current scheme where only inclination is actively controlled in the 
atmosphere and any node shift is corrected on-orbit. 

The current algorithm has available three possible control vari- 
ables, (1) the length of the cruise phase, (2) angl e-of-attack during 
exit, and (3) engine throttling during exit. These would seem to offer 
a means of controlling the node. However, it is also unclear if the 
added complexity is worth any possible fuel savings. 

7.3.2 PLANE. CHANGE T RAJE CTO RY.. IMP R O VE MEN TS 

During our analysis of single-engine failure performance, some of 
our results seemed to suggest that some performance gain could be made 
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by allowing the thrust to vary during the exit phase. In theory, this 
concept could be tested with only two or three POST runs. Likewise, 
there are other possible parameters that have been fixed during our 
plane change analysis that if allowed to vary or change more often could 
offer some performance improvements. For example, a few POST runs were 
made where alpha and bank angle were allowed to change two additional 
times during entry, resulting in a small gain of roughly 100 pounds of 
fuel. Thus, it seems possible that a "more optimum" trajectory exists. 

7.3.3 HE A T - RATE CONTROL 

Our final scheme used a lead-point for our heat-rate controller that 
will still modify the first heat-rate peak under some perturbed condi- 
tions causing fuel efficiency to suffer. A method to avoid this was 
mentioned in Chapter 4 but never tested. Specifically it might be pos- 
sible to actively control heat-rate such that the first heat-rate peak 
is not eliminated or significantly modified. The resulting fuel savings 
could be as high as 8% of the total fuel used or roughly 450 pounds for 
a 20 degree plane change. 

7.3.4 MINIMUM FUEL 

One advantage of the predictor/corrector guidance scheme is the 
availability of fairly accurate predictions of end-conditions like fuel 
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remaining. This information could be used to make adjustments to the 
intended trajectory (i.e., reduced the amount of desired plane change) 
as necessary to insure the ERV returns to orbit with the minimum 
requi red fuel . 

7.3.5 APOGEE ALTITUDE CONTROL 

Several of the test cases in Chapter 6 for multiple off-nominal con- 
ditions, failed to guide to the apogee altitude within the desired accu- 
racy criterion. All of the cases were for a heat-rate of 80 BTU/ft 2 -sec 
and most of them involved off-nominal aerodynamic conditions. It might 
be possible to improve the apogee accuracy by providing an improved 
estimate of the aerodynamic characteristics to the predictor/corrector 
rout i ne. 
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APPENDIX A 


POST 


A.1 INTRODUCTION 

This appendix provides the details on the primary tool that was used 
to conduct the plane change analysis. The formulation of the plane 
change problem involves many nonlinearities, beginning with the 
equations of motion and heat transfer and skin temperature computations. 
In this situation, the value of lineari 2 ed analysis is quite limited 
and, at times, even misleading. Therefore, most of the analysis and 
synthesis work was accomplished with the aid of digital simulations. 
Specifically, this appendix will provide details on the Program to Opti- 
mize Simulated Trajectories (POST) and an example of trajectory gener- 
ation for a nominal 20 degree plane change as an illustration of POST. 

A. 2 POST DESCRIPTION 

The following was taken almost verbatim from Reference 1. POST, a 
generalized point mass, discrete parameter targeting and optimization 
program, provides the capability to target and optimize point mass tra- 
jectories for a powered or unpowered vehicle operating near a rotating 
oblate planet. Developed for the vehicle analysis branch, space direc- 
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torate, NASA Langley Research Center by the Martin Marietta Corporation, 
it was intended to provide the capability to simulate and optimize 
ascent and reentry trajectories for shuttle type vehicles. However, the 
generality of the program also allows various other types of vehicles to 
be simulated. POST has been used successfully to solve a wide variety 
of atmospheric flight mechanics and orbital transfer problems. The gen- 
erality of the program is evidenced by its N-phase simulation capabili- 
ty, which features generalized planet and vehicle models. This flexible 
simulation capability is augmented by an efficient discrete parameter 
optimization capability that includes equality and inequality con- 
straints. POST was originally written in Fortran IV for the CDC 6000 
series computers. However, it is also operational on the IBM 370, and 
Uni vac 1108 computers. 

One of the key features of POST is an easy to use NAMELIST-type 
input procedure. This feature significantly reduces input deck set-up 
time (and costs) for studies that require the normal large amount of 
input data. In addition, the general applicability of POST is further 
enhanced by a general-purpose discrete parameter targeting and optimiza- 
tion capability. This capability can be used to solve a broad spectrum 
of problems related to the performance characteristics of aerospace 
vehicles. 
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The basic simulation flexibility is achieved by decomposing the tra- 
jectory into a logical sequence of simulation segments. These trajecto- 
ry segments, referred to as phases, enable the trajectory analyst to 
model both the physical and the nonphysical aspects of the simulation 
accurately and efficiently. By segmenting the mission into phases, each 
phase can be modeled and simulated in a manner most appropriate to that 
particular flight regime. For example, the planet model, the vehicle 
model, and the simulation options can be changed in any phase to be com- 
patible with the level of detail required in that phase. 

Every computational routine in the program can be categorized 
according to five basic functional elements. These elements are: The 
planet model, the vehicle model, the trajectory simulation model, the 
auxiliary calculations module, and the targeting and optimization mod- 
ule. The planet model is composed of an oblate spheroid model, a gravi- 
tational model, an atmosphere model, and a winds model. These models 
define the environment in which the vehicle operates. The vehicle model 
comprises mass properties, propulsion, aerodynamics and aeroheating and 
a navigation and guidance model. These models define the basic vehicle 
simulation characteristics. The trajectory simulation models are the 
event-sequencing module that controls the program cycling, table 
interpolation routines, and several standard numerical integration tech- 
niques. These models are used in numerically solving the translational 
and rotational equations of motion. The auxiliary calculations module 
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provides for a wide variety of output calculations. For example, conic 
parameters, range calculations, and tracking data are among the many 
output variables computed. The targeting and optimization module pro- 
vides a general discrete parameter iteration capability. The user can 
select the optimization variable, the dependent variables, and the inde- 
pendent variables from a list of more than 400 program variables. An 
accelerated projected gradient algorithm is used as the basic optimiza- 
tion technique. This algorithm is a combination of Rosen's projection 
method for nonlimear programming and Davidon's variable metric method 
for unconstra i nted optimization. In the targeting mode, the minimum 
norm algorithm is used to satisfy the trajectory constraints. The cost 
and constraint gradients required by these algorithms are computed as 
first differences calculated from perturbed trajectories. To reduce the 
costs of calculating numerical sensitivities, only that portion of the 
trajectory influenced by any particular independent variable is reinte- 
grated on the perturbed runs. this feature saves a significant amount 
of computer time when targeting and optimization is performed (1). A 
more complete description and addtional information can be found in the 
formulation, utilization, and programmers manuals (1), (2), and (3). 

A. 3 TRAJECTORY GENERATION EXAMPLE 

The baseline plane change trajectories were primarily developed by 
R.W. Powell and J.C. Naftel of NASA Langley Research Center. However, 
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the author with heip and assistance from Powell and Naftel made modifi- 
cations and some improvements to the POST baseline trajectory input deck 
for his specific application. The following 14 pages are presented as 
an example of a POST input deck necessary to produce the nominal 20 
degree plane change trajectory with a heating rate of 125 BTU/ft 2 -sec. 
The complete ERV CL and CD data is included as a part of the input deck. 
The data is set up in three types of tables corresponding to three cor- 
relation parameters: (1) altitude (ALTITO) , (2) viscous (VINV) , and (3) 
mach number (MACH). Each correlation parameter table is broken into a 
number of sections for CL and CD (e.g., ALTITO = 300,000' and 537,000'). 
Within each section the CL or CD values are catalogued for ang 1 es-of-at- 
tack 0 to 60 degrees (1 degree increments). 

Following the input code is a brief description of the input condi- 
tions, the various phases of the trajectory, and the independent vari- 
ables that were used while optimizing the amount of fuel required for 
the plane change along with the constraints or dependent variables that 
were used in generating the plane change trajectory with POST. Most of 
the notation used in this appendix is the same notation that is used by 
POST. However, where possible both the POST and the more generic nota- 
tion used in the main body of this thesis is shown. 
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P 0 S T - MMC 3D TARGETING AND OPTIMIZATION PROGRAM. REVISED 052181 , PROBLEM NO. 1 

COPYRIGHT 1981 MARTIN MARIETTA CORPORATION 

LSSEARCH -* 1-- 

PERT= 01 , . 1 ,3*. 01 , .00001 ,6* 01, " 2-- 

PCTCC- .01 , -- 3” 

SRCHM s 4 , -- 4-- 

0PT=1, -- 5-- 

0PTVAR=6HWE IGHT , -- 6-- 

0PTPH= 1000 , *• 7-- 

NDEPV=3 , -- 8-- 

DEPVR=6HINC .6HALTA , 6HXMAX4 , -- 9” 

DEPVAL=48.5,160,1.0, " 10” 

DEPTL*. 10,5, .01 . -- 11” 

DEPPH= 1000, 1000, 1000, ” 12” 

NINOV= 12 , ” 13” 

INDVR=5HCRITR , 5HCRITR , 6HALPPC 1 .6HBNKPC 1 , ” 14” 

GHALPPC1 , 6HDPVLS 1 , 5HCRITR , SHDBANK , 6HDALPHA .5HCRITR , 6HALPPC 1 , 6HBNKPC 1 , -- 15” 

INDPH=50, 55, 100, 100, 105, 105, 180, 180, 180, 200, 200, 200, -- 16” 

U= - 9. 27. 27, -64, 25, .004, 40. -66, 20, 170, 8, -50, ” 17-- 

$ -- 18-- 

LSGENDAT -- 18-- 

PRNC = 1 , -- 19" 

NPC(9) = 1 ,NPC(22) = 1 ,ETAPC( 1 )=0, -* 20” 

TITLE=OH*ERV CRUISING SYNERGISTIC PLANE CHANGE 0-D0T=125», -- 21-- 

EVENT - 1 , -- 22-- 

SRE 1 77 . 40 , ” 23” 

LREF = 25 . , ' -- 24-- 

FESN= 1000, ” 25” 

NPC( 2 ) =4 , -- 26-- 

WGTSG= 1 1430, -- 27-- 

NPC( 1 ) =3 , -- 28-- 

NPC( 1 5 ) — 1 , ” 29-- 

NPC(3)=5. -- 30- - 

NPC( 5 ) =0 , -- 31” 

NPC( 8 ) =0 , ” 32” 

NPC( 1 2 ) =0 , -- 33” 

NPC( 16 ) = 1 , -- 34-- 

ISPV = 295 . 0 , ” 35-- 

DT = 100 , -- 36-- 

IGUID(1 ) = 3,0,1, ” 37 — 

ALTP= 160, -- 38-- 

ALTA= 160 , -- 39-- 

PINC= 1000 , -- 40” 

NPC( 25 ) = 3 , -- 4 1 - - 

PRNT (91 ) *5HXMAX 1 , 5HXMIN2 , 5HXMAX3 , 5HXMIN4 , 5HXMAX4 , 5HBETAI , 6HTIMRF 1 , 6HVIDE AL , -- 42” 

6HVINV ,6HINC , 6HPST0P , ” 43-- 

M0NX-6HHEATRT ,6HALTIT0 , 6HDYNP , 6HETA , -- 44-- 

INC=28 . 5 , -- 45 — 

TRUAN=150. ,PGCLAT=-19 , ” 46-- 

BETA=0 , ” 47-- 

$ --48-- 

LSTBLMLT TVC1M=3 ,AE 1M = 3 ,CLM= 5 ,COM= . 5 ,CLHM= . 5 ,CDHM= . 5 , CLVM= . 5 ,CDVM= . 5 $ ” 48” 

LSTAB TABLE =5HTVC IT ,0 , 1 100$ ” 48-- 

LSTAB TABLE=6HAE IT ,0,1.4 $ -- 48-- 

L$T AB -- 48-- 

T ABLE S 6HCLHT ,2 .5HALPHA .6HALTIT0 ,51.2,1,1,1,1,1,1,1,1, ” 49-- 

3. 00000000 E +05, -- 50” 

0. , 1 .54130435E-02, 1 . 00000000 E+ 00 , 2 . 90433333E -02 , -- 51-- 

2. 00000000 E +00, 4 . 26 16 1832E -02 , 3 . 00000000 E +00, 5 . 6 14 15108E -02 , -- 52-- 

4.00000000E+00, 6 . 96 270748 E -02 , 5 . 00000000 E +00 , 8 . 30790323E -02 , -- 53-- 

6. 00000000 E +00, 1 .03075910E-01 , 7 . OOOOOOOOE+OO , 1 . 23070807E -01 , -- 54-- 

8 . OOOOOOOOE +00 , 1 .43063966E-01 , 9 . OOOOOOOOE+OO , 1 . 63055590E -01 , -- 55-- 

1 .OOOOOOOOE+OI , 1 .83045852E-01 . 1 . 10000000E+01 , 2 .09544373E-01 , -- 56-- 

1 .20000000 E +01 , 2 . 36086535 E -01 , 1 . 30000000E+01 , 2 . 62668672E -01 , -- 57-- 

1 .40000000 E +01 , 2 . 892875 18E -0 1 , 1 . 50000000E+01 , 3 . 15940152E-01 , -- 58-- 

1 .60000000E+01 . 3 46924729E-01 , 1 . 70000000 E +01 , 3 . 78013237E -01 , -- 59-- 

1 80000000 E +01 , 4 .09 199247E -01, 1 . 90000000E+01 , 4 . 40476852E -01 , -- 60” 

2 . OOOOOOOOE +01 , 4.71 8406 14E -0 1 , 2 . 10000000E+01 , 5 . 0592 1 147E -01 , -- 61-- 
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ORIGINAL PAGE r- 
OF POOR QUALITY 


2.20000000E+01 , 
2 . 40000000E +0 1 , 
2 . 60000000E+0 1 . 
2. 80000000 E +01 . 
3. 00000000 E +01 , 
3. 20000000 E +01 , 
3. 40000000 E +01. 
3 . 60000000E+01 , 
3 . 80000000E+01 . 
4. 00000000 E+01 . 
4 . 20000000E+01 . 
4 40000000 E+01 , 
4. 60000000 E+01 . 
4 . 80000000E+0 1 . 
5 . OOOOOOOOE+O 1 . 
5 . 37000000E+05 . 
0 

2 . 00000000 E +00 . 
4 . OOOOOOOOE +00 , 
6 . OOOOOOOOE +00 . 
8. OOOOOOOOE +00. 
1 . OOOOOOOOE+O 1 , 
1 . 20000000E+01 , 
1 40000000E+01 • 

1 60000000E+01 . 
1 . 80000000E+01 . 
2. 00000000 E+01 . 
2 . 20000000E+01 , 

2 40000000E +0 1 . 
2 6 0000000 E+01 > 
2 . 80000000E +0 1 , 
3 . OOOOOOOOE+0 1 . 
3 . 20000000E+0 1 . 
3 . 40000000E+0 1 . 
3 . 60000000E+01 > 
3 . 80000000E+01 , 
4. OOOOOOOOE +01 . 
4. 2 0000000 E+01 . 
4 40000000E+01 . 
4 . 60000000E+0 1 > 
4 . 80000000E +0 1 . 
5. OOOOOOOOE+O 1 , 

$ 

l$TAB 

TABLE=6HCDHT .2 
3 ■ OOOOOOOOE+OS > 
0. 

2 ■ OOOOOOOOE+OO , 
4 OOOOOOOOE+OO. 
6. OOOOOOOOE+OO, 
8 OOOOOOOOE+OO. 
1 . OOOOOOOOE+O 1 . 
1 .20000000E+01 , 
1 40000000E+01 , 

1 60000000E+01 , 
1 . 80000000E+01 , 
2. 00000000 E+01 , 
2.20000000E+01 , 

2 40000000E+01 , 
2 . 60000000 E+01 , 
2. 80000000 E+01 , 
3.00000000E+01 , 
3 . 20000000E+0 1 , 

3 40000000 E +0 1 , 
3. 60000000 E+01 , 
3 80000000 E+01 , 
4. 00000000 E+01 , 
4. 20000000 E+01 , 
4 . 40000000E+01 , 
4. 60000000 E+01 , 


5 39782977E-01 , 

6 074753 13E "01 , 
6.75698509E-01 , 
7 . 44588487E *01 , 
8. 15700000E-01 , 
8.87540000E-01 , 
9 . 59380000E -01 , 
1 . 02692000E+00, 
1 .09016000E+00, 
1 15340000E+00, 
1 . 20220000E+00, 
1 .25100000E+00, 
1 .28949000E+00, 
1 .31767000E+00, 
1 .34585000E+00, 

9. 20449438 E -03, 
1 27396078E -02 , 
1 70659970E -02 , 
2. 15988761 E -02 , 
2.33543284E-02, 
2 . 97500000E -02 , 
3 . 85700000E -02 , 
4 . 73900000E -02 , 
5.69150000E-02, 
6.71450000E-02 , 
7.73750000E-02, 
8 . 83450000E -02 , 
9.931 50000E -02 , 
1 . 10300000E-01 , 
1 .21300000E-01 , 
1 32300000E-01 . 
1 42510000E-01 . 
1 -52720000E-01 , 
1 . 62160000E-01 , 
1 70830000E -01 , 
1 79500000E-01 , 
1 859 10000E -01 , 
1 .92320000E-01 , 
1 .97300000E-01 , 
2 .00850000E-0 1 , 
2.04400000E-01 , 


2 . 30000000E+01 . 
2 . 50000000E+01 , 
2 . 70000000E+01 . 
2.90000000E+01 , 
3 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 90000000E+01 . 
4. 10000000 E+01, 
4 . 30000000E+01 , 
4 . 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3 OOOOOOOOE+OO, 
5. OOOOOOOOE+OO. 
7 .OOOOOOOOE+OO, 
9 .OOOOOOOOE+OO, 
1 10000000E+01 , 

1 30000000E+01 , 
1 .50000000E+01 , 
1 . 70000000E+01 , 
1 . 90000000E+01 , 

2 10000000E+01 , 
2 . 30000000E+01 , 
2 . 50000000E+01 , 
2 . 70000000E+0 1 , 
2 . 90000000E +0 1 , 

3 10000000E+01 , 
3.30000000E+01 , 
3 . 50000000E+01 , 
3. 70000000 E+01, 
3 . 90000000E+01 , 
4. 10000000E+01 , 
4 . 30000000 E+01 , 
4. 50000000 E+01, 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


5 . 73567038E -01 , 
6 . 4 1499793E -01 , 
7. 10064539E -01 , 
7.79261646E-01, 
8 . 5 1620000E -01 , 
9 . 23460000E -01 , 
9 . 95300000E -01 . 
1 .05854000E+00, 
1 . 12178000E+00, 
1 17780000E+00, 
1 . 22660000E+00 , 
1 . 27540000E+00 , 
1 .30358000E+00, 
1 .33176000E+00, 


1 08389264E -02 , 
1 483 16963E *02 , 
1 94088435E-02, 
2 . 24467 1 4 1 E -02 , 
2 . 63200000E -02 , 
3 . 41600000E -02 , 
4 . 29800000E -02 , 
5. 18000000 E- 02, 
6. 20300000 E- 02, 
7 . 22600000E -02 , 
8 . 28600000E -02 , 
9 . 38300000E -02 , 
1 . 04800000E -01 , 
1 . 1 5800000E -01 , 
1 . 26800000E -01 , 
1 . 37405000E -01 , 
1 47615000E-01 , 
1 .57825000E-01 , 
1 . 66495000E -01 , 
1 . 75 1650G0E -01 , 
1 . 82705000E -01 , 
1 .8911 5000E -01 , 
1 . 95525000E-01 , 
1 .99075000E-01 , 
2 . 02625000E -01 , 


,5H ALPHA ,6HALTIT0,51,2 , 1 , 1 , 1 , 1 , 1 , 1 . 1 , 1 , 


9 73826087E -02 , 
1 022122 1 4E -01 , 
1 074707 48E -01 , 
1 16979753E -01 , 
1 .3051 1 33 1 E -01 , 
1 44054585E -01 , 
1 67867737E -01 , 

1 9 1448482E -01 , 
2 . 227 1 137 1 E -01 , 
2 . 6 1490092E -01 , 

2 99887372E-01 , 
3 . 5490 1625E -01 , 
4 1087 1250E-01 , 
4 . 74809028 E -01 , 
5 . 46673604E -01 , 
6. 19675000E-01 , 

7 15715000E-01 , 

8 1 1755000E -01 , 
9. 15945000E -01 , 
1 02828500E+00 , 
1 . 14062500E+00 , 
1 .26684500E+00, 
1 3 9 306 500 E +00 , 
1 .52469000E+00, 


1 .OOOOOOOOE+OO 
3 .OOOOOOOOE+OO 
5 .OOOOOOOOE+OO 
7 .OOOOOOOOE+OO 
9 .OOOOOOOOE+OO 
1 . 10000000 E+01 

1 30000000E+01 
1 . 50000000 E+01 
1 .70000000E+01 
1 . 90000000E+01 
2 . 10000000 E+01 
2 . 30000000E+01 
2. 50000000 E+01 
2 . 70000000E+01 

2 90000000E+01 
3 . 10000000E+01 
3 . 30000000E+01 
3 . 50000000E+01 
3 . 70000000E+01 
3 . 90000000E+01 
4. 10000000E+01 
4 . 30000000E+01 
4 . 50000000E+01 
4 . 70000000E+01 


9.97333333E-02, 
1 . 04797 122E -01 , 
1 . 10219355E-01 , 
1 . 23743898E -01 , 
1 .37281670E-01 , 
1 55992797E-01 , 

1 79684720E -01 , 
2 . 03 163258E -01 , 

2 42 15 1422E -01 , 
2 . 80733524E -01 , 
3 . 2707 1201E-01 , 
3 . 82984404E-01 , 
4.38574793E-01 , 
5 . 10838124E-01 , 
5 . 82326 1 4 1 E -01 , 
6.67695000E-01 , 
7.63735000E-01 . 
8.59775000E-01 , 
9.721 15000E-01 , 
1 . 08445500E+00 , 
1 . 20373500E+00 ■ 
1 . 32995500E+00 , 
1 .45617500E+00, 
1 .59320500E+00, 


62- 

63- 

64- 

65- 

66 - 

67- 

68 - 

69- 

70- 

71- 

72- 

73- 

74- 

75- 

76- 

77- 

78- 

79- 

80- 
81- 
82- 

83- 

84- 

85- 

86 - 

87- 

88 - 

89- 

90- 

91- 

92- 

93- 

94- 

95- 

96- 

97- 

98- 

99- 
100 - 
101 - 
102 - 

103- 

104- 

104- 

105- 

106- 

107- 

108- 

109- 

110 - 
111 - 
112 - 

113- 

114- 

115- 

116- 

117- 

118- 

119- 

120 - 
121 - 
122 - 

123- 

124- 

125- 

126- 

127- 

128- 

129- 

130- 
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4 . 80000000E +0 1 . 
5. 00000000 E +01, 
5 . 37000000E+05 , 
0. 

2. 00000000 E +00, 
4. 00000000 E +00, 
6. 00000000 E +00, 
8. 00000000 E +00, 
1 . 00000000 E +0 1 , 
1 .20000000E+01 , 
1 .40000000 E +01 , 
1 .60000000 E +01, 
1 . 80000000 E +01 , 
2. 00000000 E +01 , 
2.20000000E+01, 
2. 40000000 E +01 , 
2. 60000000 E +01, 
2. 80000000 E +01, 
3 . 00000000 £ +0 1 , 
3. 20000000 E +01, 
3 . 40000000E +0 1 , 
3 . 60000000E+01 , 
3 . 80000000E+01 , 
4 . 00000000E +0 1 , 
4 . 20000000E+01 , 
4 40000000E+01 , 
4 . 60000000E+01 , 
4 . 80000000E+01 , 
5. 00000000 E +01, 

$ 

L$TA8 

T ABLE =6HCLVT ,2 
5. 00000000 E- 03, 
0. 

2 . OOOOOOOOE+OO , 
'4. 00000000 E +00, 
6 . 00000000 E+ 00 , 
8 . 00000000 E + 00 , 
1 . OOOOOOOOE +0 1 , 
1 .20000000E+01, 
1 .40000000 E +01, 
1 .60000000E+01 , 
1 .80000000E+01, 
2. OOOOOOOOE +01. 
2 . 20000000E +0 1 , 
2 40000000E+01 , 
2 . 60000000E+01 , 

2 80000000E+01 , 

3 OOOOOOOOE +01, 
3 . 20000000E+01 , 
3 . 40000000E+01 , 
3 60000000 E +0 1 , 

3 80000000E+01 , 
4. OOOOOOOOE +01 , 
4 . 20000000E+01 , 

4 40000000 E +0 1 , 
4. 60000000 E +01, 
4 80000000E+01 , 
5. OOOOOOOOE +01, 
1 . 03000000E -02 , 
0. 

2. OOOOOOOOE +00, 
4. OOOOOOOOE +00, 
6. OOOOOOOOE +00, 
8. OOOOOOOOE +00, 
1 .OOOOOOOOE +01, 
1 . 20000000E+01 , 
1 40000000E +0 1 , 

1 60000000E+01 , 

1 . 80000000 E +01, 


1 . 66172000E+00, 
1 . 79875000E+00, 

6.84452809E-01, 
7 . 24948627E -01 , 
7. 73231049E-01 , 
8 . 44880153E -0 1 , 
9 . 66479 159E -01 , 
1 .08697500E+00, 
1 .21427500E+00, 
1 .34157500E+00, 
1 . 46999400E+00, 
1 .59953200E+00, 
1 .72907000E+00, 
1 . 85840200E+00, 
1 .98773400E+00, 
2. 1 1653500E+00, 
2 . 24480500E+00 , 
2 . 37307500E+00, 
2 . 50044500E+00, 
2 62781500E+00, 
2 75382000E+00 , 
2 87846000E+00, 
3. 0031 0000 E +00, 
3. 12228000E+00, 
3 . 24146000E+00, 
3 . 35755500E+00, 
3 . 47056500E+00, 
3 . 58357500E+00, 


, 5HALPHA .6HVINV 

1 .46464286E-02, 
4 . 25473300E-02 , 
7 . 02738964E -02 , 
1 .04420102E-01 , 
1 45075276E-01, 
1 .85721212E-01 , 
2 . 38744784 E -01 , 
2.91953236E-01 , 
3 . 5004 1208 E -0 1 , 
4. 13206959E -01 , 
4 . 76757 143E -0 1 , 
5 . 44962558E -0 1 , 

6 13558 139E-01 , 
6.82654131E-01 , 

7 52433468E-01 , 
8 . 27800000E -01 , 
9 .00010000E -01 , 
9.72220000E-01, 
1 . 03987000E+00, 
1 . 10296000E+00 , 
1 16605000E+00 , 
1 . 21374000E+00, 
1 . 26 1 43000E+00 , 
1 .29905500E+00, 
1 .32661500E+00, 
1 . 354 1 7500E+00 , 

1 4892857 1 E -02 , 
4.2591 6877E -02 , 
7.01326577E-02, 
1.041 38937E -01 , 
1 . 44689774E -01 , 
1 . 85230435E -01 , 
2 . 37959444E -01 , 
2 . 908828 13E-01 , 
3 . 48780309E-01 , 
4. 1 183362 1 E-01 , 


4 . 90000000E+0 1 , 


1 73023500E+00, 


1 . OOOOOOOOE+OO , 
3 . OOOOOOOOE+OO , 
5. OOOOOOOOE+OO, 
7 . OOOOOOOOE+OO , 
9. OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 . 30000000E+01 , 
1 . 50000000E+01 , 

1 70000000E+01 , 
1 . 90000000E+01 , 
2. 10000000E+01 , 
2 . 30000000E+01 , 
2 . 50000000E+01 , 

2 70000000E+01 , 

2 90000000E+01 , 
3 . 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 

3 70000000E+01 ■ 

3 90000000E+01 , 

4 10000000E+01 , 
4 . 30000000E+01 , 
4 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


,51,8,1,1,1,1,1, 

1 .OOOOOOOOE+OO, 
3. OOOOOOOOE+OO, 
5. OOOOOOOOE+OO, 
7 .OOOOOOOOE+OO, 
9. OOOOOOOOE+OO, 
1 10000000E+01 , 
1 .30000000E+01, 
1 50000000E+01 , 

1 70000000E+01 , 
1 .90000000E+01 , 

2 10000000E+01 , 
2 30000000E+01 , 
2 . 50000000E+01 , 
2 70000000E+01 , 

2 90000000E+01 , 

3 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 

3 90000000E+01 , 

4 10000000E+01 , 
4 . 30000000E+01 , 
4 . 50000000E+01 . 
4 . 70000000E+01 . 
4 . 90000000E+01 . 


1 .OOOOOOOOE+OO, 
3. OOOOOOOOE+OO, 
5. OOOOOOOOE +00, 
7. OOOOOOOOE +00, 
9. OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 .30000000E+01 , 
1 .50000000E+01 , 
1 . 70000000E+01 , 
1 .90000000E+01 , 


7 . 03390696E -01 , 
7 . 483901 13E -01 , 
7 . 99 140136E -01 , 
9 . 06022390E -01 . 
1 . 02663500E+00 , 
1 . 15062500E+00, 
1 .27792500E+00, 
1 .40522500E+00, 
1 . 53476300E+00, 
1 .66430100E+00, 
1 .79373600E+00, 
1 . 92306300E+00 , 
2 . 05240000E+00 , 
2. 18067000E+00, 
2 . 30894000E+00, 
2 . 43676000E+00 , 
2 . 564 13000E+00 , 
2 69150000E+00, 
2 . 8 16 14000E+00 , 
2.94078000E+00, 
3 . 06269000E+00 , 
3. 18187000E+00, 
3 . 30105000E+00 , 
3.41 406000E +00 , 
3 . 52707000E+00 , 


,1,1, 

2 . 86227845E -02 , 
5. 642875 15E -02, 
8 . 40882353E -02 , 
1 .24748993E-01 , 
1 .65399262E-01 , 
2. 12207849E-01 , 
2.65327825E-01 , 
3. 18617669E-01 , 
3.81 573 1 17E -01 , 
4 . 44936796E -01 , 
5. 10854598E-01 , 
5 . 79 199794E -01 , 
6 . 48030103E -01 , 
7. 1 7458692E -01 , 
7 . 89700000E -01 , 
8.63905000E-01 , 
9.361 15000E -01 , 
1 . 00832500E+00 , 
1 .07141500E+00, 
1 . 13450500E+00 , 
1 . 18989500E+00, 
1 .23758500E+00, 
1 .28527500E+00, 
1 .31283500E+00, 
1 . 34039500E+00, 


2 . 87657296E -02 , 
5 . 63 785969 E -02 , 
8 . 38588235E -02 , 
1 .24415788E-01 , 
1 .64961229E-01 , 
2. 11568432E-01 , 
2 . 643989 19E -01 , 
3. 17407519E-01 , 
3 . 80257749E -01 , 
4.43502191E-01 , 


-- 131 - - 
-- 132- - 
-- 133* 
-- 134- - 
-- 135- - 
-- 136-- 
-- 137- - 
-- 138-- 
-- 139- - 
-- 140* - 
-- 141 - - 
-- 142 - - 
-- 143 - - 
-- 144- - 
-- 145- - 
-- 146 - - 
-- 147- - 
-- 148 - - 
-- 149- - 
-- 150* - 
-- 15 1 - * 
-- 152-- 
-- 153- - 
-- 154- - 
-- 155- - 
-- 156-- 
-- 157- - 
-- 158- - 
-- 159- - 
-- 160- - 
-- 160- - 
-- 16 1 - - 
-- 162-- 
-- 163-- 
-- 164- * 
-- 165- - 
-- 166- - 
-- 167 - - 
-- 168 - - 
-- 169- - 
-- 1 70- - 
-- 17 1 - - 
-- 1 72 - - 
-- 1 73 - - 
-- 1 74- - 
-- 175 — 
- - 1 76- - 
-- 177- - 
-- 178-- 
-- 1 79 - - 
-- 180* - 
-- 18 1 - - 
-- 182* * 
-- 183- * 
-- 184- - 
-- 185- - 
-- 186-- 
-- 187 - - 
-- 188- - 
-- 189- - 
-- 190- - 
-- 19 1 - - 
-- 192- - 
-- 193- - 
-- 194- - 
-- 195- - 
-- 196- - 
-- 197- - 
-- 198 - - 
-- 199- - 
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ORIGINAL PAGE 13 

OF POOR QUALITY 


2. OOOOOOOO E +01 , 
2. 20000000 E +01 , 
2. 40000000 E +01 , 
2. 60000000 E +01 , 
2. 80000000 E +01 , 
3. OOOOOOOO E +01 , 
3. 20000000 E +01 . 
3. 40000000 E +01 , 
3 . 60000000E+01 , 
3 . 80000000E+01 . 
4. OOOOOOOO E +01 ■ 
4. 20000000 E +01 . 
4 . 40000000E+01 . 
4 . 60000000E+01 , 
4 . 80000000E+01 . 
5 . OOOOOOOOE +0 1 . 
1 . 48000000E -02 , 
0. 

2. OOOOOOOOE +00, 
4 . OOOOOOOOE +00 , 
6 . OOOOOOOOE +00 , 
8. OOOOOOOOE +00, 
1 .OOOOOOOOE+OI , 
1 .20000000E+01 , 
1 40000000E+01 , 
1 60000000E+01 , 
1 .80000000 E +01 , 
2. OOOOOOOOE +01 , 
2. 20000000 E +01 , 
2 . 40000000E+01 , 
2 . 60000000E +0 1 , 
2 . 80000000E+01 , 
3. OOOOOOOOE+OI , 
3 . 20000000E+01 , 
3 . 40000000 E +01 . 
3 . 60000000E+01 , 
3.80000000E+01 , 
4. OOOOOOOOE +01 , 
4 . 20000000E+01 . 
4. 40000000 E +01 , 
4.60000000E+01 , 
4 . 80000000E+01 , 
5. OOOOOOOOE +01 , 
2. 08000000 E- 02, 
0. 

2 . OOOOOOOOE+OO , 
4 OOOOOOOOE+OO, 
6. OOOOOOOOE +00, 
8. OOOOOOOOE +00, 
1 .OOOOOOOOE+OI , 
1 .20000000E+01 , 
1 .40000000 E +01 , 
1 .60000000E+01 , 

1 80000000E+01 , 

2 OOOOOOOOE+OI, 
2 . 20000000E+01 , 
2. 40000000 E +01 , 
2.60000000E+01 , 
2. 80000000 E +01 , 
3. OOOOOOOOE+OI , 
3.20060000E+01 , 
3 . 40000000 E +01 , 
3. 60000000 E +01 , 

3 80000000 E +01 , 
4. OOOOOOOOE+OI , 
4 . 20000000E+01 , 
4 . 40000000E+01 , 
4 . 60000000E+01 , 
4.80000000E+01 , 
5. OOOOOOOOE+OI , 


4.75253163E-01 , 
5.43145521E-01 , 
6 . 1 14408 17E-01 , 
6 . 80353476E -01 , 
7 . 50032232E-01 , 
8.25275000E-01 , 
8 . 97415000E-01 , 
9 . 69555000E-01 , 
1 . 037 14000E+00 , 
1 . 10017000E+00, 
1 . 16320000E+00 , 
1 .21 132000E+00, 
1 . 25944000E+00 , 
1 . 29737000E+00 , 
1 .3251 1000E+00, 
1 . 35285000E+00, 

1 .48928571E-02, 
4 . 25225253 E *02 , 
7. 00035747 E -02, 
1 .03926982E-01 , 
1 44368098 E *01 , 
1 .84800000E-01 , 
2 . 37542553 E *01 , 
2 . 90472643E *01 , 

3 48269013E-01 , 

4 1 1 1238 1 9E -01 , 
4 . 74357 143E -01 , 
5. 422178265-01 , 
6 . 104793 18E -01 , 
6 79290655E-01 , 
7 . 488 12 100E -01 , 
8 . 23900000E -01 , 
8 . 95990000E-01 , 
9.68080000E-01 , 
1 . 03566000E+00 , 
1 . 09873000E+00 , 
1 16 180000E+00 , 
1 .21008000 E +00, 
1 .25836000E+00, 
1 .29641000E+00, 
1 32423000E+00 , 
1 35205000E+00, 

1 48857 143E -02 , 
4. 26460606 E -02, 
7 . 02564706 E -02 , 
1 . 041 86043E -01 , 
1 . 445 15 184E -0 1 , 
1 . 84839738 E -0 1 , 
2 . 375 10927E -01 , 
2 . 90380495E -01 , 

3 480968 1 9E -01 , 

4 10840679E-01 , 
4 . 73957 1 43E -01 , 
5.41 777826E -01 , 
6.09999318E-01 , 
6 . 78723589 E -01 , 
7 . 48 103 126E -01 , 
8.23050000E-01, 
8.951 10000E-01 , 
9.671 70000E -01 , 
1 . 03471000E+00, 
1 . 09773000E+00 , 
1 16075000E+00, 
1 . 20923000E+00, 
1 .25771000E+00, 
1 .29589000E+00, 
1 .32377.000E+00, 
1 .35165000E+00, 


2 1 0000000 E +01 ■ 
2 . 30000000E+01 , 
2 . 50000000E+01 , 

2 70000000E+01 , 
2 . 90000000E+01 , 

3 1 0000000 E +0 1 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 . 90000000E+01 , 
4. 10000000E+01 , 
4. 3 0000000 E +01, 
4 . 50000000E+01 , 

4 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3. OOOOOOOOE +00, 
5. OOOOOOOOE +00, 
7. OOOOOOOOE +00, 
9 OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 .30000000E+01 , 
1 . 50000000E+01 , 

1 70000000 E +01 , 
1 .90000000E+01 , 

2 1 0000000 E +0 1 , 
2 . 30000000E+01 , 
2 . 50000000E+01 , 
2 . 70000000E+01 , 
2 . 90000000E +0 1 , 
3 . 1 0000000 E +01 , 
3 . 30000000E+0 1 , 
3 . 50000000E+01 , 
3. 70000000 E +01, 
3. 90000000 E +01 , 
4. 1 0000000 E +01 , 
4 . 30000000E +0 1 , 
4 . 50000000E +0 1 , 
4 . 70000000E+01 , 
4 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3. OOOOOOOOE +00, 
5. OOOOOOOOE +00, 
7 OOOOOOOOE+OO, 
9 OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 .30000000E+01 , 
1 50000000 E +01 , 
1 70000000E+01 , 

1 90000000E+01 , 

2 10000000E+01 , 

2 30000000 E +0 1 , 
2 . 50000000E+01 , 
2 . 7 OOOOOOOE +0 1 , 
2 . 90000000E+01 , 

3 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3. 70000000 E +01 , 
3 . 90000000E+01 , 

4 1 OOOOOOOE +0 1 , 
4 30000000E+01 , 
4 . 50000000E+01 , 
4 . 70000000E +0 1 , 
4 . 90000000E+01 , 


5 . 09193591E -01 , 
5 77230735E-01 , 
6 . 45768041E -01 , 
7. 15111494E-01 , 
7 . 87275000E -01 , 
8 . 6 1345000E -01 , 
9 . 33485000E -01 , 
1 .00562500E+00, 
1 .06865500E+00, 
1 13168500E+00, 
1 . 18726000E+00 , 
1 .23538000E+00, 
1 .28350000E+00, 
1 3 1 124000E+00 , 
1 33898000E+00 , 


2. 8729705 IE -02, 
5 . 62785680E -02 , 
8 . 3 702 1 505 E - 02 , 
1 .24148841E-01 , 

1 64585063E-01 , 

2 1 1 145700E -01 , 

2 63986169E-01 , 

3 16998496E-01 , 

3 79646234E-01 , 

4 42695922E-01 , 
5.08280279E-01 , 
5 . 76287000E -01 , 

6 44787010E-01 , 

7 1 3969013E -01 , 
7.85940000E-01 , 

8 59945000E-01 , 

9 32035000E-01 , 
1 00412500E+00, 
1 .067 19500E+00, 
1 . 13026500E+00 , 
1 18594000E+00 , 
1 23422000E+00 , 
1 .28250000E+00, 
1 3 1032000E+00 , 
1 .33814000E+00, 


2 . 87880965E -02 , 
5 . 646692 12E -02 , 
8 . 401 93548 E -02 , 
1 .24351260E-01 , 
1 .64677967E-01 , 
2.11 148206E-01 , 

2 63923105E-01 , 

3 16879323E-01 , 
3 . 79419351 E -01 , 

4 42355049E-01 , 
5.07860279E-01 , 
5 . 75827000E -01 , 
6 . 44287010E -01 , 
7 13332239E -01 , 
7 85160000E-01 , 
8 . 59080000E -01 , 
9.31 140000E-01 , 
1 .00320000E+00, 
1 .06622000E+00, 
1 12924000E+00 , 
1 . 18499000E+00 , 
1 .23347000E+00, 
1 28 195000E+00 , 
1 .30983000E+00, 
1 3377 1000E+00 • 


200 - 

201 - 

202 - 

203- 

204- 

205- 

206- 

207- 

208- 

209- 

210 - 
211 - 
212 - 

213- 

214- 

215- 

216- 

217- 

218- 

219- 

220 - 
221 - 
222 - 

223- 

224- 

225- 

226- 

227- 

228- 

229- 

230- 

231- 

232- 

233- 

234- 

235- 

236- 

237- 

238- 

239- 

240- 

241- 

242- 

243- 

244- 

245- 

246- 

247- 

248- 

249- 

250- 

251- 

252- 

253- 

254- 

255- 

256- 

257- 

258- 

259- 

260- 
261- 
262- 

263- 

264- 

265- 

266- 

267- 

268- 
269- 
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2 . 45000000E -02 , 

0. , 1 49043165E-02, 
2 . OOOOOOOOE+OO » 4 .25933840E -02 , 
4. 00000000 E +00, 7.01096489E-02, 
6. OOOOOOOOE+OO, 1 ,04000741 E -01 , 
8 .OOOOOOOOE+OO, 1 .44358409E-01 , 
1 .00000000 E +01 , 1 .84710526E-01 , 
1 .20000000 E +01 , 2 . 3741S705E -01 , 
1 40000000 E +01 , 2 . 903 17019E-01 , 
1 .60000000 E +01 , 3 . 47993097 E -01 , 
1 -80000000E+01 , 4. 10631327E-01 , 
2. 00000000 E +0 1 , 4 73645392E -01 , 
2 . 20000000E+01 , 5 41402120E-01 , 
2. 40000000 E +0 1 , 6.09531081E-01 , 
2. 60000000 E +0 1 , 6 . 78 192470E -01 , 
2 . 80000000E+0 1 , 7 . 47538288E-01 , 
3. 00000000 E +0 1 , 8 .22275000E -01 , 
3. 20000000 E +0 1 , 8 .94325000E-01 , 
3. 40000000 E +0 1 , 9 .66375000E-01 , 
3 . 60000000E+01 , 1 .03393000E+00, 
3 . 80000000E+01 , 1 . 09699000E+00 , 
4 . OOOOOOOOE+O 1 , 1 . 16005000E+00, 
4 . 20000000E +0 1 , 1 .20853000E+00, 
4.40000000E+01 , 1 .25701000E+00, 
4 . 60000000E+0 1 , 1 . 29523000E+00 , 
4 . 80000000E+0 1 , 1 . 323 19000E+00 , 
5. OOOOOOOOE+O 1 , 1 351 15000E+00, 
3 . 73000000E -02 , 

0. , 1 48971223E -02 , 
2. OOOOOOOOE+OO, 4 . 283941 70E -02 , 
4. OOOOOOOOE+OO, 7 . 05992072E -02 , 
6. OOOOOOOOE+OO, 1 0457 1655E -01 , 
8. OOOOOOOOE+OO, 1 . 448536 1 7E -01 . 
1. OOOOOOOOE+O 1, 1 85129694E-01 , 
1 .20000000E+01 , 2.375931 14E-01 , 
1 .40000000 E +01 , 2 . 90249170E -01 , 
1 .60000000E+01 , 3 . 477 18108E -01 , 
1 80000000E+01 , 4 10179134E -01 , 
2. OOOOOOOOE+O 1 , 4 . 73006826 E -01 , 
2 . 20000000E+01 , 5 40789620E -01 , 
2. 40000000 E +0 1 , 6 . 088 16258E -01 , 
2 . 60000000E+01 , 6 . 77340073E-01 , 
2 . 80000000E+0 1 , 7 465 1 2333E -01 , 
3 -OOOOOOOOE+O 1 , 8 . 20400000E -01 , 
3 . 20000000E+0 1 , 8 . 92350000E -01 , 
3. 40000000 E +01, 9 . 64300000E -01 , 
3 . 60000000E+01 , 1 . 03 183500E+00 , 
3 . 80000000E+01 , 1 09495500E+00 , 

4 OOOOOOOOE+O 1, 1 . 15807500E+00, 
4. 20000000 E +0 1 , 1 . 20675500E+00, 
4 . 40000000E+01 , 1 .25543500E+00, 
4 . 60000000E +0 1 , 1 . 29379000E+00, 
4 . 80000000E+01 , 1 . 32 182000E+00,' 
5. 00000000 E +01, 1 .34985000E+00, 

5 70000000E-02, 

0. , 1 -52669065E-02, 
2. OOOOOOOOE+OO, 4 . 27977440E-02 , 
4. OOOOOOOOE+OO, 7 .01564892E -02, 
6. OOOOOOOOE+OO, 1 .04128818E-01 , 
8. OOOOOOOOE+OO, 1 . 447957 13E -01 , 
1 .00000000 E +01, 1 .85445852E-01 , 
1. 20000000 E +01, 2 . 378 10337E -01 , 
1 .40000000 E +01 , 2 . 90362002E-01 , 
1 .60000000 E +01 , 3 47704080E-01 , 
1 80000000 E +01 , 4. 10021703E-01, 
2. 00000000 E +01, 4.72708874E-01, 
2. 20000000 E +01, 5 . 40566855E -01 , 
2 . 40000000E+01 , 6.08385743E-01, 
2 . 60000000E+01 , 6 .76708562E-01 , 
2. 80000000 E +01, 7 . 45670728E -01 , 
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3. 00000000 E +01 , 
3 . 20000000E+01 , 
3. 40000000 E +01. 
3 . 60000000E+01 , 
3 . 80000000E+01 , 
4. 00000000 E +01 , 
4. 20000000 E +01 , 
4 . 40000000E+01 , 
4 . 60000000E+01 . 
4. 80000000 E +01 . 
5. 00000000 E +01. 
7 . 44000000E -02 , 
0. 

2 . OOOOOOOOE +00 , 
4.00000000E+00, 
6. OOOOOOOOE +00. 
8 .OOOOOOOOE +00. 
1 .OOOOOOOOE +01 , 
1 .20000000E+01 . 
1 .40000000E+01 , 
1 .60000000E+01 . 
1 .80000000E+01 , 
2 . 00000000 E + 0 1 , 
2 . 20000000E+01 , 
2 . 40000000E+01 , 
2 . 60000000E+01 . 
2 . 80000000E+01 , 
3.00000000E+01 . 
3 . 20000000E+01 > 

3 40000000E+01 . 
3 . 60000000E+01 . 
3 . 80000000E+01 , 
4 . OOOOOOOOE +01 , 
4 . 20000000E+01 . 
4 . 40000000E+01 . 

4 60000000E+01 . 
4 . 80000000E+01 , 
5. OOOOOOOOE +01 , 

$ 

l$TAB 

T ABLE =6HCDVT .2 
5. OOOOOOOOE -03. 
0 

2 . OOOOOOOOE+OO . 
4. OOOOOOOOE +00. 
6. OOOOOOOOE +00, 
8. OOOOOOOOE +00, 
1 .OOOOOOOOE +01 , 
1 .20000000E+01 , 
1 40000000E+0 1 , 

1 60000000 E +0 1 , 
1 . 80000000E+01 . 
2. OOOOOOOOE +01 , 
2 . 20000000E+01 , 

2 40000000 E +0 1 , 
2. 60000000 E +01 , 
2 80000000 E +01 , 
3.00000000E+01 . 
3 . 20000000E +0 1 , 
3 . 40000000 E +01 , 
3.60000000E+01. 
3 . 80000000E+01 , 
4. OOOOOOOOE +01 , 
4 . 20000000E+01 , 
4 40000000E+01 , 
4 . 60000000E+01 , 
4 . 80000000E+01 , 
5. OOOOOOOOE +01 , 
1 03000000E -02 . 
0. 


8 18075000E -01 , 
8 . 89955000E -01 , 
9 . 6 1835000E -01 , 
1 .02932500E+00 , 
1 .09242500E+00, 
1 . 15552500E+Q0 . 
1 . 20438500E+00 , 
1 . 25324500E+00 , 
1 .29173000E+00, 
1 .31984000E+00, 
1 .34795000E+00, 

1 54 130435E -02 . 
4.2616 1832E-02 , 
6 . 96270748E -02 , 
1 .030759 10E -01 , 
1 . 43063966E -01 , 
1 . 83045852E -01 , 
2 . 36086535E -01 , 
2 . 892875 18E -01 . 
3 . 46924729E -01 , 
4.09199247E-01 , 
4 7 18406.1 4E -01 . 
5 . 39782977E -01 , 
6 . 074753 1 3E -0 1 , 

6 75698509E-01 , 

7 44588487E-01 , 
8. 1S700000E-01 . 
8 . 87540000E -01 . 
9 . 59380000E -01 , 
1 . 02692000E+00 , 
1 09016000E+00, 
1 15340000E+00, 
1 20220000E+00, 
1 .2S100000E+00, 
1 .28949000E+00, 
1 .31767000E+00, 
1 34585000E+00, 


, 5HALPHA , 6HVINV 

4 61507143E-02, 

5 33939295E -02 , 
6 . 09642 1 1 7E -02 , 
7 . 237 1 7058E -02 ■ 
8 . 74465792 E -02 , 
1 02536364E -01 , 
1 .28157208E-01 , 
1 53509668E -01 , 
1 . 85489352E -01 , 
2 . 239 1 5035E -01 , 
2.61911224E-01 . 

3 1 7964993E -0 1 , 
3 . 7342949 1 E -01 , 
4 . 370953 19E -01 , 
5 . 089 15933E -01 , 
5 . 85725000E -01 , 

6 82605000E-01 , 
7 . 79485000E -01 , 
8.84770000E-01, 
9 . 98460000E -01 , 
1 . 11215000E+00, 
1 .24043000E+00, 
1 .36871000E+00, 
1 . 502 10000E+00 , 
1 . 6 4060000 E +00, 
1 .77910000E+00, 

4 6 5000000 E -02, 


3. 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 . 90000000E+01 , 
4. 1 0000000 E +01 , 
4 . 30000000E+01 , 
4 . 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3 OOOOOOOOE+OO, 
5 OOOOOOOOE+OO, 
7 OOOOOOOOE+OO, 
9 .OOOOOOOOE+OO, 
1 . 1 0000000 E +01 , 
1 .30000000 E +01 , 
1 . 50000000E+01 , 
1 . 70000000E+01 , 

1 90000000E+01 , 

2 10000000E+01 , 

2 30000000E+01 , 
2. 50000000 E +01 , 
2. 70000000E+01 , 
2 . 90000000E+01 , 
3. 10000000E+01 , 
3 . 30000000E+01 , 

3 50000000E+01 , 
3 70000000E+01 , 

3 90000000E+01 , 

4 10000000E+01 , 
4. 30000000 E +01 , 
4 . 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


,51,8,1,1,1,1.1, 

1 .OOOOOOOOE+OO, 
3 .OOOOOOOOE+OO, 
5 OOOOOOOOE+OO, 
7 .OOOOOOOOE+OO, 
9. OOOOOOOOE +00, 
1 10000000E+01 , 
1 .30000000E+01 . 
1 . 50000000E+01 , 

1 70000000E+01 , 
1 . 90000000E+01 , 
2. 10000000E+01 , 
2 . 30000000E+01 , 

2 50000000 E +01 , 
2 . 70000000E+01 , 
2 90000000E+01 , 
3 . 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 . 90000000E+01 , 
4. 10000000E+01 , 
4. 3 0000000 E +01 , 
4 . 50000000E+01 , 
4 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 


8 . 54015000E -01 , 
9 . 25895000E-01 , 
9 . 97775000E -01 , 
1 . 06087500E+00, 
1 . 12397500E+00, 
1 . 17995500E+00, 
1 . 2288 1500E+00, 
1 27767500E+00, 
1 30578500E+00, 
1 . 33389500E+00, 


2 . 90433333E-02, 
5 . 61415108E -02, 
8 . 30790323 E -02 , 
1 . 23070807E -01 , 
1 . 63055590E-01 , 
2 . 09544373E -01 , 
2.62668672E-01, 
3. 15940152E-01, 
3 . 78013237E-01 ■ 
4 . 40476852E -01 , 
5 .05921 147E -01 ■ 
5 . 73567038E -01 , 

6 41499793E -01 , 

7 10064539E -01 , 
7 7926 1646E -01 , 
8 . 5 1620000E -01 , 
9 . 23460000E -01 , 
9 . 95300000E -01 , 
1 . 05854000E+00 , 
1 . 121 78000E+00 , 
1 17780000E+00, 
1 .22660000E+00, 
1 27540000E+00 , 
1 30358000E+00 , 
1 . 33 1 76000E+00 , 


, 1 , 1 . 

4 . 972372 1 6E -02 , 
5 71450416E-02, 
6 . 484 1 1 765E -02 , 
7. 9907038 IE -02, 
9 . 49898290E -02 , 
1 . 1 5383294E -01 , 
1 .40864190E-01, 

1 66098496E-01 . 

2 04759 134E -01 , 
2.42963689E-01, 
2 . 89909726E -01 , 
3.45800284E-01, 
4 . 0086536 1 E -01 , 
4 73 108038E -01 , 
5 . 46885000E-01 , 
6 . 34 165000E -01 , 
7 . 3 1045000E -01 . 
8 . 27925000E -01 , 
9 . 416 15000E -01 , 
1 . 05530500E+00, 
1 . 17629000E+00, 
1 . 30457000E+00 , 
1 43285000E+00, 
1 .57135000E+00i 
1 70985000E+00, 


5 .02267470E -02 , 


340- 

341- 

342- 

343- 

344- 

345- 

346- 

347- 

348- 

349- 

350- 

351- 

352- 

353- 

354- 

355- 

356- 

357- 

358- 

359- 

360- 

361- 

362- 

363- 

364- 

365- 

366- 

367- 

368- 

369- 

370- 

371- 

372- 

373- 

374- 

375- 

376- 

377- 

378- 

378- 

379- 

380- 

381- 

382- 

383- 

384- 

385- 

386- 

387- 

388- 

389- 

390- 

391- 

392- 

393- 

394- 

395- 

396- 

397- 

398- 

399- 

400- 

401- 

402- 

403- 

404- 

405- 

406- 

407- 

408- 
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2 .OOOOOOOOE+OO, 
4 . OOOOOOOOE+OO , 
6 . OOOOOOOOE+OO , 
8 OOOOOOOOE+OO, 
1 . 00000000 E +0 1 , 
1 . 20000000E +0 1 , 
1 40000000 E +0 1 , 
1 .60000000E+01 , 
1 . 80000000E+01 , 
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3 . 30000000E+01 , 
4 . 00000000E +0 1 , 
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3. OOOOOOOOE +01 , 
3 . 20000000E+01 , 
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2.39718792E-01 , 
2.93204541E-01 , 
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683- 

1 65485293E-01 , 

-- 
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-- 
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1. 60000000 E +01, 
1. 80000000 E +01, 
2. 00000000 E +01, 
2. 20000000 E +01, 
2. 40000000 E +01, 
2. 60000000 E +01, 
2 . 80000000E+01 , 
3. 00000000 E +01 , 
3. 20000000 E +01, 
3 . 40000000E+01 , 
3 . 60000000E+01 , 
3. 80000000 E +01. 
4. 00000000 E +01, 
4 . 20000000E+01 , 
4 . 40000000E+01 , 
4 . 60000000E+01 , 
4 . 80000000E+01 , 
S. 00000000 E +01, 
1 . 00000000 E +01 , 
0. 

2 . OOOOOOOOE+OO , 
4. 00000000 E +00, 
6. OOOOOOOOE+OO, 
8. OOOOOOOOE+OO, 
1 .OOOOOOOOE+OI, 
1 . 20000000E+0 1 , 
1 . 40000000E+01 , 
1 .60000000E+01, 
1 . 80000000E+01 , 
2 . OOOOOOOOE +0 1 , 
2 . 20000000E+01 , 
2 . 40000000E+01 , 
2 . 60000000E +01, 
2 . 80000000E +0 1 , 
3 . OOOOOOOOE+O 1 , 
3 . 20000000E+01 , 
3.40dbOOOOE+C1, 
3 . 60000000E+01 , 
3 . 80000000E +0 1 , 
4. OOOOOOOOE +01 , 
4. 20000000 E +0 1 , 
4 40000000 E +01 , 
4 . 60000000E+01 , 
4 . 80000000E +0 1 , 
5. OOOOOOOOE +01, 
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2 . OOOOOOOOE+OO , 
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2. OOOOOOOOE+OO, 
4. OOOOOOOOE +00, 
6. OOOOOOOOE +00, 
8. OOOOOOOOE+OO, 
1 .OOOOOOOOE+OI, 
1 . 20000000E+01 , 
1 .40000000 E +01, 
1 .60000000E+01 , 
1 . 80000000E+01 , 
2. OOOOOOOOE+OI , 
2 . 20000000E+01 , 
2. 40000000 E +01, 
2. 60000000 E +01, 
2 . 80000000E+01 , 
3. OOOOOOOOE+OI , 
3 . 20000000E+01 , 
3 . 40000000E+01 , 
3 . 60000000E+01 , 
3 . 80000000E+01 , 
4. OOOOOOOOE +01 , 


3. 500966 10E -01 , 
4. 13196102E-01 , 
4 . 769976 1 1 E -01 , 
5.45741749E-01 . 
6 . 142802 19E-01 , 
6 . 83 160000E -01 , 
7 . 529838 10E -01 , 
8 . 28400000E -01 . 
9 . 00690000E -01 , 
9 . 72980000E -01 , 
1 .04068500E+00, 
1 . 10380500E+00, 
1 . 16692500E+00, 
1 .21444500E+00, 
1 .26196500E+00, 
1 .29943500E+00, 
1 . 32685500E+00 , 
1 .35427500E+00, 

1 . 46464286E -02 , 
4 . 25473300E-02 , 
7 . 02738964E -02 , 
1 .04440102E-01 , 
1 45 135276E -01 , 
1 . 858212 12E -01 , 
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2.20588705E-01, 
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6.74845000E-01, 
7 71915000E-01 , 
8 . 77890000E -01 , 
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2 . 30000000E+01 , 
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-- 702 -- 
-- 703- - 
-- 704- - 
-- 705 - - 
-- 706- - 
-- 707 - - 
-- 708 -- 
-- 709- - 
-- 7 10 - - 
-- 7 1 1 - - 
-- 7 1 2 - - 
-- 7 13- - 
-- 7 14- - 
-- 7 1 5 - - 
-- 7 16- - 
-- 7 1 7 - - 
-- 7 18 - - 
-- 719- - 
-- 720- - 
-- 72 1 - - 
-- 722-- 
-- 723-- 
-- 724-- 
-- 725-- 
-- 726-- 
-- 727-- 
-- 728-- 
-- 729-- 
-- 730- - 
-- 73 1 — 
- - 732-- 
-- 733-- 
-- 733-- 
-- 734-- 
-- 735-- 
-- 736-- 
-- 737-- 
-- 738-- 
-- 739-- 
-- 740- - 
-- 741 -- 
-- 742-- 
-- 743-- 
-- 744-- 
-- 745-- 
-- 746-- 
-- 747-- 
-- 748-- 
-- 749-- 
-- 750- - 
-- 75 1 - - 
-- 752-- 
-- 753-- 
-- 754-- 
-- 755-- 
-- 756-- 
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ORIGINAL PAGE IS 
Of POOR QUALITY 


4.20000000E+01 , 
4. 40000000 E +01 , 
4 . 60000000E+01 , 
4. 80000000 E +01 , 
5.00000000E+01 , 
4. 00000000 E +00, 
0. 

2. 00000000 E +00, 
4. 00000000 E +00, 
6. 00000000 E +00, 
8 . 00000000 E+ 00 , 
1. 00000000 E +01, 
1 .20000000 E +01 , 

1 40000000E+01 , 
1 .60000000E+01 . 
1 .80000000E+01 , 

2 OOOOOOOOE+OI , 
2 . 20000000E+01 , 
2 • 40000000E+01 , 
2 . 60000000E+0 1 , 

2 - 80000000E+01 , 
3 . OOOOOOOOE+O 1 , 
3 . 20000000E +0 1 , 

3 40000000 E +01 , 

3 S0000000E+01 , 
3 . 80000000E+01 , 

4 OOOOOOOOE+OI , 
4 . 20000000E +0 1 , 
4 40000000E+01 , 
4 . 60000000E +0 1 , 
4 . 80000000E+01 , 
5. OOOOOOOOE+OI , 
6 ■ OOOOOOOOE+OO , 
0- 

2- OOOOOOOOE+OO, 
4 OOOOOOOOE+OO, 
6 OOOOOOOOE+OO, 
8. OOOOOOOOE+OO, 
1 .OOOOOOOOE+OI , 
1 .20000000E+01 , 
1 .40000000E+01 , 
1 60000000E+01 , 
1 .80000000E+01 , 
2. OOOOOOOOE+OI , 
2 . 20000000E+0 1 , 
2 . 40000000E+01 , 
2 . 60000000E+01 , 
2 ■ 80000000E+01 , 
3. OOOOOOOOE+OI , 
3 . 20000000E+0 1 , 
3 . 40000000E+01 . 
3 ■ 60000000E +0 1 , 

3 80000000E+01 , 
4. OOOOOOOOE+OI , 
4.20000000E+01 , 
4 . 40000000E+01 , 
4. 60000000 E +01 , 

4 80000000E+01 , 
5. OOOOOOOOE+O A , 
8. OOOOOOOOE+OO, 
0. 

2. OOOOOOOOE+OO, 
4. OOOOOOOOE+OO, 
6. OOOOOOOOE+OO, 
8 OOOOOOOOE+OO, 
1 .OOOOOOOOE+OI , 
1 . 20000000E+01 , 
1 .40000000 E +01 , 
1 .60000000 E +01 , 
1 80000000E+01 , 


1 .23698000E+00, 
1 . 36631000E+00. 
1 .50058500E+00, 
1 . 63980S00E+00 , 
1 .77902500E+00, 

5.01285714E-02, 
5 . 68429723 E -02 , 
6 . 39 1 135 14E *02 , 
7 . 4660898 9 E *02 , 
8 . 89053233E *02 , 
1 .03163203E-01 , 
1 .2791 1396E-01 , 
1 .52383737E-01 , 
1 . 83630995E -01 , 
2.21 446024E *01 , 

2 58792833E *01 , 

3 1418276 1 E -01 , 
3 69 19 1082E -01 , 
4 . 32534376E -01 , 
5 04125652E-01 , 
5 . 80250000E -01, 
6 . 76860000E-01 , 
7 . 73470000E -01 , 
8 . 79060000E -01 , 
9 . 93630000E -01 , 
1 . 10820000E+00, 
1 .23736000E+00, 
1 .36652000E+00, 
1 50059000E+00, 
1 .63957000E+00, 
1 . 77855000E+00 , 

4 . 95500000E -02 , 
5 . 648947 10E-02 , 
6 . 37693243E -02 , 
7 . 47755056E -02 , 
8 93225679E -02 , 
1 . 03875758E -0 1 , 
1 . 29066269E -01 , 

1 53987351 E -01 , 
1 .85640707E-01 , 
2 . 23825344E -01 , 

2 6 1 554266E -01 , 
3. 1 7297990E -01 , 

3 7258 1234E -01 , 

4 36 150653E -01 , 
5.07906087E-01 , 
5 . 84400000E -01 , 
6 . 8 1230000E -01 , 
7 . 78060000E -01 , 
8 . 83385000E -01 , 
9 . 97205000E -01 , 
1 1 1 102500E+00, 
1 .23967500E+00, 
1 36832500E+00, 
1 . 50203500E+00, 
1 .64080500E+00, 
1 .77957500E+00, 

4 74642857E -02 , 
5 . 44096725E -02 , 
6. 16830180E-02, 
7 . 28820838E -02 , 
8.78365886E-02, 
1 0280606 IE -01 , 
1 .2823 1739E -01 , 
1 .53388510E-01 , 
1 .85226780E-01 , 
2 . 235473 14E -01 , 


4 . 30000000E+01 , 
4 . 5000000QE+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3. OOOOOOOOE+OO, 
5. OOOOOOOOE+OO, 
7 .OOOOOOOOE+OO, 
9 .OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 . 30000000E+01 , 
1 . 50000000E+01 , 
1. 70000000 E +01, 
1 . 90000000E+01 , 
2. 10000000E+01 , 
2 . 30000000E+01 , 
2 . 50000000E+01 , 
2 . 70000000E+01 , 
2 . 90000000E+01 , 
3. 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 . 90000000E+01 , 
4. 10000000E+01 , 
4 . 30000000E+01 , 
4 . 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


1 . OOOOOOOOE+OO , 
3 . OOOOOOOOE+OO , 
5. OOOOOOOOE+OO, 
7 . OOOOOOOOE+OO , 
9 OOOOOOOOE+OO. 
1 . 10000000E+01 , 
1 .30000000E+01 , 
1 . 50000000E+01 , 
1 . 70000000E+01 , 
1 . 90000000E+01 , 
2. 10000000E+01 , 
2 . 30000000E +01, 
2 . 50000000E+01 , 
2 . 70000000E+01 , 
2 . 90000000E+01 , 
3 . 10000000E+01 , 
3 . 30000000E+01 , 
3 . 50000000E+01 , 
3 . 70000000E+01 , 
3 . 90000000E+01 , 
4. 10000000E+01 , 
4. 3 0000000 E +01, 
4 . 50000000E+01 , 
4 . 70000000E+01 , 
4 . 90000000E+01 , 


1 .OOOOOOOOE+OO, 
3. OOOOOOOOE+OO, 
5 .OOOOOOOOE+OO, 
7. OOOOOOOOE+OO, 
9 .OOOOOOOOE+OO, 
1 . 10000000E+01 , 
1 .30000000E+01 , 
1 .50000000E+01 , 
1 . 7 0000000 E +01 , 
1 90000000E+01 , 


1 .30164500E+00, 

-- 

757- 

1 .43097500E+00, 

-- 

758- 

1 .57019500E+00, 

-- 

759- 

1 .70941 500E+00, 

-- 

760- 


-- 

761- 


-- 

762- 

5 . 3433 1 727E -02 , 

-- 

763- 

6.03403329E-02, 

-- 

764- 

6 75449198E-02, 


765- 

8 17812121E-02, 

-- 

766- 

9 60327813E -02 , 

-- 

767- 

1 . 1 5574824E -01 , 

-- 

768- 

1 40179175E-01 , 

-- 

769- 

1 .64530075E-01 , 


770- 

2.02600462E-01 , 

-- 

771- 

2 . 401 74659E -01 , 

-- 

772- 

2 . 86355912E-01 , 

-- 

773- 

3,41 789781 E -01 , 

-- 

774- 

3.96399589E-01 , 

-- 

775- 

4 . 6843886 1 E -01 , 

-- 

776- 

5 41535000E-01 , 

-« 

777- 

6.28555000E-01 , 

-- 

778- 

7 . 25 165000E-01 , 

-- 

779- 

8 2 1 775000E -01 , 

-- 

780- 

9 36345000E -01 , 


781- 

1 . 05091 500E+00, 

-- 

782- 

1 . 1 7278000E+00 , 


783- 

1 .30194000E+00, 


784- 

1 43 1 10000E+00 , 

-- 

785- 

1 .57008000E+00, 

-- 

786- 

1 .70906000E+00, 

-- 

787- 



788- 



789- 

5 . 2969 1 566 E -02 , 


790- 

6 . 00939834E -02 , 

-- 

791- 

6 . 75048 128E -02 , 

-- 

792- 

8 . 2048 1 720E -02 , 

-- 

793- 

9 . 65984878E -02 , 


794- 

1 . 16507664E-01 , 


795- 

1 .41557683E-01 , 

-- 

796- 

1 .66360150E-01 , 

-- 

797- 

2 04793324E -01 , 


798- 

2 4274356 1 E -01 , 

-* 

799- 

2 . 89342649E -01 , 


800- 

3 . 45039520E -01 , 


801- 

3 . 99935934E -01 , 


802- 

4.721 36528E -01 , 

-- 

803- 

5 . 45600000E-01 , 

-- 

804- 

6 . 328 15000E -01 , 

-- 

805- 

7 . 29645000E -01 , 


806- 

8 . 26475000E -01 , 


807- 

9 . 40295000E -01 , 


808- 

1 0541 1500E+00, 

-- 

809- 

1 17535000E+00, 

-- 

810- 

1 3 0400000 E+ 00 , 


811- 

1 43265000E+00 , 


812- 

1 57 142000E+00, 

-- 

813- 

1 .71019000E+00, 

-- 

814- 


-- 

815- 



816- 

5 08882463E -02 , 

-- 

817- 

5 . 80122235E-02 , 

-- 

818- 

6.541 17647E-02 , 


819- 

8.03572239E-02, 


820- 

9 . 53 196760E-02 , 


821- 

1 15555479E-01 , 


822- 

1 . 40840937E-01 , 


823- 

1 . 65879323E -01 , 


824- 

2 . 04447009E -01 , 


825- 

2 . 42534451E-01 , 

-- 

826- 
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2 . OOOOOOOOE+0 1 , 2 . 6 14 14676E -01 , 2 1 0OOOOOOE +0 1 , 2 . 89349726E-01 , 

2 . 20000000E +0 1 , 3 . 1 7344993E -01 , 2 . 30000000E+01 , 3 . 45 120284E -01 , 

2. 40000000 E +01, 3. 7268949 IE -01 , 2 . 50000000E+01 , 4 .0006536 IE -01 , 

2. 60000000 E +0 1 , 4 . 363153 19E -01 , 2 . 70000000E+01 , 4 . 72348038E-01 , 

2 . 80000000E+01 , 5 .08175933E-01 . 2 . 90000000E+01 . 5 . 46 165000E-01 , 

3. 00000000 E +01, 5 . 85025000E -01 , 3 10000000E+01 , 6 . 33440000E -01 , 

3. 20000000 E +01, 6 . 8 1855000E -01 , 3 . 30000000E+01 , 7 . 30270000E -01 , 

3 40000000E+01 , 7 . 78685000E -01 , 3 . 50000000E+01 , 8 . 27 100000E -01 , 

3 . 60000000E+01 , 8 . 83965000E -01 , 3 . 70000000E+01 , 9 . 40830000E -01 , 

3.80000000E+01, 9 . 97695000E -01 , 3 . 90000000E+01 , 1 .05456000E+00, 

4. 00000000 E +0 1 , 1 . 1 1142500E+00, 4 . 10000000E+01 , 1 17568000E+00, 

4 .20000000 E +01 , 1 . 23993500E+00 , 4 . 30000000E+01 , 1 . 30419000E+00, 

4 . 40000000E+01 , 1 .36844500E+00, 4 . 50000000E+01 , 1 43270000E+00, 

4 . 60000000E +0 1 , 1 . 50199500E+00 , 4 . 70000000E+01 , 1 . 57 129000E+00, 

4.8000 OOOOE+OI, 1 . 64058500E+00 , 4 . 90000000E+01 , 1 . 70988000E+00, 

5 . OOOOOOOOE+O 1 , 1 .77917500E+00, 

1 . OOOOOOOOE+O 1 , 

0 , 4.59507143E-02, 1 . OOOOOOOOE+OO , 4 . 950372 16E -02 , 

2 . OOOOOOOOE+OO , 5 . 3 1539295E-02 , 3 . OOOOOOOOE+OO, 5 . 68850416E -02 , 

4. OOOOOOOOE+OO, 6 . 068421 17E -02 , 5 .OOOOOOOOE+OO, 6 . 4541 1765E-02 , 

6. OOOOOOOOE+OO, 7 . 207 1 7058E -02 , 7 . OOOOOOOOE+OO, 7 . 96070381 E -02 , 

8 .OOOOOOOOE+OO, 8 . 7 1465792E -02 , 9 . OOOOOOOOE+OO , 9 . 46898290E -02 , 

1 .OOOOOOOOE+O 1 , 1 .02236364E-01 , 1 10000000E+01 , 1 . 15083294E-01 , 

1 .20000000E+01, 1 .27857208E-01 , 1 . 30000000E+01 , 1 40564 190E -01 , 

1 40000000E+01 , 1 .53209668E-01 , 1 . 50000000E+01 , 1 . 65798496E -01 , 

1 .60000000E+01 , 1 85 17240 1 E -01 , 1 . 70000000E+01 , 2 . 04423723E -01 , 

1 .800000006+01 , 2 . 23559745E -01 , 1 . 90000000E+01 , 2 42587 1 84E -01 , 

2 .OOOOOOOOE+0 1 , 2 61512245E-01 , 2 10000000E+01 , 2 . 89488701 E -01 , 

2 . 20000000E+01 , 3 1 7S20457E -01 . 2 . 30000000E+01 . 3 . 453299 10E -01 , 

2. 40000000 E +01, 3 7293 1099E -01 , 2 . 50000000E+01 , 4 . 00336907E -01 , 

2 . 60000000E+01 , 4 . 36587057E -01 , 2 . 70000000E+01 , 4 72622878E -01 , 

2 . 80000000E+01 , 5 . 084566 1 7E -01 , 2 . 90000000 E +01 , 5 . 46455000E -01 , 

3. 00000000 E +01, 5.85325000E-01, 3 . 1 OOOOOOOE +0 1 , 6 . 33755000E -01 , 

3 . 20000000E+01 , 6 . 82 1 85000E -01 , 3 . 30000000E+01 , 7 .30615000E-01 , 

3 . 40000000E+01 , 7.79045000E-01 , 3 . 50000000E+01 , 8 .27475000E-01 , 

3 . 60000000E+01 , 8 . 84350000E -01 , 3 . 70000000E+01 , 9 41225000E-01 , 

3 . 80000000E+01 , 9 . 98 100000E -01 , 3 . 90000000E+01 , 1 . 05497500E+00 , 

4 . OOOOOOOOE+O 1 , 1 1 1185000E+00, 4 . 10000000E+01 , 1 17599000E+00, 

4. 20000000 E +01, 1 .24013000E+00, 4 . 30000000E+01 , 1 . 30427000E+00 , 

4. 40000000 E +0 1 , 1 3684 1000E+00 , 4 . 50000000E+01 , 1 43255000E+00 , 

4 . 60000000E+01 , 1 . 501 82000E+00 , 4 . 70000000E+01 , 1 . 57 109000E+00 , 

4 . 80000000E+01 , 1 .64036000E+00, 4 . 90000000E+01 , 1 70963000E+00 , 

5 . 00000000 E +0 1 , 1.778 90000 E +00 , 

ENDPHS= 1 $ 

LSGENDAT 

EVENT=50,CRITR=5HGCLAT,ENDPHS=1 ,ETAPC( 1 ) = 1 , NENG= 1 ,NPC( 2 1 ) = 1 , 

ALPPC( 1 >=180. , 

NPC( 2 ) - 1 ,DT = 1 , 

$ 

LSGENDAT 

EVENT = 55 ,NPC( 2 ) =4 ,DT = 100 ■ CRITR=6HALTP .VALUE =28 . 8348843 , ETAPC( 1 ) =0 , ENDPHS' 1 $ 
LSGENDAT 

E VENT= 100 ,CRITR = 6HALTIT0 .VALUE = 400000 , ENDPH5 = 1 , PRNC= 1 ,NPC( 12 )* 1 , 

P I NC = 1 00 , 

IGUID=0, 1 ,0,IGUID(6) = 1 ,0, 1 ,NPC( 1 5 ) = 1 ,DT= 1 ,NPC( 2 ) = 1 ,NPC( 5 ) =4 , 

ALPPCt 1 ) =37 .00 10556, BNKPC( 1 > = -69.2957708, 

NPC(8)=-2,NPC(5)=4,DTIMR=1 ,8ETA=0. , 

$ 

LSGENDAT 

EVENT= 105, CRITR=6HHEATRT,PINC=20, VALUE = 1 15, IGUID(6) = 1 ,0,4, 

NPC(9)=1 ,NPC(22)=1,PRNC=1 , 

NDEPVS=1 ,DEPVRS=5HVELAD,DEPTLS= .00001 ,INDVRS=6HETAPC1 , 

US=. 1 ,MAXITS=10, 

DGF( 3 ) =6HHEATRT , IDGF(3 ) = 1 ,KDG( 3 ) = 30,KRG(3 >=2000, 

$ 

LSTBLMLT $ 

LSTAB TABLE=5HGDF3T , 1 .5HTDURP , 5 , 1 , 1 , 1 , 

0,115,15,120.30,124,60,125,1000,125 $ 

LSTAB T ABLE=6HGNMN3T ,0, - 180S 
LSTAB TABLE=6HGNMX3T,0,0$ 


-- 827-- 
-- 828-- 
-- 829-- 
-- 830- - 
-- 831 -- 
-- 832-- 
-- 833-- 
-- 834-- 
-- 835-- 
-- 836-- 
-- 837-- 
-- 838-- 
-- 839-- 
-- 840“ - 
-- 84 1 - - 
--- 842-- 
-- 843-- 
-- 844-- 
-- 845-- 
-- 846-- 
-- 847-- 
-- 848-- 
-- 849-- 
-- 850- - 
-- 85 1 - - 
-- 852-- 
-- 853-- 
-- 854-- 
-- 855-- 
-- 856-- 
-- 857-- 
-- 858-- 
-- 859-- 
-- 860- - 
-- 86 1 - - 
-- 862-- 
-- 863-- 
-- 864-- 
-- 865-- 
-- 866 -- 
-- 867-- 
-- 868 -- 
-- 869-- 
-- 870- - 
-- 870- - 
-- 87 1 - - 
-- 872-- 
-- 873-- 
-- 874-- 
-- 874-- 
-- 875-- 
-- 875-- 
-- 876-- 
-- 877-- 
-- 878-- 
-- 879-- 
-- 880- - 
-- 88 1 - - 
-- 88 1 - - 
-- 882-- 
-- 883-- 
-- 884-- 
-- 885-- 
-- 886 -- 
-- 887-- 
-- 887-- 
-- 887-- 
-- 888 -- 
-- 888 -- 
-- 888 -- 
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ORlGfNAL PAGE 13 


L$TAB TABLE=6HGN0M3T, 0,-70. OOOOOOO. ENDPHS=1 $ -- 888- 

L$GENDAT -- 888- 

EVENT=180,CRITR=6HINC , ENDPHS= 1 ,NDEPVS=0 , -- 889- 

IGUI0(6)=3,O,3,DBANK=O,DALPHA=12 $ -- 890- 

LSGENOAT -- 890- 

EVENT = 190 ,CRITR=5HTDURP , PINC= 1000 .VALUE* 10 , ETAPC( 1 ) = 1 , ENDPHS 5 1 $ -- 891- 

LSGENOAT -- 891- 

£ VENT =200 , CRITR = 6HALTA , ENDPHS = 1 ,DT = 10 . IGUID ( 6 ) = 1 ,0 , 1 , -- 892- 

£TAPC( 1 )=0, -- 893- 

PRNC= 10, -- 894- 

$ -- 895- 

LSGENDAT -- 895- 

EVENT=210, CRITR=6HALTIT0 .VALUE = 400000 . NPC( 2 ) =4 ,DT = 100 ,NPC ( 5 ) =0 ,NPC ( 8 ) s 0 , -- 896- 

ENDPHS* 1 $ -* N 897- 

LSGENDAT --' 897- 

EVENT =220 , CRITR=6HGAMMAI , VALUE =0 ,NPC( 2 )= 1 , DT= 1 , ETAPC( 1) = 1 , -- 898- 

IGUID=3 ,0, 1 ,ALPPC( 1 )=0,BETPC( 1 )=0 ,BNKPC( 1 )=0, ENDPHS = 1 $ -- 899- 

LSGENDAT -- 899- 

E VENT = 1000 ,CRITR=6HVELI , VALUE =25354 . 1388 , ENDPHS= 1 , ENDPRB= 1 , ENDJOB= 1 $ -- 900- 
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A. 3.1 INPUT FILE DESCRIPTION 


A. 3. 1.1 INPUT CONDITIONS 


The following is a list of the input conditions for the 
plane change trajectory. 


Initial inclination = 28.5 degrees 


Initial Orbit altitude = 160 nautical mile circular orbit 


True anomaly (Truan) =150 degrees 


Vehicle Model: 


Reference Surface Area (SREF) = 177.4 square feet 


Reference length (LREF) = 25 feet 


Initial weight (WGTSG) = 11,430 pounds 


nom i na 1 


Engine ISP (ISPV) = 295.0 seconds 


Coefficients of lift (CL) and drag (CD) as a function of angle- 
of-attack and a correlation parameter of either altitude (ALTI- 
TO) , viscous interaction (VINV) , or mach number (MACH). 


Spherical, rotating earth model. 


1976 standard atmospheric model 


A. 3. 1.2 TRAJECTORY PHASES 


The trajectory has 11 different phases called events and are 
bered: 1 , 50, 55, 1 00, 1 05, 180, 1 90,200,210,220, and 1000. Each 

represents a change in either the vehicle and/or the environmental 
ditions. The following is a list of the events and their purpose, 
of these phases are shown on Figure A-4. 

Event Humber Purpose 


1 

50 

55 

100 

105 

180 


Sets initial conditions 
Starts deorbit burn at best point 
Stops deorbit burn at best point 
Turns on atmosphere at 400,000 feet 
Starts level-off based on heat-rate 
Starts exit based on best inclination 


num- 

phase 

con- 

Most 
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190 


Completes the exit set-up by going 100% throttle 
200 Cuts-off engines based on desired apogee altitude 

210 Turns off atmosphere at 400,00 feet 

220 Starts circularization burn 

1000 Stops circularization burn when a 160 nautical mile 

circular orbit is obtained. 


A. 3.1 


POST has been given a number of parameters which it can vary to 
optimize the problem. There are 12 variables that are what we call 
exterior parameters and two interior parameters. These variables are 
listed below with the associated phase/event numbers and a brief 
descr ipt ion. 
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Exter i or 
Independent 
Var i abl es 
Names 

Event 

Number 

Description 

GCLAT 

50 

Geocentric latitude for deorbit burn 

ALTP 

55 

Perigee altitude to stop burn 

ALPPC1 

100 

Entry angl e-of-attack 

BANKPC1 

100 

Entry bank angle 

ALPPC1 

105 

Level-off angle-of-attack 

DEPVLSI 

105 

Cruise acceleration 

INC 

180 

Exit i nc 1 i nat i on 

OBANK 

180 

Ex i t bank ang 1 e 

DALPHA 

180 

Exit angle-of-attack 

ALTA 

200 

Apogee altitude for engine cut-off 

ALPPC1 

200 

Angle-of-attack after engine cut-off 

BANKPC1 

200 

Bank angle after enging cut-off 

Interior 
Independent 
Var i abl es 

Event 


Names 

Number 

Description 

ETAPC1 

105 

Throttle % to maintain cruise accelei 

BANK 

105 

Bank angle to maintain heating rate 


A. 3. 2 OUTPUT DESCRIPTION 


As might be expected, POST has the capacity to generate a large 
amount of output data and plot information. Included below are samples 
of the output data and generated plots from our sample trajectory. 

A. 3. 2.1 OUTPUT DATA 

Figures A-1 thru 3 present a sample from the output data file for 
the nominal 20 degree plane change. Figure A-1 presents the 270 parame- 
ters offered as nominal output. Specifically, the parameters are for 
the start of phase 105 (cruise phase). Figure A-2 presents the same 
parameter for the end of the trajectory (last frame of phase 220 and 
first and only frame of phase 1000). Figure A-3 provides a summary of 
the trajectory to include the independent variable values (u) , the end 
condition errors (e) , the optimization parameter value (PI) and a tar- 
geting parameter (P2) . 

A. 3. 2. 2 TRAJECTORY PLOTS 

Shown in Figures A-4 thru 12 are the baseline trajectory plots gen- 
erated by POST for a heat-rate of 80 and 125 BTU/ft 2 -sec. The plots for 
QD0T = 125 were generated with the trajectory generation input deck 
shown in Section A. 3. The Figures present a time history of altitude 
above the planet, relative atmospheric vehicle velocity, angle-of-at- 
tack, bank angle, Engine throttle percentage, heat-rate, heat load, and 


210 


ERV CRUISING SYNERGISTIC PLANE CHANGE 0-D0T=125 PROBLEM NO 



cs O 

■*- 

9 

O CO 


o •*- 

CO 

CO 

CM 

CO 

CM 

T- 




CM P- 

r- rf 



CO 



o o o 

o 

o o o o o o o 

o o o o o 



o o o o 



O 



4 + 

4 

4 

4 

4 

1 

4 4 

4 

4 

I 

1 

i 

4 

i 



4 4 

4 4 



4 



UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

a 

a 

a 

a 

a 

a 



a a 

a a 



a 



r~ co 

CO 

p- 

p** 


t- 

o o 

■»— 


CO 

CO 

o 

ID 

CO 



cm in 

P- CD 



r- 



p- 


rr 

T- 

ID 


O O ID 

ID 

O’ 


in 

in 

CO 



CO ID 

r- O 



CP 



CM CM <D 

Or 

rr 

n 


O O -O 

•O’ 

co 

CM 

ID 

03 

o 



rr co 

in id 



CM 



<J> CD 


a> 

CO 

n 

P- 

O O co 

CO 

O ID 

in 

ID 

ID 



O CM 

rr rr 



CD 



■*- co 

CM CO O 

CM 

o 

O O f- 

p- 

CO 

CO 

p- 

CO 

CM 



CO CM 

O - 






CT> ID 

O rr 

in 

ID 

o 

o 9 

CD 

CD 

rr 

CO 

O CD 

P~ 



r- CO 

in O 



ID 



O rp CM 

in 

CD 

ID 

ID 

o o 

P- 

P- 

CM 

in 

O 


CO 



CM CD 

ID 00 



CD 



■*" CM 

f- 

CM 

in 


** 

O LO 

in 

in 

CM 

ID 

OOP- 



co in 

in in 






rr •»- 

CM 

CO 

co 

CO 


O CM 



CO 

p- 

ID 


p» 



- - 

CO CM 



CM 







CM 


3 ►- 

<x 

p- 








o 







t— < 




a 


< o 

X 

o 

o 

a 

O 


«— * 



< 

a cj 






Z CD 




z 

o 

> 

a 

b - 

00 

co 

CO 

»-4 

< 



Q- OC 

m oc 






uj Z 

t~-t 


t-4 

a. 

z 

o a 

— i 

—1 

— 1 

b- 

3 

PM 

b~ 



C3 a 

z 






h- o 

X 

> 

PM 

a. 

tn 

Z -1 

<x 

< 

o 

•— « 

<1 

CO 

a 



OC O 

•-« CJ 



a 



<L a 

< 

< 

< 

a 

< 

*-• < 

o 

O 

oc 

a 

>- 

<t 

CO 



< a. 

*- > 



o 



r- <N 

CO 

rr 

rr 

co 

CM 

o o 

o 

o 

V- 

o 


O 

o 




r» CO 



CM 



o o 

O 

O 

O 

O 

o 

5 o 

o 

o 

o 

o 

o 

O 

o 



O r- 

o o 



O 



t 4 

4 

4- 

4 

4 

4 

4 4 

4 

4 

4 

4 

4 

4 

4 



4 4 

4 4 



4 



UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

a 

a 

a 

a 

a 

a 



a a 

a a 



a 



TT CO 

p- 

P- 

CO 

p» 

O 

o o 

O 

O 


r** 

o 

CO 

o 



co — 

ID CO 



rr 



r~ IN 

CO 

CO 


ID 

co 

O Q 

o 

o 

CO 

O’ 

CO 

o 

c > 



O o- 

CM CD 



O 



co 01 

to 

CD 

O 

■O’ 

CO 

O O 

o 

o 

o 

CD 

ID 

in 

() 



O O’ 

co p- 



CM 



CM CO 

rr 

r- 

CD 

■o- 

■*— 

o o 

o 

o 

CO 

ID 

o 

r-~ 

o 



mOO-- 



ID 



*- in 

O 

CM 

CM 

CM 

in 

o o 

o 

o 

CD 

O 


i" 

o 



in -r* 

in a* 



rr 



CM O 

'T 

r- 

in 

ID 

p- 

o o 

9 

o 

P~ 

P~ 

CM 

O 

o 



in in 

O O 






in co 

CM 

ID 

co 

ID 


o o 

o 

o 

CD 

o 

in 

r- 

o 



p- CO 

r- O 



ID 



p~ rr 

CO 

CD 

co 

r- 

CO 

o o 

o 

o 

CD 

co 

■*“ 

in 

o 



ID O 

O O 



CO 



r- co 


- 

- 

CO 

CM 

6 b 

o 

o 

CM 

r- 

p~ 

r- 

- 



id in 

CO CM 



CM 









X 1- 

►— 

►- 








o z 




a 







o 

a 

o o 

o 

o 





O’ 



o o 

Cl. oc 



<t 



in (J 




z 

o 

a q 

Q 

Q 

OC 

oc 

oc 

►— « 

X 



to SE 

b- t-4 



a 



uj Z 

*— < 

t—t 

•— • 

cr 

OC 

a a 

»— 

y 

3t 

b- 

a 

>- 

< 



oc a 

Z CJ 



a 



OC O 

X 

>- 

PM 

a 

a 

Z -1 

UJ 

z 

<1 

*—* 

o 

CO 

z 



a z 

t-H > 






Q- a 

> 

> 

> 

a 

3t 

»-* < 

00 

00 

>• 

a 

oc 

<. 

X 



Q. <X 

b- o 



> 



p- o 

p- 

in 

ID 

CM 

CM 

o — 

in 


CM 



O 

o 



r- CM 

O O’ 



CM 



o o 

O 

O 

o 

o 

O 

o o 

■»— 

o 

O 

O 

o 

O 

o o 


o o o o 



o 



• 4 

4 

4 

4 

4 

4 

4 4 

1 

4 

4 

4 

4 

4 

4 

4 


4 4 

4 4 



4 



UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

a 

a 

a 

a 

a 

a 

a 


a a 

a a 



a 



■*- rr 

id 


T— 

ID 

O 

o o 

CD 

CD 

ID 

rr 

co 

CD 

O co 


co in 

— co 



c— 



O in 

to 

rr 

ID 

CM 

co 

o o 

p- 

O 

CD 

p^ 

in 

CD 

O CD 


ID co 

CO ID 



CO 



CO CM 

ai 

m 

CO 

■»— 

co 

o o 

▼— 


CO 

ID 


in 

O O 


O CD 

ID CO 



o 



in 

to 

CD 

CD 

o 


o o 

CO 

o 

in 

P^ 

CO 

■*— 

O O 


O O ID ID 



CO 



cr> 

CM 

co 

CO 


in 

9 o 

r* 

o 

CD 

CO 

O 

in 

O ID 


cd r- 

00 ID 



CM 



co •»- 

■*— 

CO 

CO 

ID 

r- 

o o 

O 

o 

rr 

CM 

in 

CM 

O ID 


CD ID 

O O 



O 



co 


<D 

r- 

CM 


o o 


o 

ID 

in 

•*— 

co 

O CM 


CM ID 

P- 00 



CD 



rr <d 


in 

O 

ID 

CD 

O in 

o 

o 

P~ 

O 

ID 

CM 

O rr 


o- in 

r- CO 



•O’ 



CM CM 

CM 




CM 

O CM 


p** 


in 

in 

■o- 

o 

CO 


CO CO 

in cm 


■# 

CD 







i 





• 






# 




* 



* 






z 



C.3 







* 




-» 



* 





a 

o 

— i < 


z 





O’ 


* 


a 






in < 




z 

o 

3 X 

<t 

< 

*—« 

»-* 

►— < 

»— < 

z 




> a 


in 



o 

Z -j 




a 

t— • 

z a 

b- 


a 



X 

►— < 

o 


U Z 

CD > 


a 

OC 


o 

UJ U 

fM 

♦— » 

i— i 

a 

UJ 

>- -l 

UJ 

z 

o 

<i 


CO 

z 

z 

t— < 

Z -x 

OC CL 


in 

a 


o 

o id 

X 

>- 

rsi 

u 

* 

CH <. 

00 

00 

a 

>• 

a 

< 

X 


cq 

to a 

< <c 


in CD 


in 

88 

V- 

f- 

■r- 

co 

o 

o o 

o 


co 

CO 

o 

in 


CM 

o 

r- CM 


IN 

a 



O 

O 

o 

o 

o 

o 

o o 

o 

o 

O 

O 

O 

o 

O 

O 


o o o o o 


O 



4 4 

4 

4 

4 

4 

4 

4 4 

4 

4 

4 

4 

4 

4 

4 

i 

o 

• 4 

4 4 

4 

>- 




UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

a 

a 

a 

a 

a 

a 

a 


a ul 

ui a 

a 

K- 

a 



o rr 

rr 

r- 

r- 

p- 


O CD 

o 


in 

in 

in 


ID 

P- 

*- 

co o 

rr o* 

O 

to 

O 


UJ 

o in 

co 

CO 

co 

O 

in 

o 

o 

ao 

CM 

CD 

co 

co 

O’ 

in 

O- 

CO CO 

CO CO 

O 

U 

in 


in 

O cm 

P- 

CM 

CM 


in 

O O’ 

o 

O 

CM 

co 

o 

CO 

CO 

in 

to 

ID O’ 

in in 


O 

CO 


< 

o rr 

10 

<D 

CD 


CM 

o o 

o 

00 



CD 


CO 

r- 

a 

CO O 

IN a» 

CO 

a 

co 


X 

O — 

cd 

O 

o 

” 0 - 


o o 

o 


CO 

b- 

CD 

P- 

ID 

CO 

a 

O’ ID 

P- ID 

p** 

a 

CD 


CL 

o •*- 

o 


■»- 

CO 

CD 

O CD 

o 

in 

v- 

ID 

O 

ID 

CO 

ID 

a 

o CD 

p- O 

> 

CM 



o - 

CO 

ID 

ID 

ID 

CM 

o CD 

o 

CM 

O 

O 

CO 


•r- 

ID 


00 CM 

CO r- 

in 


O 


■* 

o <d 

in 

CO 

CO 

r- 

CM 

O CO 

o 

CM 

CO 

ID 

in 

in 

CO 

CM 

* 

rr Q 0 

P- CM 

o 

* 

CO 


* 
















* 







* 

O CM 

in 

in 

in 

co 

CO 

o in 

o 

*•" 

CM 

CO 

in 

•O’ 

ID 

CM 

■* 

■— r- 

r- CO 

co 

* 

rr 




H 

OC 

<x 

o 


— 1 





b~ 












a. >- 

— J 

— J 

— i 

z 


-J 





C 

o 

CO 



z z 

z a 

Z 





a: < 

UJ 

UJ 

UJ 

OC 

t— 

3 



O 

b- 

a 

z 

X 

> 


a <t 

O a 

o 





X -J 

> 

> 

> 

z 

o 

Z DO 

CO 

CO 


a 

X 

> 

<t 

z 


U 3 

a > 

a 


a 



a a 

PM 

PM 

PM 


o 

a x 

>- 

PM 

OC 


a 

a 

z 



CJ oc 

CD CD 

CL 


b- 



H- CD 

<t 

< 

< 

o 


*-• < 

< 

< 

o 

a 

b- 

cc 

X 

> 


a H- 

0. CL 

< 


< 



O r- 

-r> 


<«- 

CM 


■*“ CM 

o 

CO 

f- 


CM 


in 

CM 


r- P«* 

O CM 

CM 





O O 

o 

O 

o 

O 

o 

o o 

o 

O 

o 

O 

O 

O 

O 

o 


o o 

O O 

o 


o o 


4 4 

• 

• 

1 

t 

4 

» 4 

4 

4 

• 

• 

4 

4 

4 

4 


4 4 

4 4 

4 


4 

4 


UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

a 

a 

a 

a 

a 

a 

a 


a a 

a a 

a 


a 

a 


O co 

O 

CO 

CO 

CD 

ID 

CO CD 

o 

co 

O’ 

CO 

O 

ID 

CM 

o 


in o- 

id co 

ID 


CO 

co 


O O 

in 

-T 

’sr 

CM 

OO 

ID ID 

o 

CD 

CO 

ID 

O 

O 


O 


O CO 

ID rr 

CD 


CD 

CD 


O co 

co 


Tf 

CM 

*T 

P- O 

o 

CM 

p~ 

ID 

o 

CO 

r- 

CM 


O’ rr 

ID P- 

r- 


CO 

CD 


O CD 

O 

in 

in 

r- 

CM 

O o 

o 

CD 

O’ 

O’ 

o 

CM 

CD 

ID 


co •*- 

<D r* 

p~ 


ID 

in 


O "* 

rr 

t- 


ID 

in 

ID ID 

o 

O 

ID 


o 

CO 

CO 

O 


00 ID 

O CO 



r- 



O rr 

cr> 

<D 

CD 

in 

•O’ 

ID — 

o 

O’ 

CM 

o 

o 

CD 

CM 

O’ 


O ■’*- 

00 CD 

■»— 


p- 

CD 


O - 

O 

CD 

(D 

CD 

T— 

in O 

o 

O’ 

in 

CM 

in 

ID 

ID 

CD 


CD O 

CM CO 

o 


CO 

P- 


O -• 

CM 

T 


CO 


CO ID 

o 

CM 

O 

CM 


CM 


CO 


ID CO 

cm ao 

rr 


p- 

ID 


O CM 

in 

in 

in 

CM 

- 

t- m 

o 

O’ 

CM 

CO 

- 

CM 

CM 

CM 


CO r- 


CO 


rr 

00 



t-H 

oc 

< 


h- 






b- 



a 


X 

b~ 






m a 

<t 

-s 

< 

a 

X 






OC 


CM 

< 


<x 

< z 

z 





UJ <x 

Z 

z 

z 

< 

o 

-i co 

CO 

CO 



b- 

X 

z 

a 


< 3 

a <c 

o 


oc 



S oc 

Z 

z 

z 

> 

l-U 

< x 

>- 

PM 



<1 

o 

►— t 

O 


I- Z 

CJ <1 

a 


a 

t-H 


CJ 

<1 

< 

< 

PM 

UJ 

1 - < 

<t 

< 

o 

a 

a 

< 

z 

i— ♦ 


a a 

CD a 

Z 


> 

a 


*- CD 

o 

cd 

o 

< 


UJ LU 

a 

a 

CJ 

U 

X 

z 

X 

> 


< tn 

a. z 

<t 



CD 


co m 

rr 

rr 

■o- 


CO 

CO 

O 

O 

CM 

w- 

o 

•«> 

CM 

CM 

CO 

CM CO 

P» CM 

CM 

CO 

CO 

CM 


O o 

o 

o 

o 

O 

O 

o o 

O 

O 

o 

o 

o 

o 

O 

o 

o o o 

o o 

O 

O o o 


+ + 

4 

4 

4 

i 

4 

• 4 

4 

4 

1 

• 

4 

4 

4 

4 

4 

4 4 

4 4 

4 

4 

4 

4 


UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a a 

a a 

a 

a 

a 

a 


CM CM 

to 

O 

O 

ID 

■o 

CO O’ 

O 

O 

ID 

O 

O 

CD 

O 

CM 

CM 

CD co 

rr in 

in 

CM 

co 

CO 


in p- 

in 

CO 

CO 

00 

ID 

ID ID 

O 

O 

00 

o 

O 

O’ 

O 

ID 

in 

CM CO 

O ID 

co 

in 

CD O 


cm r- 

o* 

rr 

■o- 

O 

O’ 

P^ O' 

o 

o 

CO 

p- 

O 

CO 

o 

CD 

CM 

ID CO 

O in 

CD 

CM O r r 


— a 

CM 

•»— 

■*— 


O’ 

o o- 

o 

o 

CD 

in 

o 

CO 

o 

■*— 


■»“ CD 

O CM 

O 


in 

00 


in co 

CO 

to 

ID 

CO 

ID 

ID ID 

o 

o 

O 

O 

o 

CD 

o 

CD 

in id p** 

O f- 

p- 

in 

rr 

CD 


to CM 

r- 

in 

in 

•o 

ID 

ID ID 

o 

o 

p- 

VD 

o 

CO 

o 

ID 

ID 

CD — 

in o 

CD 

ID 

CM 

co 


m ci 

CO 

in 

in 

CO 

CM 

in cm 

o 

o 

CD 

co 

o 


in 

CO 

in 

CM CM 

•*" CO 

CD 

in 


co 



CO 

CM 

CM 

CM 

O’ 

CO O’ 

o 

o 

CD 

p^ 

o 

CO 


ID 


ID CO 

r* - 00 

in 


•*" 

p- 


CO CM 

CM 

CM 

CM 

CM 

CM 

P" CM 

o 

o 

•o 

CO 

o 

ID 

- 

CD 

CO 

00 CO 

CM r- 

co 

CO 

CM 

rr 


o 





►— 











> 

Q 






b- 




a 

in 









a 


o 

< z 

a 





UJ t-» 

►h 

at 

< 

< 

D 

CO 

00 

03 




a 

X 

oc 

a 

CL OC 

OC <1 

> 

a 


HU 


Z b- 

— J 

-j 

-J 

z 

oc 

< X 

>- 

PM 




z 

< 

Z 

Z 

b- a 

O CJ 

Z 

Z 

OC 

a 


t-» _i 

UJ 

UJ 

UJ 

< 

X 

P- b- 

*— 

b~ 

<1 

z 

>- 

>- 

Z 


to 

a z 

CL CJ 

< 

t— i 

a 

> 


t— <t 

> 

> 

> 

o 

h~ 

UJ u. 

a 

a 

U 

u 

o 

Q 

X 

b~ 


< a 

< a 

a 

b - 

Q 



211 


igure A-l. Sample of POST output data for cruise phase 
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inclination for the nominal conditions. The time reference for the 
plots (t=0) is atmospheric interface at 400,000 feet. The phase/event 
divisions and numbers are shown on the altitude time history. 
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Figure A-4. 8aseline altitude profiles 
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Figure A-5. Baseline vehicle velocity relative to the atmosphere 
profi les. 
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Figure A-6. Baseline angl e-of-attack profiles. 
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Figure A-7. Baseline bank angle profiles 
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Figure A-10. Baseline heat-load profiles. 
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APPENDIX B 


DRAPER ERV COMPUTER SIMULATOR 
AND GUIDANCE ALGORITHM 

B.1 SIMULATOR DESCRIPTION 

The ERV simulator is a three degree-of-f reedom (point mass) simula- 
tor used to simulate the atmospheric phase of an orbital plane change 
maneuver. A fourth order Runge-kutta integrator is used to integrate 
the vehicle state. A spherical earth gravity model and a 1962 standard 
atmospheric model are used. Bank angle, angle-of-attack, and throttle 
are controllable with bank angle and angle-of-attack being rate limited. 
Thrust is modeled as a force along the ERV's longitudinal axis through 
the center of gravity. Complete ERV CL and CD data including viscous 
interaction effects are used. 

B .2 SIMULATOR AND GUIDANCE ALGORITHM COMPUTER PROGRAM 

The computer program for the CSDL ERV simulation and guidance algo- 
rithm actually consists of a control program and a plot program that 
collects and uses 14 other programs to generate the trajectory data and 
plot information. While some of the programs were developed or modified 
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by the author for his own application, most of the code was modified and 
or developed by R. Richards based on previous in-house computer simula- 
tors at C.S Draper Laboratory (19) . 

The software is written in MAC, a language developed at and used 
only at Draper Laboratory. It is an excellent language for dealing with 
vector and matrix quantities and operations. However, because of the 
unique language in which the software is written and because the total 
collection of code is (nearly 300 pages, the actual code is not included 
(The computer code can be obtained though by contacting the CSDL NASA 
Department, Guidance and Navigation Division, 10E). 

To provide at least a basic understanding of the ERV computer simu- 
lator and guidance program, a listing and brief description of the pro- 
grams used are shown below. Figure B-1 provides a top level flowchart 
showing how the programs are grouped and connected. 

1. Control Program : Provides a top-level input/output control channel 
for the various subprograms. Allows for easy changes in vehicle and 
atmospheric conditions (e.g., density, aerodynamic, and engine perform- 
ance) by the operator. 

2. Plot Program : Uses data collected by the control program to plot 
various trajectory parameters (e.g., altitude, velocity, and heat-rate). 
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3. Driver Program : Reads data from the control program and calculates 
the initial vehicle state to be used by the other subprograms. Also 
collects information from the other programs and prints the desired out- 
put. 

4. ERV Simulator Program : Uses the initial ERV state information along 
with atmospheric, vehicle performance, and guidance information to esti- 
mate the current vehicle state. 

5. Guidance Program : Provides the guidance information produced by the 
guidance algorithm described in Chapter 5m Outputs include the commanded 
angl e-of-attack , bank angle, and throttle setting. 

6. Pred i ctor /Corrector Program : Predictor/corrector algorithm for the 
guidance program. Outputs include the estimated final inclination 
change, apogee altitude, and final vehicle weight. 

7. Conic Program : Performs the conic calculations for the 
predictor/corrector and simulator program. 

8. Kepler Program : Solves the Kepler problem for the conic program. 


226 


9. Atmosphere Program; Provides the atmospheric model data (e.g., den- 
sity, temperature, and dynamic pressure) for the various subprograms. 


10. CLCD Driver Program: Provides the ERV aerodynamic data to the sim- 

ulator program. Interpolates the vehicle's CL/CD data based on angle- 
of-attack and a correlation parameter of either altitude, viscous, or 
mach number . 

11. VINV Programs : Calculates the viscous correlation parameter values 

for the CLCD Driver program. 

12. DATA. DATA 1 .2 and 3: Contains the tables of ERV aerodynamic data. 
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Figure B-l. ERV top-level computer simulation and guidance algorithm flowchart. 


















